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INTRODUCTION, 



T am of the opinion that the present work, entitled ''Tlie Harmonic Formula of 
Fourier and Bessel and Its Application to the Study of the Diurnal Variation of the At- 
mospheric Pressure in Manila during the Period 1890-1909/' prepared by Rev. Antonio 
Galán, S. J., a temporary member of the staff of the Manila Observatory, will be very 
well received by those meteorologists who think that such studies and discussions should 
be made in all the principal observatories where sufficient material has been collected 
for a number of years. 

I hope that in a short time we may be able to find time and opportunity to work 
up, in a manner similar to that done by Father (Jalan for the barometric data, the 
observations collected in the Observatory on the temperature, humidity, and other pe- 
riodic elements of meteorology. 

And I hope also that the stay in Manila of leather Galán, who is appointed to con- 
tinue the splendid work of Fathers Vines, Gangoiti, and Gutierrez Lanza of the Belen 
Observatory, Habana, may be as useful to him as my stay in Hal)ana during the first 
few months of 1892 was to me. 

José Algué, 
Divdctor, Weather lifneafL 

Manila Observatory, J ¡me 19, I9l/^. 



THE HARMONIC FORMULA OF FOURIER AND BESSEL 



The harmonic formula of Fourier and Bessel is usually presented in one or other 
of the two following forms: 

(I) ?/=:M+Pi cos x+q, sin x-^p^ cos 2x + q2 sin 2x-^p^ cos Sx+q,, sin Sx-}- .- 

(II) y = M+u, sin(U,+x) +u^ ^in{U,+2x) -\~u, sin(i7,+3a;) + 

It is evident that these two formulae are one and the same under different form, for 
developing (II) we have: 

y = M-\rM, sin U, cos x+u, cos U^ sin x-{-u.^ sin C/, cos 2x+u., cos U. sin 2x~\-. .....etc. 

and making in the expression — 

ui sin Ui=^pi ui cos Ui=qi 

iC2 sin Í72=p2 ti2 cos U2=q2 

etc. etc. 

we find that (II) is reduced to (I). 

The practical applications of this formula are many and important. 

Fourier's theorem (say Thomson and Jait) is not only one of the most beautiful results of modern 
analysis, but may be said to furnish an indispensable instrument in the treatment of nearly every 
recondite question in modern physics. To mention only sonorous vibrations, the propagation of electric 
signals along the telegraph wire, and the conduction of heat by the earth's crust, as subjects in their 
generality intractable without it, is to give but a feeble idea of its importance. 

In the present study we limit ourselves to the application which is made of the for- 
mula in meteorology ; and we will endeavor to show what each one of its elements repre- 
sents, and point out its utility and practical application.' 

The greater number of the phenomena studied in meteorology are periodic, and 
the curves which represent them are consequently periodic curves — ^that is to say, 
curves which represent the movement of a point which passes an indefinite number 
of times over the same path in the same time. 

Fourier proved that any periodic curve could be resolved into a series of harmonic 
curves of periods 1, |, i, etc., of the given curve, and that only one combination of these 
elementary curves was possible to reproduce a specified curve. This corresponds with 
the fact observed by Helmholtz that the same composite sound is always resolved into 
the same elementary sounds. 

Prescinding for the moment from the practical utility which this resolution of 
curves may have in meteorology, let us examine how the formula which we are study- 
ing is nothing more than the analytical expression of periodic curves in terms of their 
harmonic components. 

A harmonic curve is one which represents graphically the harmonic motion of a 
point. 



^ One of the best practical treatises which we know on the employment of Fourier's series in mathe- 
matical physics is "An Elementary Treatise on Fourier's Series and Spherical Harmonics" by W. E. 
Byerly, published by Ginn & Co., New York. 

The original work of Fourier himself **Théorie Analytique de la Chaleur" is perhaps as up to 
date as any other modern publication on this subject. Mellor*s "Higher Mathematics for Students of 
Chemistry and Physics," Longmans, London and New York, may also be consulted. 
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Harmonic motion is that to which the projection of a 
point, moving uniformly . around a circumference on to a 
diameter, is subjected. The circumference is called the 
circle of reference. 

If the point M (fig. 1) moves around the circumference 
uniformly, the motion of its pro- 
jection P on the diameter BD is 
harmonic; while M describes the 
circumference once in the direc- 
tion ABCDA.Vae point P traverses 
the diameter in the direction 
OBODO, 

As a help to those who may 
not be familiar with this class of motion, the following defini- 
tions of terms are given : 

Period is the time required for the point M to describe 
the circumference, or for P in its double journey of the 
diameter. 

Elongation is the variable distance between the center 
and the moving point P. 

Amplitude is the maximum elongation, or maximum 
distance which the point P can be separated from O; as is 
evident, this distance is the radius of the circle of reference. 
and A are generally taken as the starting points of the 
motion. The zero of time is arbitrary; it may be taken on 
any point in the circumference. 

Phase is the angular distance between the moving point 
and the starting place of the motion. In speaking of the 
point which moves on the diameter, it is necessary to refer 
it to the circle of reference in order to determine its phases; 
for example, in figure 1 if the point is found in P, its phase 
is described by the degrees of the arc MA or of the angle MOA. 
Epoch of a definite phase is the instant at which this 
phase takes place; v. g., if the arc MA is 60"" and it is two 
hours since the moving point was at the zero of time, we say 
that 2 o^clock is the epoch of the phase of 60"'. 

The initial phase is the angular distance between the zero 
of time and the starting place of the motion. In the figure, 
if we take the point N as the zero of time, the initial phase of 
the motion is given by the arc NA. 

Suppose that the point M requires a day to describe the 
circumference — that is to say, that the period of its motion 
is twenty-four hours. Starting from the zero of time N (fig. 2) , 
let us divide the circumference into 24 equal parts; draw the 
straight line NX equal to the length of the circumference; 
this line will also be divided into 24 equal parts, each one 
of which represents the length of an arc of 15% or one hour 
of time. Suppose the circle of reference to move along the 
straight line NX parallel to the vertical diameter with a motion 
equal to that of the moving point in the circumference; at the 
end of each hour the center will pass through each of the 

Fig. 2. 
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HARMONIC FORMULA OF FOURIER AND BESSEL. 7 

successive points in which the straight line was divided, and the point which moves 
on the diameter will describe the sine-curve of the figure. 

If it is desired to describe another harmonic curve of a period | of the former, the 
parts into which the circle of reference is to be divided must contain double the num- 
ber of degrees, for it is clear that, if the circumference is to be described in twelve 
hours, the point must pass through 30'' each hour. In like manner for curves of 

periods ^, i, of the first, the respective circles of reference must be divided 

into parts containing triple, quadruple, etc., the number of degrees; so that the periods 
of these curves are inversely to the number of degrees contained by the parts into 
which their circles of reference are divided. 

Any ordinate of the curve represents the elongation (distance of the moving point 
P from 0) in the corresponding moment. Call these ordinates y/, y/% etc. Then 
in the rectangular triangle POM (fig. 1) we have as a general expression of the 
ordinates — 

OP—y^OM sin a=u sin {U+x) 

where u is the radius of the circle of reference, or the amplitude of the curve, U the 
initial phase, and x the parts of the circumference through which the point M has 
moved at a given moment. 

With these considerations in view, we get for the different ordinates of the char- 
acteristic harmonic curves of any periodic curve the following expressions: 

l^iisl ordinate. Second ordinate. Third ordinaie. 

First curve .... n.smU, k, sin(í7i-(- ^) • ii, sin(t/,+2x) 

Second curve .. i/., sin U, u^út\{U.+2x) íl¿ ^m{Ur^-\Ax) 

Third curve ... t^.^ sin U,^ i/, sin(C7.5-|-3a') u,, sin(í7..,-f6aO 



Making an algebraic sum of homonymous ordinates we have the following expres- 
sions for the corresponding ordinates of the resultant periodic curve : 

Resultant* 

First ordinate = ?/, sin C/^+i^. sin U.+u.^Hin Í7.¡+ 

Second ordinate = u, sin(U^-{-x) -{-u^ sin(í7.+2ír) +?i.. sin(í7.^+3a') + 

Third ordinate == ¿/^ sin(f7,-f2x) -\-u. sin(í7. + 4ír) -\-u.^ mn{U:^-{-6x) -f 

which, as is readily seen, are nothing more than particular expressions of the general 
formula — 

¿f, sin(i/, + .T) -^iL, sin (U .~{-2x) ~\-u., sin(í7.,4-3^0 + 

in which the variable element x successively takes the corresponding values. 

Deducing BesseFs formula by this process, what was said in the beginning becomes 
clear, viz, that it is the analytical expression of a periodic curve in terms of its harmonic 
components ; each one of the terms representing one of these components. 

The constants u and U are usually called harmonic coefficients {u lineal, and U 
angular) , for they are the elements which characterize the harmonic motion in question ; 
u gives the amplitude, and consequently the importance of the motion, for as may be 
seen in the figures, according as u decreases, the motion of the point P on the diameter 
becomes slower; U (fig. 2) fixes the epoch of the different phases; an increase in its 
value accelerates the epoch of all the phases, and, on the contrary, a decrease retards it. 

It is of great importance to remember this property of the initial phase; for 
example, if in a station whose longitude was IS"" W. of Greenwich it was required to 
calculate in Greenwich mean time instead of in local mean time, the zero of time is 
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advanced one hour, and it is hence necessary to diminish the value of the initial phases 
15^ in L/„ 30"^ in L/^, 45'' in C/3, etc. Hence, it is important that in all cases the time 
employed should be explicitly stated. 

The quantity M which appears in the formula represents the mean of the observa- 
tions which enter into the discussion. What is developed in the trigonometrical series 
of the formula is the difference between each observation and the mean of all of them ; 
consequently, in order that the observations (or the ordinates of the observation curves) 
be represented according to this formula, it is necessary to add the quantity M to the 
development of these differences. 

What has been said is no doubt sufficient to give a general idea of the mechanism 
of the formula, and in the application which will be made of it in the study of the 
diurnal variation of the atmospheric pressure in Manila we will endeavor to make this 
idea more concrete. 

APPLICATION OF THE FORMULA TO THE ANALYTICAL STUDY OF THE DIURNAL 

VARIATION OF THE ATMOSPHERIC PRESSURE IN MANILA 

DURING THE PERIOD 1890-1909. 

In the preparation of the data which were to enter directly in the calculation the 
following method was followed : 

The hourly means of each month were taken from the official publications of the 
Observatory, and the mean of the zero hour, or first midnight, of the representative 
day of each month was calculated separately; this mean is the mean value of the read- 
ings of the first midnight of each day of the month. As the zero hour, or first mid- 
night, of a day is the last midnight of the preceding, in any month (January for 
example) the readings of the 0^ will be: For the 1st, that of the 12 p. m. of the last 
day of December of the preceding year ; for the 2d, that of 12 p. m. of January 1, etc., 
so that the values that enter into the mean of this hour are the same as those of 12 
p. m., except the reading for 12 p. m. of the last day of each month, which is substituted 
for the corresponding one of the previous month. 

It is evident, therefore, that, having calculated the mean of the 12 p. m. readings, 
it is sufficient, in order to calculate that of the 0^' readings, to subtract the 12 p. m. 
reading of the last day of the month in question from the corresponding reading of 
the preceding month, divide this difference by the number of days of the month, and 
add the result algebraically to the mean of the 12 p. m. readings. 

Having obtained these 25 hourly means, the general hourly means of the twelve 
months for the whole period have been calculated by grouping together Hke months, 
and taking the mean of the means obtained for each hour of a given month in different 
years. These general hourly means have been approximated to thousandths and are 
given in Table I. 

Table II is the same as I with the aperiodic element removed. The following will 
give an idea of this correction : 

Noncyclic or aperioddc variaUon, — In the explanation of the formula we saw that 
the elementary harmonic curves consist of 25 ordinates corresponding to the diiferent 
hours of the day, from 0^ to 24% both inclusive. The extreme ordinates are equal for 
the same curve, consequently the ordinates corresponding to 0'' and 24^^ in the resultant 
curve must be equal, in that they are the algebraic sum of the corresponding ordinates 
in different elementary curves. 

Hence, it will be understood why it is necessary to calculate the mean of the zero 
hours; for one must be quite sure that the curve to be analyzed as a periodic one has 
the essential characteristic of periodicity, viz, that its extreme ordinates be equal. In 
case this is not so, then some element has been introduced into the march of the phenom- 
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enon which perturbs its periodicity. This perturbation is what is called noncyclic or 
aperiodic variation. 

What is ordinarily the cause and origin of this variation will be understood by a 
familiar illustration. 

The temperature varies in the different seasons of the year, rising in spring and 
summer and falling in autumn and winter; so that in May, for example, although it 
may not be said in general that every day is warmer than the preceding, nevertheless 
the general tendency is toward a rise of temperature; and if we wish to calculate the 
mean diurnal variation of temperature for this month in an average year by assigning 
to each hour of the day the mean for that hour derived from the whole thirty-one days, 
we shall find that the second midnight of the representative day is somewhat warmer 
than the first midnight. 

In the case of barometric observations we come across an analogous phenomenon. 
It is well known that the atmospheric pressure is subject to an annual oscillation over 
and above the diurnal oscillation. In studying, therefore, the characteristics proper to 
the monthly diurnal variation of atmospheric pressure, we must begin by eliminating 
the aperiodic element. 

This elimination is usually performed by applying to each hour the correction — 

D{\2-h) 
24 

in which D is obtained by subtracting the value of the 0^ from that of the 24'\ and h 
represents the hour in which the correction is to be applied. This correction is logical 
only on the hypothesis that D comes in at a uniform rate throughout the twenty-four 
hours. 

To prove the legitimacy of this hypothesis presents difficulties ; nevertheless, when 
the investigation is limited to the twenty-four hours it is usually accepted, for it seems to 
be the only workable hypothesis. 

Let us see how the formula implicitly supposes this hypothesis. If the variation is 
supposed to be uniform during the twenty-four hours, it must be rerpresented by a 
straight line. The correction formula will, therefore, be the equation of the straight 
line — 

Take as origin of rectangular coordinates the zero hour and transfer as usual the 
hours to the axis X. In the formula (1) x represents the successive hours; y the cor- 
rection to be applied to each. When the correction has been applied to the extreme 
hours their values must coincide; therefore, the semidifference of their values is to be 
added to one and subtracted from the other. The correction for these two hours, called 
D their difference, will be — 

forthe24>^.„ ±/?-6+24a 



Eliminating 6, we have^ 



for the Qh. 



2 

D 

2 
D 
2A 



and substituting the values obtained for a and h in the equation (1) we get — 

^D , D _D (12-x) 
^ ' 2 24^^ 24 

Hence, we see that the hypothesis which requires the aperiodic variation to be uniform 
during the twenty-four hours is necessary and sufficient for the deduction of the correc- 
tion formula. 
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An example will demonstrate the influence which this correction has on the curves. 

Let us take the month of November, which in our case has the largest correction. 
The smooth curve (fig. 3) represents the variation without correction; the dotted one 
represents the same variation, but with the correction applied. 

From the form and position of these curves it is clearly seen that the elimination 
of the aperiodic element is reduced to revolving the curve through a small angle with 
12 a, m. as the center, the characteristic undulations remaining the same. 

Once this correction applied, the diurnal means were calculated for Table II by 
dividing the 24 hourly means of each month, from 0^ to 23,^ by 24, and Table III was 
formed by subtracting the diurnal mean of each month from each hourly mean. 

In this last calculation we have approximated only up to two decimal places; but, 
as in the previous calculations, the approximation was carried to three decimal places, 
the inaccuracy of the data in the last table cannot be greater than 0.005 and the sum of 
the positive values does not differ absolutely from that of the negative by more than two 
or three hundredths at most, which would not have been obtained if in the earlier cal- 
culations the approximation had only been carried to two decimal places. 

In the analytical study which follows we have used local mean time. 

All that is necessary for the determination of each of the elementary waves, into 
which the diumal barometric variation is to be resolved, are the lineal harmonic coeffi- 
cients M, and the angular harmonic coefficients Í7 — that is to say, the constants of For- 
mula II, for these are the ones which characterize those waves. 




Fig 3 

These constants cannot be found directly; they must be determined by means of 
the constants p and (/ of Formula I which are calculated directly. 
The equations of these constants are, as we saw before: 

p li sin U) /^x 
q ti cos U]^^ 
Hence, we obtain the relations : 

sin i7 , TT P • TT P TT Q 

r^tan C/ ^^; sm Í7^^^ ; cos [/ 
cos U q tí u 

Fvoxti the first of these equations we can deduce the angular constant U in terms of 
p and q; and from the last two in terms of i> or q and u; this last is obtained by the 
equation : 

which is easily obtained from the equation (1) by squaring and adding. All the diffi- 
culty consists in determining the constants p and (/. 

The method which Bessel elaborated for this, and which we will explain at the end 
of this study, is long and tedious, and, moreover, requires a constant concentration in 
the work if errors are to be avoided; this is owing to the great number of angles to 
be reduced, and the frequency with which the signs change in passing from one qua- 
drant to another. Father Dechevrens, S. J., says of this method: *'La méthode na- 
turelle n'est que pour les hábiles." This Father has published two pamphlets, one en- 
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titled **Méthode simplifiée dite des facteurs pour le calcul des series de Fourier et de 
Bessel appliquées a la Météorologie," in which by means of tables he simplifies consider- 
ably the calculation. The other pamphlet is entitled **Note complémentaire a la méthode 
simplifiée du calcul, etc." ^ in which he explains the manner in which the tables are 
formed. 

The first three of these tables serve for the determination of the coefficients p and 
q in series of 4, 6, 8, 10, 12, 24, and 36 observations. 

We think that it will not be without interest for those who may not be accustomed 
to this class of analysis, to work out in detail an example by Bessel's method, and see 
in what the utility of Father Dechevrens' tables consists ; and, even at the risk of being 
too detailed, we will follow step by step the method of Bessel in the determination of the 
constants p and q for January. 

The equations which give for the present case the values of these constants for the 
first eight terms of Formula (I) are:'-^ 

12pi=2am cos mx \2q\=%am sin mx 

12p2=^Sara cos m (2a;) 12qf2=2am sin m (2a;) 

12pr~=^(im cos m(3a:) 12g3=2ani sin m(3a:) 

12p4==Saa, cos m (4x) 12g4=Sam sin m {4x) 

The sign 5 indicates that the quantity which it affects is to be developed into a 
series of terms of the same form, which follow according to the series of natural num- 
bers 0, 1, 2, etc., so that we write; 

2a,n cos mx =ao cos + ai cos x+at cos 2x +a3 cos 3x + , 

X represents the value of one of the parts into which the circle of reference of the first 
elemental curve is divided. As the period of this curve is twenty-four hours, its circle 
of reference is divided into 24 parts, and consequently a:=15°. The periods of the other 
curves are i, i, { of that of the first, and therefore the divisions of their respective 
circles of reference equal 

2ir; 3x; 4a;; or 30 ; 45 ; 60 ; 

a,» represents the values of Table III corresponding to the different hours; thus for 
January — 

ao^=46, ai===17, a2^— 20, etc. 

According to this, in order to determine p, in January we have — 

12pi--46 cos f 17 cos 15-20 cos 30-45 cos 45' 

All that is then required is to find the values of these cosines in the trigonometrical 
tables and work out the operation. The work, which would be enormous if the multi- 
plication had to be made in each term, can be simplified in the following manner. 
The value of the cosines in the different hours are : 



0" ... 


cos 0° = 


1 


12" 


zz - 


"1 


1 . ... 


cos 15°=: 


.966 


13 


—. 


- .966 


2 


cos 30°=: 

COS 45° = 

cos 60° = 


.866 
.707 
.5 


14 





- .866 


3 ... 


15 





- .707 


4 ... 


16 


~. 


- .5 


5 -- 


cos 75° = 

cos 90° = 


.259 



17 ... . 


. 


- .259 


6 .... 


18 


...,. — 





7 ... 


.... cos 105° = -cos 75° = 


- .259 


19 


rr 


.259 


8 .... 


cos 120° = -cos 60° = 


- .5 


20 


~ 


.5 


9 


, "ZZ 


- 707 


21 


rr 


.707 


10 .... 


m 


- .866 


22 


.. ■. .. m 


.866 


11 


, = 


- .966 


23 


= 


.966 



Tipografía della Pace di Filippo Cuggiani, Roma, Via della Pace, num. 35. 
' The mathematical proof of these equations will be given later on. 
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Each of those numbers appears as a common factor of several hours, 
therefore : 



We have 



12pi=(0-~12)l+(l-ll--13+23).966+(2-10— 14+22).866+(3-9-15+21).707 

+ (4-8-16+20). 5+ (5-7— 17+19). 259 

Substituting for the hours which appear in parentheses the corresponding values taken 
from the Table III we get : 

12?>i=(46-15)l+(17-79+63+59).966+. 

In the determination of g, we obtain I 

0*' sin o°=o 

1 ^ sin 15^= .259 

2 sin 30°= .5 

3 - - sin 45"= .707 

4 sin 60"= .866 

5 - sin 75"= .966 

6 sin 90° = 1 

7 sin 105^=sin 75°= .966 

8 = .866 

9 - = .707 

10 = .5 

11 , = .259 



►y the same method: 
12*' 


.... - 


13 - 


.... =- .259 


14 . . .- 


=- .5 


15 


=- .707 


16 


— — .866 


17 


=- .966 


18 


-~1 


19 


=-~ .966 


20 


=:- .866 


21 


=- .707 


22 


=- .5 


23 


=- .259 



12gi 



Taking into account the common factor of several hours we write: 
-(1+11-13-23). 259 ¡(2+10-14-22). 5+ (3 +9—15-21). 707 



f (4+8-16-20).866+ (5i- 7— 17-19).966+ (6—18)1 

and substituting their values for the hours we find q,. 

For the determination of p.^ and q^ we proceed in the same way, remembering that 
the divisions of the circle of reference equal SO"", we say therefore: 



P2 



0" cos 0° 

1 cos 30° 

2 cos 60° 

3 cos 90° 

4 cos 120° = -cos 60° 

5 ....cos 150° = -cos 30° 

6 - 

7 

8 - 

9 ^.. 

10 

11 ^-..- ^ 
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12»' 


.866 


13 


.5 


14 





15 


- .5 


16 


-~ .866 


17 


-1 


18 


-~ .866 


19 


- .5 


20 





21 


.5 


22 


.866 


23 



= 1 

= .866 
= .5 
= 
= - .5 
= - .866 
= -1 
= - .866 
= - .5 
= 
= .5 
= .866 



12pg^(0-6+12-18)l+(l-5 7 i 11+13-17 -19+23). 866+ (2-4-8 

+104 14— 16-20+22). 5 

Q2 



0\...... 
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sin 
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.5 


13 


2 ....... 


sin 


60° = 


.866 
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sin 


90° = 


1 


15 
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sin 120° = sin 


60° = 


.866 


16 


5 
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.5 


- .5 

-" .866 
--1 

- .866 

- .5 


17 
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18 


7 .. ... 




__ 


19 
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20 
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n 


10 
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11 . .... 




— . 


23 



zr 


.5 


= 


.866 


zn. 


1 


zr: 


.866 


= 


.5 


= 





= - 


~ .5 


=: - 


~ .866 


m - 


-1 


= - 


~ .866 


= - 


- .5 
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12g2=(l+5-7~ll+13+17-19— 23).5-f(2— 4-8-10"rl4+16-20-22).866 

+ (3-9+15-21)1 
The circle of reference of the third curve is divided at every 45^ hence : 

Pa 
0" cos 0°= 1 12" 



1 „ cos 45"= .707 

2 cos 90°= 

3 cos 135° = -cos 45° = - .707 

4 =-1 



13 

14 
15 
16 



= -1 
= - .707 
= 
= .707 
= 1 



5 
6 
7 
8 
9 
10 
11 



= ~ .707 i 17 = .707 



I 18 
.707 i 19 

1 20 



= .707 
= 
= - .707 



21 
22 
23 



= 
= - .707 
= -1 
= - .707 
= 
= .707 



12p3= (0-4 +8-"12+ 16-20) 1 + (1-3—5+ 7+9—11—13+15+17—19-21 + 23) .707 

«3 



0'^ : sin 0° = 

1 sin 45° = 

2 sin 90° = 

3 sin 135°=sin 45° = 

4 = 



.707 



5 
6 
7 
8 
9 
10 
11 



.707 
= 
= - .707 
= -1 
= - .707 
= 
= .707 
= 1 
= .707 



12\ 

13 . 

14 - 

15 - 

16 . 

17 . 

18 . 

19 . 

20 . 

21 . 

22 . 

23 . 



= 
= - .707 
= -1 
= - .707 
= 
= .707 
= 1 
= .707 
= 
= - .707 
= -1 
= - .707 



fore: 



12g3=(l+3-5-7 t 9+11-13— 15+17+19-21— 23).707+(2— 6+10— 14+18— 22)1 
The divisions of the circle of reference of the fourth curve equal 60'' ; we have there- 



Pé 



0" cos 0° 

1 cos 60° 

2 cos 120° = -cos 60° 

3 

4 

5 

6 

7 

8 

9 

10 

11 



1 


12' 


.5 


13 


- .5 


14 


-1 


15 


- .5 


16 


.5 


17 


1 


18 


.5 


19 


- .5 


20 


-1 


21 


- .5 


22 


.5 


23 



= .5 

= - .5 
= -1 
= - .5 
= .5 

= 1 
= .5 

= - .5 
= -1 
= - .5 
= .5 



12p4=(0— 3+6-9+12— 15 M8— 21)1 + (1—2-4+5+7— 8-10+11 f 13-14— 16 + 

17+19— 20— 22+23) .5 

= 

= .866 

= .866 

= 

=- .866 

="- .866 

= 

= .866 

= .866 

= 

=- .866 

=- .866 



0' 

1 

2 

3 


sin 0°= 

sin 60°= .866 

sin 120°=sin 60°= .866 

= 


12' 
13 
14 
15 


4 


=- .866 


16 


5 


=- .866 


17 


6 .. ... 


= 


18 


7 


= .866 


19 


8 


= .866 


^0 


9 


— 


?1 


10 


=~ .866 


?? 


11 


=- .866 


23 
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12g4-:(l+2-4— 5+7+8-10— 11+13+14— 16— 17+19+20— 22— 23). 866 
We have descended to so many particulars because we wished to show that Father 
Dechevrens' tables for the determination of these coefficients could be thoroughly relied 
on, in that they are founded on BesseFs method, and, moreover, the operations most ex- 
posed to errors are already verified. The formulae given by Father Dechevrens are: 

pr-(0-12).083+ (1+23-11-13). 08+ (2+22-10-14). 072 

+ (3+21-9"-15).059+ (4+20-8-16). 042 + (5+19-7— 171.022 
^r--(6— 18). 083+ (1-^11-13-23). 022+ (2+10-14-22). 042 

+ (3+9-15-21).059+ (4+8-16-20). 072+ (5 + 7-17-19). 08 

i>2--(0 + 12-6-18).083+(l + ll+13+23~5-7~17-19).072 

+ (2+ 10+14+22-4-8-16-20). 042 
g2= (3+15-9-21). 083+ (1+5+13+17-7-11-19-23). 042 

+ (2+4+14+16-8-10-20-22). 072 

f p8 ^-(0+8+16-4-12-20). 083+ (1+7+9 fl5+17+23-3-5-ll-13-19-21).059 
^^^ Us (2+10+18-6— 14 -22).083+ (1 + 3+9 4 114 17+19-5-7-13-15—^ 



II 



rp4=.(o+6+12+18-3-9-15-21). 083 + (1 + 5+7+11+13 



IV- +17 + 19+23-2-4-8-10-14-16-20-22). 042 

U4=-(l+2+7+8+13+14+19+20-4-5-10~ll-16-17-22-23).072 

Comparing them with those found by BesseFs method we see that they are iden- 
tical; the coefficients of the parentheses are the same as Bessel's divided by 12. 

By means of these tables we have determined the constants which form Table IV. 

With these constants we proceed to the determination of the harmonic coefficients by 
means of the formulae which, as we have seen, join them. We give below the ones we 
made use of, and the method followed : 

^=tg Í7. '. log p— log g log tg U 
^ "sin Í7**' log^ logP~~^li>gsin U 
In this example for January we should, according to this, have: 

logp, 1.127429 
log gi- 1.738463 
log pi-logQi 1.388966-^log tg ÍA .'. Uv AS" 45' 37" 

log pi=-L 127429 
log sin J7i^.376320 
log pi-log sin CA--1^751109 -log in .\ t^i -0.5638 

One must remember in which quadrant U lies. This rule will help to verify it easily. 

If p and q are positive, U belongs to the first quadrant ; if p is positive and q negative, 
to the second quadrant; if both are negative, to the third; and if p is negative and q 
positive, then U is in the fourth quadrant. 

Thus to determine U. in this month we find directly in the tables the angle IS"" 
23'; but as p^ is positive and q^ negative, it must be reduced to the second quadrant, 
which gives 

1/2=164^ 37' 

In this way the values of m and U given in Table V have been deduced for each 
month. 

Although this calculation actually suffices for the discussion of the elemental curves, 
nevertheless we have calculated separately these four curves for each month. 
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The values of the ordinates of each hour together with those of the resultant, 
appear in Table VI, to which we have added another column containing the differences 
between the calculated resultant variation and the observed value. The smallness of 
these differences seems to be a good proof of the correctness of the calculation. 

Once those waves calculated, we thought that it would not be altogether without 
interest to add the 12 plates in each of which appear the elemental harmonic curves of 
a month, together with the resultant. In this manner it will be easier to appreciate 
the importance of each of the curves. From them it will be seen quite clearly why it 
was enough to calculate the first four terms of the formula; the series is rapidly con- 
vergent, so that the fourth sinusoid is almost confounded with a straight line. 

Although Father Dechevrens also gives tables for the verification of this last 
calculation, we thought it would be much shorter, if not easier, to work it out by the 
natural method. 

As an example, we will calculate the four components for January. The general 
expression for the ordinates of the component of twenty-four hours is 

III sin{Ui~\-x) 

consequently we have the general expression 

0.564 sin(13 46'+^) 

In order to determine the angle in the parentheses for the different hours, it is sufficient 
to note that in the first wave, the point which passes round the circle of reference, has 
described an arc of W at the end of first hour; 30'' at the end of the second; 45 "^ at the 
end of the third, etc., so that for the different hours we have : 



Oh sin 13 46^= .238 

1 sin 28 46^= .481 

2 sin 43 46^= .692 

3 sin 58 46^= .855 

4 „ sin 73 46'= .960 

5 sin 88 46^=1 

6 __^_ sin 103 46'=sin 76 14'= .971 



7h sin 118 46'= sin 6! 14'= .877 

8 sin 133^46'= sin 46 14'= .722 

9 sin 148 46'rrr sin 31 14'= . 519 

10 sin 163 46'= sin 16 14'= .280 

11 „-.. sin 178^46'= sin 1 14'=r .022 

12 sinl93 46'=-sinl3 46'=--.238 



After the twelfth hour, the values begin to repeat in the same order, but with 
opposite sign. It then suffices to multiply each one of the natural fuinctions, which, 
appear in the last member, by the coeflftcient 0.564 in order to obtain the values which 
represent the diflierent ordinates of this wave. 

Second wave. — The general expression is: 

U2BÍn{U2+2x) 
which is reduced to 

0.978 sin (164 37'+2a^) 
The variable element 2x increases 30° every hour; hence 

Oh ..„_. .^ . .... sin 164 37'= sin 15 23'= .265 

1 - ^ „ . „ „ „ . sin 194" 37'r=: sin 14 37'= - . 252 

2 _.. . „ sin 224 37'=-sin44 37'=--.702 

3 - _ - sin 254 37'=--sin 74^ 37'= . 964 

4 _ _ . . sin 284 37'=-8Ín 75 23'=--.968 

5 . „ „ _ _ . sin 314 37'=-8in 45 23'=-. 712 

6 _-„.-. sin 344 37'=-sin 15 23'= -.265 

After the sixth hour the values begin to repeat but with opposite sign, so that in 
order to obtain the ordinates of this curve, all that is necessary is to multiply these 
values by 0.978, and go on repeating this series of products till the twenty-third hour, 
being careful to change signs at each repetition. 
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The other two waves are obtained by an analogous process, taking into account the 
value of the variable element in each one. 

Fourth wave, 
0'' ..sin 130*^ 1T= sin 49" 43'= .763 

1 ..sin 190° ir=~-sin 10" 17' = ~-.179 

2 ..sin250° 17'=:-sin70° 17' = -.941 

3 ..sin 310° 17'= -sin 49° 43'= -.763 





Third 


wave. 




0»» 




sin 11° 1' = 


.191 


1 




sin 56° 1' = 


.829 


2 , 


.sin 101° r= 


sin 78° 59' = 


.982 


3 


-sin 146° 1' = 


sin 33° 59' = 


.559 


4 


.sin 191° 1' = 


sin 11° 1'=- 


-.191 



This is the natural manner of proceeding in the calculation; in reality it is much 
simpler, if we note that after the sixth hour in the first wave, after the third hour in 
the second, and second hour in the third wave the angles become the complementary of 
the former angles. Hence it is sufficient to determine the angles up to the fifth hour, 
inclusive, in the first wave ; up to the second hour in the second, and the first hour in the 
third wave ; and without paying attention to the rest, take at once from the trigonome- 
trical tables the value of the sine of each of those angles and their complementarles, 
which are found at the end of the same line. For example, for the zero hour in the first 
wave of January we find the value of the sine of 13 "" 46% namely 0.238, and at the end 
of the same line we find the value 0.971 which is that of the sine of the complementary 
angle 76'' 14% which we write in front of the sixth hour. 

In order to find the ordinate for the zero hour of January by means of Father De- 
chevrens' tables, we have to apply the following formula : 

(0-12) ,083+ (1 13) .08 f (2 44) .072+ (3-15) .059+ (4-16) . 042+ (5 47) .022+ (7-19) X —.022 

+ (8 20) X--.042+ (9 21) X—.059+ (10-22) X— .072+ (11-23) X— .08 

substituting the values for the hours in parentheses, and performing the proper operations. 
In the present case after having made the subtractions in the parentheses, we have 



.31X.083+.8X.08+1.05X.072+1.1X.059+1.04X.042+.89X.022~-.84X.022- 

—. 79X. 059— .58X.72—.2X. 08- 0.1341 



.87X.042 



By the natural method we can obtain the same result, after the values of the trigo- 
nometrical linra have been found, by one single multiphcation : 

0.238X0.564-0.1342 

Of course, Father Dechevrens' method is less exposed to errors, as it reduces the 
work to simple arithmetic ; but we think that a little practice in trigonometrical calcula- 
tion would be sufficient to save time by the natural method. 

Epoch of the critical phases {maxima and minima). — The determination of the 
epochs of the critical phases offers no difficulty. As every ordinate is expressed by 

tiBiniUnx) 

the necessary and sufficient condition in order that this value be a maximum is that the 
sine, which appears in the general expression, be equal to unity, which leaves as condi- 
tion for the maximum 

f/ i a: 90 or 450^ 

according as U is less or greater than 90'\ This being supposed, let us remember that 

the variable element x has for the successive hours 0, 1, 2, 3, /i, in the first 

curve, the values 0% 15^, 2x15^ 3x15%. ...- ^XlS"^; and in the same way for 

the second, third, and fourth curves, this variable has for any hour the value /ix30^, 
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/tX^S"", /ixSO'^, respectively. We will now determine the epochs of the first maximum 
for the four January curves. 

For the first we have as condition: 



for the second: 

for the third: 
and fourth : 



15 15 



164" 37' 1-/1X30" -=450\\/i=^^^^^ ^^"-.^Qh 31. 



30 



ir l'+/^X45^-=90\\/t==P 46"^ 



130^ 17'+/ix6(r=:450^*./i=5*^ 20"^ 



When we have determined the epoch of the first maximum, we have determined 
the epochs of all the other critical phases, because they are twelve hours distant in the 
first curve, six in the second, four in the third, and three in the fourth. The epochs 
of the first maximum for the four curves in the different months are given below. 



Epoch of the first maxvmum» 



Month. 



January 

February - - . 

March 

April 

May 

June 

July - 

Augrust 

September 
October _... 
November . 
December . 



First 


! 
Second 


wave. 


wave. 


H. m. 


H. m. 


5 5 


9 31 


5 47 


9 41 


5 23 


9 45 


5 3 


9 44 


4 10 


9 44 


3 11 


9 49 


2 57 


9 52 


3 9 


9 61 


3 25 


9 41 


4 16 


9 30 


4 20 


9 25 


4 26 


9 26 



Third 
wave. 



H. m. 

1 46 

1 44 

1 10 

49 



16 I 
20 



6 27 
6 52 
1 23 



1 40 



Fourth 
wave. 



H. m. 



DISCUSSION. 

In what has preceeded we have resolved the diurnal oscillation of pressure in 
Manila into a series of simple harmonic oscillations, the first two of which, viz, of 
twenty-four and twelve hour periods, are of incomparably more importance than those 
of eight and six hour periods, so that the principal characteristics of the diurnal oscilla- 
tion of pressure are represented almost exclusively by the first two. 

But this analysis of the curves is in itself nothing more than an artifice in calcula- 
tion, for every periodic curve can be treated by it. Hence, the component harmonic 
curves may have a purely mathematical meaning without corresponding to anything in 
reality. 

Helmholtz in treating of this class of analysis says: **It is only a mathematical 
fiction; admirable, because it renders calculation easy, but not necessarily correspond- 
ing with anything in reality." Nevertheless, speaking in general, the practical utility 
of this analysis cannot be denied. 

A periodic phenomenon can only be the result of periodic causes of periods the 
submultiple of that of the phenomenon; otherwise the periodicity would be disturbed. 
By means of this analysis, the curve is resolved into the only possible elementary ones 
with periods which are submultiples of the original one, and hence it is most probable 
that the partial causes which influence the phenomenon will appear either totally, 
or in their principal characteristics in the manner in which the elementary curves vary. 

127182 2 
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If it is required to ascertain whether or not two phenomena are related, the study of 
the characteristic harmonic curves of both is the only process by which we can come to 
this knowledge; the fact that there may appear no trace of mutual relation in the 
observation curves does not mean very much, for the influence which the one may 
exercise on the other may easily be disguised by other causes which are acting at 
the same time. 

Even though we may deduce nothing more from the comparison of two phenomena 
by means of this analysis than the persuasion that there is no relationship between 
them, our labor would not have been thrown away. 

Coming to the particular case of the barometer Dr. Julius Hann ^ says : 

I was myself convinced that all the attempts to explain the diurnal barometric oscillation by 
means of the daily variations of the meteorological elements at any one place, as attempted by 
Kreil, Blanford, Renou, and others, could lead to no conclusion; and I have published a series of 
papers giving a precise description of the phenomenon as manifested over the whole earth, at sea 
level as well as at all elevations for which observations exist, and I have endeavored to give the 
results in such a form as would be suited for the basis of a physico-mathematical theory. With 
this object I have represented all the results of observations in periodical functions, and have 
calculated the amplitudes and phase epochs of the individual waves, whatever was their period, 
and which, when combined together, produce the complex result of the daily barometric curve, such 
as presents itself to direct observation. * * * This is the only method which has enabled us to 
make material progress in the comprehension and explanation of the whole problem. 

What we propose to do principally is to compare the results obtained for Manila 
with those of other regions, and hence we will put the following discussion in such 
a form as to make this comparison easy : 

First wave.— One of the most fundamental facts that have been established in the 
science of meteorology is that the lower strata of the atmosphere receive their heat 
from the surface of the earth by means of conduction, and then transmit it to the layers 
immediately above by convection. This movement of convection is chiefly ascending; 
the air expands laterally in the higher atmosphere and causes a fall of the barometer. 
The opposite takes place with the cooling of the surface of the earth ; the air in contact 
with it is cooled, contracts, thus causing descending currents and a rise of the barom- 
eter; therefore by virtue of this variation of temperature there must be produced 
a maximum and a minimum in the atmospheric pressure about the hours of minimum 
and maximum temperature. 

The exchange of heat between the surface of the earth and the lower layers of 
the atmosphere is of a periodic character, the principal period being one of twenty- 
four hours. We say principal period, because this variation is not simply harmonic; 
so that if it were to be analyzed separately it would not be represented by a simple 
harmonic motion, but, in addition to a fundamental component of a twenty-four hour 

period, there would result others less important of twelve, eight, six, hour 

periods. 

The first wave in the analysis of the barometric curve is one of a twenty-four hour 
period, and invariably presents the maximum and minimum about the hours of minimum 
and maximum temperature. This general character of the curve makes it likely that it 
represents the principal character of the variation of pressure, caused by the exchange 
of heat between the surface of the earth and the lower layers of the atmosphere. 

Both the amplitude and the epochs of the critical phases of this wave present special 
characteristics in Manila, which it is well to note. The amplitude is less and the epoch 
of the maximum occurs earlier than in the greater number of the stations for which 



» Further contributions to the foundation of a theory of the daily barometric oscillation. Q, 
Journal of the Roy. Met. Soc, London, 1899, p. 40. 
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this analysis has been made. Hann made note of these peculiarities for Mauritius, 
Singapore, and Batavia in his "Untersuchungen über den taglichen Gang des Barometers/' 
Mr. Eliot includes in the following table these three stations, together with several 
others, studied by him in India, in order to compare them with Trivadrum, the only 
station in India which presents these characteristics. To these we add Manila for 
comparison. 



Station. 



Trivandrum 

Madras 

Bombay _. 

Rangoon 

Chittagong 

Calcutta (AHpore) 

Kurrachee . 

Aden .- _.. 



North 
latitude. 


Ampli- 
tude Ul 
for mean 


Epoch of Phase f/i 
max- j for mean 




of year. 


imum. 


oi year. 


o / 




H. w,. 




8 31 


0.413 


4 36 


21 2 


13 4 


.595 


6 1 


359 38 


18 54 


.472 


7 56 


330 53 


16 46 


.685 


6 20 


355 3 


22 21 


.582 


6 50 


347 45 


22 32 


.667 


6 56 


346 5 


24 47 


.497 


7 5 


343 49 


12 45 


.779 


7 26 


338 24 



Station. 



Allahabad 

359 38 i; Magpur i 21 

Lahore . „_ 

Mauritius , 

Singapore 

Batavia 

Manila 



North 


Ampli- 
tude Ul 


Epoch of 


Phase Ul 


latitude. 


for mean 


max- 


for mean 




of year. 


imum. 


of year. 


o / 




H. ?7?. 


o / 


25 26 


0.774 


7 38 


335 33 


21 9 


.767 


7 4 


343 52 


31 34 


.667 


8 32 


322 2 


20 6S 


.310 


2 54 


46 35 


1 17 


.5:n 


4 18 


25 85 i 


6 US 


,626 


4 19 1 25 17 


14 35 


.544 


4 15 ! 25 58 

i 



Mr. Eliot terminates his discussion of this first wave with the following words: 

The first component of the diurnal oscillation at Trivandrum, hence, exhibits characteristic features 
differing largely from those of other stations in India, but resembling strongly those at Batavia and 
Singapore and (to a less extent) Mauritius. The chief features are its small amplitude and the early 
epoch of the maximum phase. These are probably due mainly, if not entirely, to its maritime position 
and low latitude. 

The similarity which Mr. Eliot makes mention of as existing between Trivandrum, 
Batavia, and Singapore appears much more marked in the case of Manila, both in the 
amplitude and initial phase as well as the epoch of the maximum, and consequently 
Manila is an example which confirms the conjecture concerning these properties in coast 
stations of low latitudes. 

Annual variation of u^. — The annual variation of u^ and the annual variation of the 
absorption of heat in the lower layers of the atmosphere have between them the relation 
of cause and effect, as was pointed out by Hann ("Untersuchungen, etc/^ and Angot 
(**Annales du Bureau central Meteorologique de France,'' 1887, p. B. 297), and confirmed 
by the work of Mr. Eliot in Indian Meteorological Memoirs, Vol. XII, p. 374. 

A little consideration will show that the same takes place in Manila. In the '*Climato- 
logia" published by Rev. José Coronas in the '*E1 Archipiélago Filipino" it is shown that the 
months of February, March, and April are the driest and least cloudy months of the year, 
while July, August, and September are the rainiest and cloudiest. These circumstances 
are extremely favorable during the first three months for the absorption of heat in the 
lower layers of the atmosphere, and very unfavorable during the last three months, in 
that the clouds intercept a large amount of the calorific energy, and of that which passes 
to the surface of the earth, some is converted into latent heat of evaporation. 

A glance at the values of u^ makes it clear that the three months during which it 
reaches its maximum are February, March, and April, in which the absorption of the 
calorific energy in the lower layers of the atmosphere must be a maximum ; and it reaches 
its minimum values in July, August, and September when the absorption is least. 

Mr. Eliot mentions in the publication already cited that during the transition period 
between the rainy and dry seasons a relative maximum of the value of u^ usually takes 
place. This relative maximum appears in Manila in October which may be considered as 
the border month between the two seasons. 

The second wave, — There are two characteristics which are noteworthy in this wave ; 
first, that it is much more important than the others in all the months and in all the 
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regions where it has been studied, and, secondly, the constancy of the epoch of the critical 
phases. 

In the following table are given the epochs of the first maximum found by Mr. EUot 
for several stations of India : 



North 
latitude. 



8 31 

10 50 

13 04 

18 54 

24 47 

34 10 



Station. 



Trivandrum 

Trichinopoly 

Madras - -. 

Bombay- .--- 

Kurrachee ,„.,._..-,, 
Leh 

Mean value for tropical India . 

14 35 I Manila-.--—- 



Janu- 
ary. 

H m. 


Feb- 
ruary. 


March. 


April. 


May. 


June. 


July. 


Au- 
gruBt. 


Sep- 
tember. 


Octo- 
ber. 


Novem- 
ber. 


Decem- 
ber. 


H. m. 


H, 


m. 


H. m. 


H, m. 


H. 


w. 


H. m. 


H. m- 


H. m. 


H, w. 


H. m. 


H. m. 


9 30 


9 38 


9 


84 


9 36 


9 37 


9 


42 


9 47 


9 41 


9 29 


9 16 


9 14 


9 22 


9 50 


9 59 


9 


59 


9 58 


9 52 


9 


52 


9 53 


9 53 


9 48 


9 38 


9 39 


9 47 


9 45 


9 52 


9 


50 


9 49 


9 48 


9 


53 


9 57 


9 53 


9 39 


9 28 


9 27 


9 33 


9 43 


9 51 


9 


53 


9 54 


9 54 


10 


2 


10 6 


9 57 


9 45 


9 25 


9 23 


9 36 


9 41 


9 58 


9 


58 


9 67 


9 58 


10 


10 


10 9 


9 51 


9 49 


9 28 


9 29 


9 35 


9 48 


10 3 


10 


fi 


9 58 


9 52 


9 


50 


9 58 


10 22 


9 47 


9 49 


9 49 


9 49 


9 42 


9 53 


9 


50 


9 47 


9 47 


9 


54 


9 58 


9 55 


9 41 


9 29 


9 28 


9 36 


9 31 j 9 41 


9 


45 


9 44 


9 44 


9 


49 


9 52 


9 51 


9 41 


9 30 


9 25 


9 26 



This wave with its double daily oscillation is the one which produces the char- 
acteristic tides of the barometer. It is therefore of the greatest possible interest to 
discover its cause. 

The discussions of Hann, Angot, and Lord Kelvin on this point are too voluminous 
to reproduce here in all their details, so we will try and summarize them as shortly 
and clearly as possible. 

The phenomenon, represented by this wave must have the two properties, viz, that 
it have greater influence on the barometer than any other, and, secondly, that it be 
independent of local influences. 

From the outset, it is evident that the cause of this phenomenon can have no other 
origin than the influence of the sun on the terrestial atmosphere, "But,'' says Lord 
Kelvin, *4t cannot be owing to the gravitational tide-generating influence of the sun, 
since we would have a much greater influence of the same kind due to the moon, 
whereas in reality the lunar barometric tide is insensible or nearly so." It appears 
therefore certain that it must be attributed to the variation of the temperature. But 
this variation in the lower layers of the atmosphere is subject to all possible variations 
caused by local circumstances of the different stations, whereas the periodic action of 
solar rays on the higher layers, recurring day by day, must necessarily produce periodic 
movements of great regularity. 

This easily explains the regularity of this wave. What does not appear so clearly 
is the reason of its diurnal double oscillation and its great amplitude. The oscillation 
of the temperature in the higher layers of the atmosphere must be analogous in its 
general form to the oscillation of the temperature in the lower layers, since the causes 
which produce them are the same — the insolation during the day and the radiation 
during the night. The first causes the maximum to take place shortly after the sun 
has reached its maximum altitude, and the second produces the minimum shortly before 
sunrise. How, therefore, can this thermic oscillation of one maximum and one minimum 
give rise to the double barometric oscillation which we are studying? 

To solve this problem, Hann has recourse to the only means which exists for 
ascertaining the relation between two phenomena. He resolves the diurnal variation 
of the temperature into its harmonic components, and finds that of the two principal 
waves, the first (of twenty-four hour period) appears constantly with an amplitude 
greater than the second (of a twelve hour period) opposite to what takes place with 
the pressure. For example — 

Latitude 6" 0', equatorial Pacific... O.'^ST O.ni 

Latitude 16' 0', tropical Pacific l.^'dd 0.°81 

Latitude 33^ 0', subtropical Pacific 0.*'97 0.'*25 

Latitude 48^ 12', Temperate Zone (Vienna) 2.*'86 0.**51 
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In spite of this disagreement in the harmonic components of the two phenomena, 
Hann did not hesitate to assign as cause of the double oscillation of the barometer 
this double oscillation of the temperature. 

When this assertion of the eminent meteorologist was published, two difficulties 
were raised by scientists: First, that the thermic variation analyzed was that of the 
lower layers of the atmosphere, whereas the point at issue was the variation in the 
higher layers and of this it was questionable whether the oscillation was of exactly the 
same character as in the lower. The second difficulty was that this double daily tem- 
perature wave was only a result of calculation, a mathematical fiction, which could not 
serve in the explanation of a real phenomenon. There is in reality, no daily variation 
of temperature with two maxima and two minima; and if we, in spite of this, obtain 
a double daily temperature wave, because we insist on representing the daily march of 
temperature by a series of sines, this forced mathematical form can never serve to 
explain an observed phenomenon, whereas for this some real natural process must be 
sought for as a cause. 

The first difficulty is solved by what was said before concerning the variation of the 
temperature in the higher layers of the atmosphere; the causes which produce it are 
the same ; therefore it must be of the same general form as that observed in the lower. 
The solution of the second difficulty requires a little more explanation. It is known 
that if a limited mass of liquid or gaseous fluid be set in simple pendulumlike oscilla- 
tions, the amplitude of these oscillations are governed by the conditions proper to the 
fluid (its dimensions, temperature, etc.), so that if the mass of fluid receives a single 
impulse, it will always move with that class of undulations which is possible to it by 
reason of its peculiar conditions. These oscillations are commonly called free oscillations. 
If the impulse recurs periodically, then the mass is forced to take up oscillations of 
that period, even though these new oscillations do not coincide with any of the forms 
of oscillation which belong to the free waves. If this periodic impulse is not of a simple 
harmonic character — ^that is to say, if it cannot be represented by a simple sine curve — 
then it can be resolved, according to Fourier, into its harmonic components of periods 
1, I, ^, etc., so that the effect due to the impulse would be the same as that produced by 
several simple harmonic impulses acting together. It is quite immaterial to the question 
whether or not these component impulses exist in practice; it is sufficient if the real 
impulse can be considered as the resultant of them. 

If the period of one of these component impulses coincides with that of the free 
waves of the mass of the fluid, then the resulting forced oscillation will attain a dis- 
proportionately great amplitude; so that although the amplitude of this impulse be rel- 
atively small compared with that of the other .component impulses, the effect produced 
by it is much larger; whence it appears that if the resulting complex oscillation be 
analyzed by means of a sine series, the harmonic curve corresponding to the free oscil- 
lations will be of large amplitude, while the corresponding wave in the analysis of the 
impulse will be of small amplitude, and yet, that is no reason why the second cannot be 
the cause of the first. 

To come now to the case before us : The curve which represents the diurnal variation 
of the temperature in the higher strata of the atmosphere must necessarily be asymmet- 
rical, because its two causes, insolation and radiation, which depend on the duration of 
day and night, are asymmetrical. If, therefore, on analyzing it, we find that it is repro- 
duced almost exclusively by the combination of the two harmonic components of periods 
of twenty-four and twelve hours, then we may consider the real variation which we ac- 
tually observe as the resultant of two harmonic variations of the temperature ; one of a 
twenty-four hour period and the other of twelve; and it does not really matter whether 
or not there actually exist two causes, which give rise to these two variations; it is 
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enough for our proposition that the phenomenon be reproduced in the same manner as 
if they really existed. 

Doctor Margules, following the indications made by Lord Kelvin, undertook the 
difficult and laborious task of calculating the oscillations caused in the atmosphere by 
its periodic heating and cooling.^ In the first of his papers he showed that the period 
of free oscillation of our atmosphere is nearly twelve hours. 

These studies and conclusions of Doctor Margules have been received with a good deal 
JO uois-BODO aqq. ua^ OA^q puH uisBtsnq^uo joseveral interesting publications such as Paul 
Jaerisch's "Zur Theorie der hydrodynamischen Gleichungen in spharischen Koordinaten" 
(Meteorologische Zeitschrift, 1907, pp. 481-498) ; Trabert's ''Die Theorie der taglichen 
Luftdruckschwankung von Margules, und die tagliche Oszillation der Luftmassen" (Met. 
Zeit, 1903, pp. 481-501) ; W. J. Humphreys' "On the diurnal variations of atmospheric 
pressure'' (Bulletin of the Mount Weather Observatory, 1912, vol. 5, part 2, pp. 132-156) 
and the most complete of all, Horace Lamb's '* Atmospheric Oscillations" (Proc. Royal 
Society, 1911, vol. 84, pp. 511-572).^ 

We can, therefore, explain the double oscillation and the great amplitude of the second 
wave thus: When the atmosphere receives a thermic impulse (if we may make use of this 
expression) there is produced in it the same effect as would be produced by two simul- 
taneous impulses which could be reproduced by sine curves of twenty-four and twelve 
hour periods, respectively. The amplitude or intensity of the second is less than that of 
the first, but the oscillations which it causes in the atmosphere must be of great amplitude 
because its period coincides with that of the free oscillations ; while the oscillations pro- 
duced by the first are considerably reduced because their period does not coincide with 
that of the free oscillations. 

The hypothesis that the second wave is due to the absorption and radiation of heat in 
the higher strata of the atmosphere is confirmed by the comparison of the annual variation 



'"Über die Schwingungen periodisch erwarmter Luft." (Sitzungsberichte, Wien. Akad. Marz, 
1890.) 

"Luftbewegungen in iener rotirenden Spharoidschale." (Ibid. April, 1890; Jan. and Dec, 1893.) 

^ This work was already in the press when the Observatory was in receipt of Prof. Frank H. 
Bigelow's "The Thermodynamics of the Circulation and Radiation of the Earth's Atmosphere" published 
as Bulletin No. 3 of the Argentine Meteorological Office. On page 15, he proposes two difficulties 
against Lord Kelvin's theory; first, that the semidiurnal pressure waves disappear within 2,000 meters 
of the surface; and, secondly, that the heat term {Qi~-Qq) is not included in the general equation of 
motion as employed by Margules, Jaerisch, and Gold. It is not the first time that Professor Bigelow 
has impugned this theory; see, for example, his "Studies on the diurnal periods in the lower strata of 
the Atmosphere" (Monthly Weather Review, 1905, p. 93) and chapter 9, p. 458, of his "Report on the 
international cloud observations." At the same time it is not to be supposed that scientists like Lord 
Kelvin and Hann were ignorant that the semidiurnal oscillation diminishes with altitude till it dis- 
appears or becomes almost imperceptible, and yet, notwithstanding this, these scientists did not hesitate 
to establish their theory. Hence, it must be supposed the semidiurnal oscillation reappears in the higher 
atmosphere. Mr. Ernest Gold seems to indicate this in his interesting memoir "The International Kite 
and Balloon ascents" which won the first prize offered by the German Meteorological Society in 1908. 
For, on page 120, he says: "Hence changes of pressure at the earth's surface are connected with changes 
in the opposite direction in the stratosphere, which give rise to currents flowing between the stratosphere 
and troposphere." 

The second difficulty would seem to indicate that the three eminent mathematicians had over- 
looked the influence which the free heat tei-m could have on the general equation of motion, whereas 
Doctor Margules, after introducing into the differential equation of motion the free heat term, proceeds 
thus; "For our purpose it will be more convenient to consider the pressure variations as a consequence 
of the temperature variations, not as a conseciuence of the variable flow of heat. We, therefore, return 
to equation (4a)," where the term of which Professor Bigelow speaks certainly does not enter. 
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of these two factors and that of the amplitude m^. A brief explanation of this will suffice 
for our purpose.^ 

One of the most important factors in the influence which the absorption and radiation 
may have on changes of pressure, is the presence or absence of clouds in the atmosphere. 
The amount of heat absorbed by a dry and clear atmosphere is employed in expanding it, 
thus causing lateral and ascending currents, which are shown on the barometer by changes 
of pressure. On the other hand, the calorofic energy which reaches an atmosphere 
charged with clouds is converted, in great part, into latent heat of evaporation, which 
can not have any effect on the pressure. Hence, in the dry season the changes of pressure 
due to these causes will be greater than during the wet season, on account of the great 
difference of cloudiness between these two seasons. 

This difference in the values of u^, which are constantly greater during the months 
constituting the dry season than during the months of wet season, is seen quite clearly 
in Table V, as is also the difference in values in the perfectly defined separation of the 
two seasons from May to June. 

We will terminate the discussion of this wave by comparing its principal characters 
in Manila with those of the stations mentioned in connection with the first wave. 

The following table is copied from Volume X of Indian Meteorological Memoirs, page 
74, with the exception that we have converted the original observations in inches to milli- 
meters in order to facilitate the comparison : 



Station. 



Trivandrum _ 

Madras 

Bombay 

RanÉTOon 



Chittagong . 

Calcutta 

Kurrachee . 

Aden 

Allahabad . . 
Magpur 

Lahore 

Mauritius ._ 

Singapore -- 

Batavia 

Manila 



g 

C 



tnm. 

LlOO 

1.115 

.977 

1.023 

.946 

1.028 
.797 
.908 
.918 

1.010 

.628 
.733 
.995 
.959 
.919 



164 14 

168 5 

156 51 

153 3 



148 

147 
155 
159 
150 
154 

138 

163 

156 

159 

160 



5 

50 
40 
59 
13 
44 

18 

44 

18 * 

56 I 



H.m. 
9 32 
9 44 
9 46 

9 54 

10 4 

10 4 

9 49 
9 40 

9 59 j 
9 51 j 

10 23 j 

9 32 I 

9 47 

9 40 

9 40 



6 . 
.2-5 



1.4 

1.2 
1.5 

1.4 

1.3 

1.3 
1.6 
1.5 
1.2 
1.4 

1.3 

1.2 

1.3 

1.2 

1.4 



February 
March __. 
February 

March 

April 

March 
February - 
November 
March 
do 

February 

October- , 
(March 
\ April.. 

October 
/March 
lApril - 




November _ 

....do 

--do._.. 

j2 :/ - -do --.^. 

[IDecember . 

19 November 

.Ao 

October 



o 
1^ 



! July. 
Do. 
Do. 

j-Auisrust. 

{March. 
¡lAugrust. 
I July. 
j June. 



-do i April. ' 

November..! August. 
October . . . | July, 
.-..do. . j/February. 

11 July. 
November . ^ | February. 

do -.._„ July. 

do February. 

i do iJuly. 

iDecember ..!\ August. 



From this table one can appreciate in a general way the constancy of the principal 
characteristics of this wave in the diiferent stations; the almost perfect agreement in 



^ See a more detailed discussion of this point in Vol. XII of "Indian Meteorological Memoirs/' 
pages 328-347. 

The following may also be consulted: "Researches on solar heat and its absorption by the earth's 
atmosphere'^ by Langley in No. 15 of "Professional Papers of the Signal Service,'' and in Vol. I of 
"Annals of Astrophysical Observatory of the Smithsonian Institution/' and "Preparatory Studies for 
Deductive Methods in Storm and Weather Predictions" by Cleveland Abbe. 
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the ratio between the maximum and minimum amplitude, and the epoch of the year 
in which the minimum tak^ place, is especially evident, 

A few discrepancies will, however, be noted, which show that this wave is not quite 
free from the influence of local circumstances. Thus, for example, in some regions the 
amplitude has two maxima and two minima in the year, while in others, Manila among 
them, this double variation disappears. In stations like Batavia, Mauritius, and Aden 
which have the absolute maximum in October or November, there is another relative 
maximum in March or April. 

The reason of these variations is fairly obvious. The second wave, as it is given 
by analysis, must be a complex wave resulting from the interference of two with 
the same period of twelve hours. One of these (which we will call principal semi- 
diurnal wave) is the product of the variations of temperature in the higher strata of 
the atmosphere, as was explained above; the other (secondary semidiurnal wave) rep- 
resents the oscillations which the semidiurnal component of the temperature in the 
lower strata must produce in the pressure. The first of these two presents, according 
to Hann and Angot, a double variation during the year, with maxima at the equinoxes 
and minima at the solstices. The second, which is more irregular and of less amplitude 
than the first, owing to local influences, has one maximum only (during the dry season) 
and one minimum (in the wet season) . 

It is easy to see that the interference of these two waves can cause the complete or 
almost complete disappearance of the double annual variation in the wave we are dis- 
cussing; for the autumn equinox occurs usually in the wet season in the great majority 
of the stations, and the winter solstice generally occurs about the beginning of the dry 
season, so that the maximum and minimum corresponding to these epochs in the principal 
wave are held back either entirely or in great part, because they coincide with the 
epochs of the minimum and maximum amplitude of the secondary wave. The opposite 
happens in the summer solstice and the spring equinox, because they coincide with 
the wet and dry seasons respectively. 

M. Angot, in the publication already cited (p. B311 sq) proposes an ingenious 
method for the separation of these two waves. It would take up too much space to explain 
the whole proceeding, so we simply give in the following table the results obtained by his 
method for Manila.^ 

* The formula we finally obtain by this method of Angot is 

0.mm926 r^^ C08*X COS (2m+64^---2e) 

where h represents the mean pressure; S the declination of the sun; r the radius of the earth*» 
orbit; x the latitude of the place; m the different hours expressed in degrees; e the equation of 
time. The values which we have taken for these quantities are: 

cos 25 
2f r2"~ 

January +4° 48' , 0.9033 

February .,. +7" 01' 0.9722 

March +4" 19' 1.0043 

April ...-. .-....-- +0" 02' 0.9597 

May -r 46' 0.8734 

June -f 0^ 08' 0.8197 

July ......^ +2^ 41' 0.8424 

August - +1^ 50' ...-, 0.9208 

September -2'' 32' 0.9852 

Octol^r -r 00' 0.9806 

November -7^ 08' 0.9198 

December -1^ 56' 0.8746 
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Hour. 





1 

2 

3 

4 

5 

6 

7. 

8. 

9. 

10 

11. 

12 

13. 

14. 

15, 

16. 

17 

18. 

19 

20. 

21. 

22 

23 



January. 



Feb- 

ruary. 



+38 
+ 1 
—36 
—63 
—73 
—64 
—38 

— 1 
+36 
+ 63 
+73 i 
-I 64 I 
H 38 ■ 
+ 1 ' 
—36 
-«3 
—73 
-64 i 
—38 I 

- 1 Í 
4 36 
+63 

i 73 
n 64 \ 



+43 

-t- 4 
—36 
-66 
—79 
—70 
—43 

— 4 
+36 
+66 
+79 
+70 
-4 43 
+ 4 
—36 
—66 
—79 
—70 
—43 

— 4 
H 36 
-f66 
+79 
H70 



March. 



+41 



—40 

-70 

--82 
—71 
—41 


+40 
-f70 
+ 82 
+71 
+41 


-40 
-70 
^82 
—71 
—41 


+40 
+70 
+82 
+ 71 



April. 



+34 

— 5 

—44 
~-70 
—78 
—65 
-34 
j^.. 5 

+44 
+70 
+78 
+65 
+34 

- 5 
—44 
—70 
—78 
—65 
—34 
■+• 5 
4-44 
470 
+78 
165 



May. 



+29 

— 7 
-41 
-65 
-71 
—58 
—29 
+ 7 
+41 
+ 65 
+71 
+58 
+29 

— 7 
-41 
—65 
—71 
—58 
—29 
+ 7 
+41 
+65 
+71 
+58 



June. 



+29 

- 4 
-37 
—60 

-66 
-55 
—29 
+ 4 
+37 
+60 

+ 55 
+29 

— 4 
—37 
—60 

-66 
—55 
—29 
+ 4 
+37 
+60 

H 66 
4 55 



July. 



+33 

- 2 
-36 
-60 
-68 
—58 
—33 
+ 2 
+36 
+60 
+68 
+58 
+ 33 

— 2 
—36 
—60 
—68 
—58 
—33 
+ 2 
+36 
+60 
-^68 
+58 



August. 


Septem- 
ber. 


+35 


+32 


— 3 


- 9 


—40 


-48 


-66 


--73 


—75 


-80 


—63 


-64 


—35 


-32 


4- 3 


+ 9 


^ 40 


+48 


466 


+78 


+76 


-]-80 


+63 


4 64 


+35 


+32 


— 3 


~ 9 


-40 


—48 


—66 


—73 


—75 


-80 


-63 


-64 


-36 


-32 


4- 3 


+ 9 


4 40 


448 


466 


+ 73 


4 75 


-1-80 


4-63 


4-64 



October. 



+26 
—16 
-62 
—75 

-78 
-60 
-26 

+15 
+62 
+75 
478 
4 60 
+26 
—15 
-52 
—75 
—78 
—60 
—26 
+15 
+52 
4'- 76 
4 78 
460 



Novem- 


Decem- 


ber. 


ber. 


+24 


+29 


-15 


- 7 


-49 


-42 


-71 


-65 


-73 


—71 


-55 


-68 


—24 


-29 


+15 


+ 7 


+49 


-f42 


4'-71 


+65 


4-73 


4-71 


+ 56 


+58 


+24 


429 


-15 


- 7 


-49 


—42 


-71 


-65 


-73 


—71 


-55 


—58 


-24 


—29 


H-15 


4- 7 


4-49 


+42 


4-71 


-^65 


+73 


■h71 


+55 


1-58 



Note.— Black-face type indicates maximum or minimum. 



We cannot compare our results with those of other regions because we do not know 
whether this analysis has been worked out for others than the nine calculated by 
Angot, and for these Angot took only one or two months; nevertheless, we think it 
useful to have proved for Manila this character attributed to the principal semidiurnal 
wave. 

Third wave. — There are only two possible hypotheses with regard to the origin of 
this wave, viz: 

1. That it is produced by an independent action of an eight- hour period* 

2. That it is part of the effect due to the causes producing the first wave; since, as 
was indicated in the discussion of this wave, its causes cannot be represented by a simple 
harmonic motion. 

The first objection against the first hypothesis is that no independent physical action 
of an eight-hour period is known. And it may be added that the amplitude of the third 
component is very small, resembling in character that of a residual rather than the total 
effect due to an independent cause. 

In confirmation of these general reasons is the property which is noted in this wave 
in all regions, namely, that the annual variation of its amplitude and the epoch of its 
critical phases follow a law like to the one for the first wave, with this particular, that 
the sign of this variation is almost absolutely opposite, which raises the suspicion that 
the third wave is a correction of the values of the first. 

Given the slight influence which this as well as the fourth wave can have on the 
barometric variation, we believe that it is scarcely worth while discussing the matter. 
We do not wish, however, to omit a peculiarity of the third wave, noted by Doctor Hann, 
which is, that toward the equinoxes the epoch of the phases undergoes a notable inversion, 
and that this inversion coincides exactly with the epoch of the minimum amplitudes. 
This property appears well marked in Manila ; the inversion, spoken of by Doctor Hann, 
takes place from April (0^ 49«^) to May (7*^ 15"^) and from August (6^ 52»") to Septem- 
ber (1^ 23)°»; the minima of the amplitude also take place in April-May and August- 
September, as may be seen in Table V. 

CONCLUSION. 

From what has been said up to the present, it is deduced that the modern theory 
of the double daily oscillation of the barometer supposes two causes of the phenomenon. 
The one is quite free from local influences, and is precisely the one which introduces into 
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the barometric curve that double oscillation which has always been so mysterious to 
those who have examined its cause. It is represented by what we have called the prin- 
cipal semidiurnal wave, and is due to the variation of the temperature in the higher layers 
of the atmosphere. The other varies according to the local circumstances of each region, 
in that it is due to the variation of temperature in the lower layers. The curve represent- 
ing it is usually called the thermic wave. In the following table we give its values for 
all the months of the year; they are obtained by subtracting the principal semidiurnal 
wave from the observed variation (Table III) : 

Diurnal or thermic "wave" (unit 0.01), 



Hour. 





1 

2 

3 

4, 

5 

6 

7 

$. 

9 

10 

11. 

12. 

18. 

14. 

IB 

16 

17 

18. 

19 

«0 

21 

23 



January. 


Febru- 
ary. 


-f 8 


+ 6 


^ 16 


-i^ 15 


+ 16 


+ 22 


1' 18 


f 30 


+ 33 


+ 41 


-+ 44 


1 53 


+ 54 


+ 68 


+ 64 


~i-7ft 


+ 7« 


+ :« 


H" 72 


+ 75 


+ 49 


h 56 


! 15 


i 20 


-23 


~ 16 


-64 


-65 


-89 


-80 


^«2 


~~l>l 


-71 


— 77 


— 45 


™-57 


-"30 


- 44 


-20 


-39 


— 11 


-■31 


- 7 


— 23 


- 9 


— 18 


- 5 


— 8 



March. 



•I 18 
- 1- 22 
+ 28 
i 35 
+ 48 
t- 58 

4 84 
-184 

f- 74 
+ 48 

I 22 

- 11 
-50 
-79 

- 88 

- 74 

- 64 

- 57 

- 45 
-27 

- 12 
f 6 



April. 


-f 


27 


+ 


35 


+ 


34 


+ 


40 


^. 


53 




67 


+ 


72 


■■}'■ 


Ká 


,}, 


H2 


+ 


68 


f 


46 




18 


— 


9 


_- 


53 


— 


84 


— 


99 


™] 


104 


™ 


96 


— 


76 


™^ 


59 


_- 


48 


— 


25 


i- 


5 


4- 


22 



lay. 


June. 


July. 


Aug^ust. 

+ 25 


Septem- 
ber. 


October. 


Novem- 
ber. 


Decem- 
ber. 


+ 33 


f 39 


+ 32 


T 25 


+ 24 


+ 22 


+ 18 


+ 42 


-1-40 


+ 26 


+ 25 


+ 28 


+ 23 


+ 25 


+ 25 


} 35 


1 32 


+ 26 


+ 27 


r 31 


+ 27 


+ 23 


+ 17 


r38 


+ 38 


+ 30 


+ 29 


+ 31 


+ 31 


+ 16 


+ 16 


1 46 


I 38 


+32 


+33 


[33 


+ 29 


+ 22 


+ 26 


! 50 


-f 37 


+ 29 


+ 32 


+ 33 


+ 30 


+ 26 


+ 34 


+ 53 


-t-41 


+32 


- 29 


-Í-27 


+ 40 


+ 39 


+ 38 


■rol 


-h 38 


h30 


-]-■ 25 


+ 26 


-43 


+ 44 
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Note. — BIack*face type indicates maximum or minimum. 

In Plates XIII and XIV may be seen the curves which represent the thermic wave 
and the principal semidiurnal wave for several months. 

These are doubtless the two waves mentioned by M. Angot in his "Traite élémentaire 
de Météorologie/' page 104, in his explanation of the diurnal variation of the barometer. 

Without knowing exactly of what waves is being treated, it is difficult to grasp 
thoroughly this explanation, hence we have thought it good to translate it in this place; 
for, besides making it easier to understand with the curves before us, it will also serve as 
a discussion of our curves. 

*'The diurnal variation of the barometer has not as yet been completely explained. The 
most recent researches seem to indicate that it is the result of the superposition of two 
distinct oscillations^ — one which gives two maxima and two minima each day and whose 
period is consequently half a day, and which is called, therefore, the semidiurnal oscilla- 
tion or wave; the other which gives only one maximum and one minimum in the twenty- 
four hours, and which is called the diurnal wave. 

"The two daily maxima of the semidiurnal wave ought to be equal to one another, and 
the two minima likewise equal to each other. In all the seasons of the year and in all 
parts of the earth the two maxima take place about 10 o'clock, morning and night; the 
two minima about 4 o'clock, morning and afternoon. The difference between the maxima 
and the minima of the semidiurnal wave — that is to say, the amplitude of this wave ' 
varies but little for the same station during the course of a year; it is greatest at the 
equinoxes and smallest at the solstices. This amplitude does not depend on the topo- 
graphical conditions, but solely on the latitude— about 2 mm. at the equator, it then 

' M. Angot understands here by the amplitude of the wave, the amplitude of a complete oscillation; 
other authors and Angot himself, in the publication we have quoted several times, takes as the amplitude 
of the wave, the amplitude of a semioscillation, or the maximum ordinate of the curve. 
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diminishes slowly up to the Tropics, and afterwards very rapidly so that in high latitudes 
it is very small. It is probable that this semidiurnal wave is caused by the action of 
solar heat on the whole mass of the atmosphere, but up to the present time the manner 
in which it is produced has not been explained satisfactorily. 

**While the semidiurnal wave is the same for all regions situated in the same latitude, 
the diurnal wave has a very variable amplitude, which depends on the topographical con- 
ditions of the region, and which may be very different even in neighboring stations. 
In general, it is much weaker in coast stations than in continental stations, and its 
amplitude is in direct ratio to the amplitude of the diurnal variation of the temperature. 
The hours of the maximum and minimum vary greatly according to the countries and the 
seasons of the year ; the maximum takes place earlier in summer than in winter ; its mean 
epoch is about 9^ ; the minimum, on the contrary, whose mean epoch is about 16^, takes 
place earlier in winter than in summer. The period of falling, therefore, lasts on an 
average only seven hours, while the period of rising extends to seventeen hours ; the riKse 
is not uniform during this time for it is rapid immediately after the minimum, is slow 
during the night, and then increases rapidly again between sunrise and the maximum. 

"While the theory of the semidiurnal wave has not been explained as yet ' the diurnal 
wave, on the contrary, is fully explained by the local effect of the daily variation of the 
temperature. During the night, the lower layers of the atmosphere contract as they are 
cooled and, in order to fill the space thus produced by the contraction, air from all sides, 
especially from above, flows in. After sunrise the lower layers are warmed rapidly and 
tend to expand ; but the heating is produced first in the layers in contact with the surface 
of the earth, so that the layers of air between the superficies and the higher layers which 
have not been heated as yet, do not at first increase their volume appreciably; just as 
happens when the air in a closed chamber is warmed, the rise of temperature is manifested 
principally by an increase of pressure in the same way the first effect of the diurnal 
heating is an increase of pressure in the layers which are in contact with the surface of 
the ground, and the barometer begins to rise. 

''A little later, between 7 and 11 o'clock in the morning, according to the countries and 
season of the year, the heating of the lower layers is sufficient to determine a general 
ascending movement of the air above the station; the air expands, but as soon as the 
column of air situated above the heated region reaches a greater height than that which 
covers the neighboring regions, the air spreads out at the sides ; the quantity of air above 
the station diminishes, therefore, and the barometer falls. Toward the end of the day, 
when the temperature begins to fall rapidly in the lower layers, the ascending current loses 
its activity, ceases, and is finally replaced during the night by an inverse movement ; the 
barometer then rises, rapidly at first while the temperature is falling rapidly, and then 
slowly till sunrise." 

It will be noted that the explanation of the thermic or diurnal wave given by M. 
Angot coincides in great part with the first theory of the diurnal variation of the barom- 
eter proposed by Raymond, and perfected by Kreil and Renou, and which was unsatis- 
factory when it was applied to the total diurnal variation of the pressure. 

We may conclude, therefore, that although the analysis of the curves which represent 
the different meteorological elements of each region is a long and tedious one, it can, 
nevertheless, be of the greatest utility, both for the simultaneous study of various 
phenomena, the mutual relations of which we wish to investigate, as well as for the 
deeper study of the causes which may intervene in each of them. 

* We do not know whether M. Angot had any notice of the investigations of Doctor Margules, 
because neither here nor in the work which he published in "Annales du Bureau Central de France," 
does he make mention of them, notwithstanding their prime importance in the explanation of this wave. 
Doctor Hann, in publications posterior to those of Angot, laments the fact that the work of Doctor 
Margules was so little known, and that it had not even called the attention of the ''Meteorologische 
Zeitschrift." 



DEDUCTION OF BESSEUS EQUATIONS FOR THE DETERMINA- 
TION OF THE COEFFICIENTS OF FORMULA (I). 



Suppose that there are any number n of observations. 
Let tto, a^, 0^2, etc., represent the observations at 0, 1, 2, etc., hours. 
We saw on page 7 that these observations are represented by the following equa- 
tions:^ 

aiF=p+p\+P2+V%+ 

ai=p+pi cos x+qi sin x+p2 cos 2x+q2 sin 2x+p5 cos Sx-fga sin 3x+ 

(X2=p-f^i cos 2x-\-qi sin 2x-\-pz cos 4x+^2 sin 4x+p3 cos %x-\-q% sin 6x+ 



an-\=p-\-p\(^o% {n-~l)x^qiún (n~l)x+p2C0s2(n— I)x+g2sin2(n~-l)a: 

+pzCOS^{n~-l)x+qzsmZ{n—l)x+^ 



If the number of coefficients which we propose to determine be the same as the 
number of observations, we would have a system of as many equations as unknowns, 
and we could resolve it easily by the ordinary algebraic methocfe* But ordinarily the 
number of observations greatly exceeds the number of the coefficients we require; thus, 
for example, in the analysis of the barometric curve, we have twenty-four observations 
and we employ only eight of these coefficients. 

In this case the system is more than determined; and in the impossibility of finding 
the valu^ p and q, which substituted in the equations, reproduce exactly the observa- 
tions, we must content ourselves with estimating which are the most probable values 
of p and q, which give for the second members values compatible with the observations. 

In the theory of least squares it is proved that the most probable values are those 
which make the sum of the squares of the residual errors a minimum. The residual 
error is the departure from the observed value of the value obtained in substituting in 
the equations the most probable values of the unknowns. 

According to this definition, representing the residual errors by the symbols v^, v^, 
v^, etc., and the most probable values of the coefficients by Po, Pu P2f ^tc, we may write : 

i;o=-~aoH-p4-Pi+/>2+jt>3+ 

i;i=-~a]+p+pi COS x+q\ sin x-^-p^ cos 2x+q2 sin 2x+pz cob Sx+q^ sin Sx+ 

y2=—a2+p+pi cos 2x+qi sin 2x+p2 cos ^x+qt sin 4x+p3 cos 6ic+g3 sin 6a; -f- ^ 

t;3=--.a3+p+pi cos 3x+gi sin Zx+p2 cos Qx-\-q2 sin 6¿c4p3 cos ^x+qz sin 9x+ 



i;n-i=~~an-i+p+pi COS {n—l)x+q\ sin (n--l)x4p2 cos 2{n—\)ic^q2 sin 2(n—\)x 

\p% cos 3(n— l).T+g3 sin ^(n—l)x-\ _ - 

The condition which makes the values of p and q most probable is that 

vl \ v\~\-v\-^vl+ =- minimum. 

and writing instead of vl, v\, etc., their values we have 

(—ar^p^rP\+P2]-^ .__)24 (— aifpH P\ cos x+gi sin x+ )^ 

~f-(-~a2+p4 Pi co8 2x j-gi sin2x-f „)2=:a minimum. 

We know from the differential calculus that when a function is a minimum, its deriv- 
atives with respect to each of the variables are zero. Hence, deriving the former ex- 



^ The quantity p which appears in these equations represents the constant M which we saw com- 
pleted BesseFs formula; it is also to be noted that the equations on page 7 are given here according to 

formula (I). 
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pression with respect to p, p^, q^, etc., and representing these derivative by P, P^, d, 
etc., we obtain the following system — 

P=2(—a0+p+pi+p2+ )+2(— ai+p+pi cos x-{-qi sin x+ ) 

+2(--a2+p+pi cos 2x+qi sin 2x+ ) + =0 

Pi=2(— ao+p+|?i-f392+ )+2 cos xi-^ai+p+pi cos x+qi sin x+ ) 

+2 cos 2x(— a2+P+Pi cos 2x+qi sin 2x+ ) + =0 

Qi=2 sin x{—ai+p+pi co8.a:+gi sin x+ )+2 sin 2x{~a2+p+pi cos 2x+qi sin 2x+ ) + =0 

P2==2(— ao+J>+pi+P2 )+2 cos 2a;(— ai+p+pi cos x+qi sin x+ ) 

+2 cos éx(—a2+p+pi cos 2x+qi sin 2x+ ) + =0 

Q2==2 sin 2x{'-ai+p+pi cos x+qi sin x+ )+2 sin áx(—a2+p+P)^ cos 2ic+gi sin 2x+ ) + =0 



This is the system which must be solved in order to find the values of p and q, which 
will be the most probable ; because in this hypothesis the preceding equations have been 
obtained. 

After dividing by 2 we can write it in the following form : 

iP=— Oo +np+pi 1 +qt sin a; +p2 1 +^2 sin 2a; + =0 

— tti + cos X + sin 2x 

—as + cos 2x + sin 3a; 

-aa + cos 3a; + sin 4x 



JP,= 



— Ctn-l 

— aicosa; 
—as cos 2a; 

— eta cos 3a; 



+cos(ti~l)a; 



+ sin(n--l)a; 



+ cos 2a; 
+ cos 4x 
+ cos 6x 
+ cos 8a; 



+P 



— an.icos(n-l)a; 



iQi=- 



-ai sin X 
-ai sin 2a; 
-aasin Zx 



+P 



ci,i.tsin(w-l)a; 



1 
+ cosa; 
+ cos 2a; 
+ cos 3a; 



-f cos(w~l)a; 



+p^ 



1 
+ cos^ x 
+ cos^ 2x 
+ cos^ 3a; 



+ cos^{n-l)a; 



+ cos2(ti-l)a; 

+g 



sin 2a; 
+ sin 4a; 
+ sin 6a; 
+ sin 8a; 



+ sin 2(n~l)a; 



cos X sin X 

+ cos 2a; sin 2a; 

+ cos 3a; sin 3a; 



+ Pa 



+ cos X cos 2a; 
+ <:os 2a; cos ax 
+ cos 3a; cos 6x 



+ 92 



+ COS (n-1 ) a; cos 2 (n-1 ) x 

sin X +pi 

+ sin 2x 
+- sin Bx 



+ cos (n-l) X sin (n-1) x 

cos X sin 2a; 
-f cos 2a; sin 4a; 
+ cos 3a; sin 6x 



+. 



.=0 



+ COS (n-1) a? sin 2 (n-1) x 



+ sin (n-1) a; 



sin x cos X 
+ sin 2a? cos 2a; 
+ sin 3a; cos Zx 



+qt 



+ sin (n-1 ) X cos (n-1 ) a; 



sm a; 
+ sin' 2a; 
+ sin^ 3a; 



+ sin' (n-1) a; 



+P^ 



sin X cos 2x 
+ sin 2a; cos ix 
+ sin 3a; cos 6x 



+9^ 



+ sin (n-1 ) cos 2 (n-1 ) a; 



sin X sin 2a; 
+ sin 2x sin 4a; 
+ sin 3x sin 6a; 



+ sin (n-1) X sin 2 (n-1) a? 
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This system, according to the notation explained on page 11, may be written in 
the following manner: 

— Sam+np+pi 2 cosmx+qi 2 sinrnx+pz ^cos2mx+q2 2 sin 2ma3+ =0 

—2am cos mx-j-p 2 cos mx+pi 2cos* mx-j-qi 2 cos mx sinmx-i-pz 2 cos ma? cos 2mx+g2 2 cos mx sin 2mx-\~=0 
— 2amsinmx+jt?2 sin7?ta?+/>i ^ sin mar cos ma; +g] 2 sin* ma;+p22 sin ma; cos2ma:+^2 2 sin ma? sin 2ma;+ =0 
— 2amcos2ma;+p2 cos2ma;+i?i2 cos2ma;cosma;+gi2 cos2?na?sinmi»+p22cos^ 2ma; 

+^^2 2 cos 2ma; sin 2ma;+=0 
—2am sin 2mx+p 2 sin 2ma;H-pi 2 sin 2ma; cos ma;+gi 2 sin 2mx sin mx+p2 2 sin 2ma; cos 2mx 

+92 2 sin ' 2mx+=0 



When nx (the product of the number of observations multiplied by the number of 
degrees of the parts into which the different circles of reference are divided) is a mul- 
tiple of 27r, the preceeding system is reduced to the following:^ 

—Mm+np =0 

—Mm cos mx-^p^ S cos^ mx =0 

~2a„,sin mx-^-q^ 2 sin^ mx =0 \ (B) 
— 2a^ cos 2mx+P2 S cos^ 2mx =0 | 
—Mra sin 2mx+q2 2 sin^ 2mx =0 J 



But as will be proved later 



2 cos^ mx=^ 
2 sin^ mx=^ 

2 cos^ 2mx=^ . 

2 sin^ 2mx=^ 

then substituting these values in system (B) we obtain finally 

™-2a,+ni> =0\ ^2am 

~2a„. cos ma:+| pi=0 1 ^" ^ 

Z ^ pi=2a™ cos mx 

--2a, sin mx-{- 1 gi^Of ^ 

^ ) •'• I gi=2a„, sin ma; 

— Sa, cos 2ma:+y ??2=0[ ^ 

„ \ I P2=2a ^ cos 2ma; 

— 2a,„ sin 2mx+^ q2=0 1 2 

/ ^ q2=Mm sin 2ma: 



2 



These are in their general form Bessel's equations which we required ; in our par- 
ticular case of twenty-four observations, ^=12, which reduces these equations to those 

on page 11. It is also proved that the constant p or M, which completes BesseFs for- 
mula, is the arithmetical mean of the observations ; thus we have in our case : 

ao+ai+a2H- 4 023 

It only remains to prove the change from (A) to (B) and the equations (C) . 
For the first it is sufficient to prove that when nx is a multiple of 2w 

2 sin mx = 0, 2 cos mx == 0, 2 sin 2mx = 0, 2 cos 2mx = 0, etc. 
* This will be proved later on. 
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Let US take the well known trigonometrical formula: 

2 ain z sin i,=cos(.-i,) -cos («+!/) /. sin ,^ co8(.-i0-^cy(.4-y) ^^ 
By a series of applications of Formula (1) in which we write x, 2x, Sx 



(n~-l)a; instead of z and| instead of y, we have: 

X i^JT 

COSg—COS^ 

sin x= 

2sin| 

3a: 5x 

cos^ -cos-g- 

sin 2x= =— ^ 

2 sin| 

^ * 5x 7x 

cos — cos-p- 

sin doc= -^ — -- — ^ 

2 sin| 



{2n^B)x i2nS)x 



sin(w-2)a:= 



2sin| 
cosi?^->i-cosl25^^ 



sin(n~l)x= 
and taking the sum we obtain: 



2 sin| 



sin ^4-sin 2a:+sin 3x4- + sin (n-l)a;=2 sin mar = 



cos |~-cos nx cos|— sin nx sin- 

2 sin| 

if nx is a multiple of 2^^ we know that 

cos nx=l; sin na:=0 
and hence 

cos|— cos| 
5 sin mx=^ '— — »=0 

2 sin| 
2 



x (2n~l)a; 

cos s— cos^'^ — ^ ^ 



2 sin? 



In like manner it is proved that 2 sin 2mx = 0, S sin 3mar = 0, etc., so for this last 
will write in Form 
of y, which will give: 



we will write in Formula (1) 3x, Qx, 9x, , in~-l) Sx instead of z, and ^ instead 



cosf^cosl^M cos|^cos3tia:cos3?^Bin3nxsint 

2 sin Bmx^--——-———^—-- ..^=:^^^^A.____ ^¿_ ^ ¿_ 

2sin?| 2sin§ 

2 2 

and as Bnx is a multiple of 2ir, it follows that 2 sin 3wai = 0, 
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To prove that 2 cos mx = 0, S cos 2mx = 0, we make use of the trigonometrical for- 
mula 

2 cos z sm y=sin{z^y) —siniz-^y) .\ cos z= — 2 sin y — ^^ 

X 

and writing 0, x, 2x, (n—l)x instead of z, and ^ instead of y, we have 

cos 0= ;;;; — 

2sin| 
2 

cos x= — 

2 sin? 

sin|^-sm| 

cos2 ÍC =- -~ 

2sin| 

sm^~sm — 

cos 3a;= TT— ~ 

2sin^ 



sm A_-~ — ^ sm ^i — ^ ' 



cos(n— 2)a;=- 



2sin^ 



. i2n^\)x .^ (2n-~3)a: 
sm -'^ — ^-^^ sm ^^ — ^-J-^ 

cos (n— 1 ) a: == — 



2* it/ 
smg 

Taking the sum we obtain : 

ain , „^^^ 

2' 2 



sm^+sm "^ ^ 



1+cos a::+cos 2a;+ 4-cos(n-~l)a: = S cos wa;=- 



2 sin 5 



sin|+sin nx cos |— cos nx sin| 
¿a ¿á ¿á 

2 sin I 
and remembering that nx=%n^ it follows that: 

• U/ * X 

2 cos wa;= — = ==0 

Zsini 

127182 8 
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This having been proved, it is readily seen that the other terms of system (A) of 
the form 

S cos mx sin mx; % sin mx cos 2mx; S sin mx sin 2mx; 5 cos ma; sin 2mx; 5 cos mx cos 2mx 

are also reduced to zero ; because as the product of the sine and cosine of a same arc 
is equal to the half of the sine of the double arc, we obtain : 

2 cos mx sin mx = ^ 5 sin 2mx = 

the same can be proved for the other four expressions by making use of the well-known 
trigonometrical formulas 

sin ¿r cos |/= i sin (a:-f 2/) +1 sin(ic— Í/) 
cos xsiny= J sin {x-\-y) — i sin (o;—]/) 
cos x cos tj—ico^ix-i-y) +J cosix—y) 
sin X sin y = | cos (x—y) — ^ cos (x-^^y) 

for writing mx and 2mx for x and y, we have : 

5 sin mx cos 2mx^i% sin dmx—i:s, sin mx 
5 cos mx sin 2ma:=|2 sin Bmx-j-^'S^ sin wo; 
2 cos mx cos 2ma:=|2 cos 3ma:;+42 cos mx 
2 sin mo; sin 2mír=Í2 cos mo:— ^2 cos 3mx 

and as was proved before, the second members are zero. 

We now proceed to prove the equations (C) ; a simple trigonometrical expansion 
gives us that: 

2 cos^ x == l+cos2i^ 
and according to this we have : 

2cos2 = 1+cosO 

2cos^a; = 1+cos 2x 
2 cos^ 2x = 1-f cos 4x 
2 cos^ Sx = 1-f cos Qx 



2cos2(n— l)a; = l+cos2(n— l)a? 
taking the sum we get 

2 2 cos^ ma;=n+ (1+cos 2a?+cos 4a?+cos 6x+ ..-feos 2(n-l)x = n+2 cos 2mx 

but, as we proved before 2 cos 2mx = 0, hence 

2cós2ma:=^ 

In an analogous manner this same value can be proved for any term of the form 
2 cos® mkx. 

To prove that 2 sin^ ma; = | we make use of the known equation 

sin^ X = 1— cos=^ X 
and writing 

sin^ = 1— cos^ 

sin^ X = 1— cos^ X 
sin^ 2x = l-^-cos^ 2x 
sin^ Zx = 1— cos^ So: 



8inM^"-l)a; = 1— cosH^-~l)ir 
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and adding, we have 



S sin^ mx == n— S cos^ mx 



n 



but we have just proved that % cos^ mx = ^, therefore 



5 sin^ mx = ^— H==f 

We do not think it necessary to delay any longer over this work and give the whole 
of BesseFs method ; those who may be interested in it may consult BesseFs own expla- 
nation, translated into English by Mr. Robert H. Scott, and published as an appendix to 
the Quarterly Weather Report of the Meteorological Office, London, 1870. 

Table I. — Mean hourly values of atmospheric pressure at Manila for each month (1890-1909), 



Hour. 



0. 

1- 

2- 

3. 

4_ 

5 

6- 

7. 

8. 

9_ 
10. 
11. 
12. 
13. 
14. 
15. 
16- 
17. 
18. 
19. 
20. 



28. 
24- 



January. 



61.382 
61.092 
60.722 
60.464 
60.516 
60.718 
61.069 
61.546 
62. 028 
62.264 
62.128 
61. 701 
61.057 
60.282 
59.656 
59.360 
59. 462 
59. 814 
60.228 
60.698 
61.158 
61.461 
61.544 
61.496 
61.365 



Febru- 
ary. 



61.704 
61.404 
61.072 
60.854 
60.833 
61.047 
61.470 
61.943 
62.342 
62.630 
62.568 
62. 120 
61.494 
60.712 
60.058 
59.650 
59.666 
59.954 
60.362 
60.798 
61.274 
61.654 
61.839 
61. 852 
61.717 



March. 



60.924 
60.556 
60.214 
59. 982 
59.989 
60.204 
60.661 
61. 172 
61.564 
61.762 
61.630 
61.250 
60.620 
59.827 
59. 134 
68.672 
58.620 
58.867 
59.266 
59.749 
60.261 
60.742 
61.012 
61.078 
60.898 



April. 



59.^8 
59.624 
59.228 
59.029 
59.071 
59.346 
59.704 
60.193 
60.584 
60.693 
60.569 
60.143 
59.562 
68.732 
58.033 
67.615 
57.490 
67.700 
68.204 
68.762 
69.266 
59.753 
60.124 
60.170 
59.902 



May, 



June. 



July. 



58.872 i 

68.601 

58.187 

57.982 

57.991 

58.168 

68.486 

58.885 

59.166 

59.274 

59.142 

68.816 

58.333 

57. 686 

57. 118 

56.771 

66.656 

56.890 

57.402 

67.917 

58.397 

68.835 

59.144 

59. 157 

68.862 



58. 520 
58.193 
57. 789 
67.620 
57.556 
57.649 
67.950 
58,258 
58.446 
58.539 
58.456 
58.188 
57.824 
57.344 
56.910 
56.576 
56.496 
66.706 
57. 130 
57.619 
58.078 
58.474 
58.755 
58.811 
58.498 



57.880 
67.468 
67.128 
56. 936 
56.870 
66.944 
67,222 
67.546 
57. 762 
57.869 
57.819 
57 622 
67.283 
56.846 
66.444 
56.169 
56.066 
56.232 
56.598 
67.062 
67.488 
57.861 
68.158 
58.172 
57.868 



Augrust. 



57.834 
67.459 
67.102 
66.868 
56.820 
56.926 
57. 191 
57. 518 
57.830 
57.980 
57.918 
57.734 
57,354 
56.828 
66.420 
66.163 
56.087 
66.264 
56.640 
57.087 
57.528 
57.924 
58.218 
58.177 
67.852 



Septem- 
ber. 



67.934 
57.549 
67.190 
56.944 
56.885 
57.046 
57.306 
57.705 
58.117 
58.246 
58.182 
57.930 
57.428 
56.798 
56.378 
56.087 
56.128 
56.368 
56.788 
57.240 
57.688 
68.101 
68.294 
68.180 
67.921 



October. 



68.824 
68.408 
68.080 
57.883 
57.842 
68.030 
68.476 
58.915 
59.248 
59.462 
59.366 
58.928 
58.349 
57.732 
57.204 
56.944 
57.078 
57.391 
57.732 
58.196 
68.728 
59.071 
59.208 
59. 180 
68,867 



Novem- 
ber. 



69.684 
69.380 
58.970 
68.686 
68.726 
68.953 
59.388 
69.836 
60.218 
60.423 
60.274 
59.877 
59.258 
58.572 
58.071 
57.879 
57.993 
68.344 
58.708 
59.218 
59.698 
69.964 
60.054 
59.958 
59.736 



Decem- 
ber. 



60.880 
60.690 
60.162 
59.930 
69.967 
60.176 
60.512 
60.966 
61.467 
61.642 
61.480 
61.084 
60.485 
59.735 
59. 198 
58.952 
59.082 
59.480 
59.882 
60.369 
60.811 
61.094 
61.166 
61.059 
60.900 



Table II. — Same as Table I corrected for noncyclic variation. 



Hour. 



0.... 

1.... 

2,.„. 

8..-. 

4.... 

5.... 

6...- 

7.... 

8.... 

9..._ 
10.... 
U-... 
12.... 
18.... 
14.... 
15.... 
16.— 
17-._. 
18.... 
19..- 
^.... 
21.... 
22... 
28.... 
24.... 



January. 



6L374 
61.084 
60.715 
60.458 
60.510 
60.713 
61.066 
6L642 
62.026 
62.262 
62.127 
6L700 
6L057 
60.288 
69.657 
69.862 
69.466 
69.818 
60.232 
60.70» 
61.164 
6L467 
6L551 
61.604 
61.874 



60.910 



Feb- 
ruary. 



61.710 
61.410 
6L077 
60.869 
60.887 
6L061 
6L473 
61.946 
62.844 



62.120 
61.494 
60-712 
60.(^7 
69.648 
69.664 
69.961 
60.849 
60.794 
6L270 
6L649 
6L884 
6L846 
6L710 



6L22Í 



March. 



60.911 
60.644 
60.208 
69.972 
59.980 
60.196 
60.666 
6L167 
61.661 
6L769 
6L628 
6L249 
60.620 
69.828 
69.186 

il:^ 

68.872 
69.272 
59.757 
60.270 
60.752 
6L028 
6L090 
60. 9U 



April. 



69.608 
69.218 
59.015 
69.059 
69.336 
59.695 
60.186 
60.678 
60.688 
60.666 
60.141 
59.662 
68.784 
68.036 
67.620 
57.496 
67.707 
68.218 
68.772 
69.278 
69.767 
60.189 
60.186 
59.^0 



69.818 



May. 



68.862 
68.692 
68.179 
67.974 
67.984 
68.162 
68.481 
68.881 
69.164 
59.272 
69.140 
68.816 
68.388 
57.687 
67.120 
66.774 
66.668 
66.894 
67.407 
67.928 
58.404 
68.848 
69. 162 
69.166 



58.244 



June. 



68.609 
58.183 
57.780 
57.612 
57.549 
57.648 
67.944 
58.268 
68.443 
58.686 
58.464 
68.187 
57.^4 
57.346 
66.912 
66.679 
66.^9 
66.711 
67.186 
67.626 
68.0S6 
68.482 
68.764 
68.821 
68.609 



57.828 



July. 



67.874 
67.463 

66.866 
66.940 
67.219 
67.644 
57.760 
67.867 
67.818 
67.^1 
67.288 
66.847 
66.446 
66.171 
56.068 
66.284 
66.601 
67.066 
67.492 
57.866 
68.168 
68.177 
67.874 



67.^7 



AUfiTUBt. 



67.848 
67.467 
57. 110 
66.876 
66.826 
66.981 
67.186 
67.622 
67.888 
67.982 
67.919 
67.786 
67.354 
66.827 
66.419 
66.161 
66.084 
66.260 
66.686 
67.082 
67.622 
67.917 
68.210 
68.169 
67.848 



67.244 



Septem- 
ber. 



67.928 
67.648 
67.185 
66.989 
66.881 
67.042 
67. 8(^ 
57.702 
68.116 
68.244 
68.181 
67.929 
67.428 
66.799 
66.879 
66.089 
66.180 
66.871 
66.791 
67.244 
67.692 
68.106 
68.299 
68.186 
67.928 



67.364 



October. 



58.846 
58.428 
58.098 
57.899 
57.856 
58.043 
68.487 
68.924 
69.266 
69.467 
69.870 
58.980 
68.849 
67.780 
67.200 
66.989 
67.071 
67.882 
67.721 
68.188 
68.714 
69.066 
69.190 
69. 110 
68.846 



68.343 



Novem- 
ber. 



69.710 
69.864 
68.992 
68.706 
68.748 
58.968 
69.401 
69.847 
60.227 
60.430 
60.278 
69.879 
69.268 
68.670 
68.067 
67.872 
67.^4 
68.888 
68.696 
69.208 
69.681 
69.944 
60.032 
69.934 
69.710 



69.264 



Decem- 
ber. 



60.890 
60.699 
60.170 
69.938 
69.974 
60.182 
60.617 
60.970 
61.469 
6L644 
61.482 
61.085 
60.486 
69.784 
^.196 
68.960 
69.080 
69.476 
69.877 
60.368 
60.804 
61.086 
61.148 
61.060 
60.890 



60.424 



16 HARMONIC FORMULA OP FOURIIR AND BESSEL, 

Tamim III. — Diurnal variation of atmoepheHc pressure at Manila {1890-1909), corrected for noncyclic 

variation (unit, 0,01), 



Hour. 


January. 


Febru- 
ary. 


March. 


April. 


May. 


June. 


July. 


AURUBt. 


Septan- 
ber. 


October. 


Novem- 
ber. 


Decem- 
ber. 


„ 


+ ^ 
+ 17 

— 20 

— 45 

— 40 
-20 
+ 16 
+ 68 
+112 
+136 
+1^ 
+ 79 
+ 16 

— 63 

-m 
-1^ 

—144 
-109 

— 68 
-21 

+ 25 
+ 66 
+ 64 
+ 69 


+ 49 
+ 19 

- 14 

- 36 
-38 

- 17 
+ 25 
+ 72 
+112 
+141 
+186 
+ 90 
+ 27 

- 61 
—116 
-157 
~-166 

-m 

- 48 

+ 5 
+ 43 
+ 61 
+ 62 


+ 59 
-22 

- 12 

- 35 

- 34 

- 18 
+ 33 
+ 84 
+124 
+144 
+ 130 
+ 98 
+ 30 

- 50 
-119 
-166 
-170 
—145 
—106 

- 57 

- 5 
+ 48 
+ 70 
+ 77 


+ 61 
+ 80 

- 10 

- 30 

- 25 
+ 2 
+ 38 
+ 87 
+126 
+138 
+124 
+ 88 
+ 26 
-68 
—128 
—169 
-182 
—161 
—110 

- 64 

- 4 
+ 45 
+ 83 
+ 87 


+ 62 
+ 36 

— 6 

— 27 
~ 26 

— 8 
+ 24 
+ 64 
+ 92 
+108 
+ 90 
+ 67 
+ 9 
-66 
—112 
-147 
—169 
—135 
-84 
-32 
+ 16 
+ 60 
+ 91 
+ 92 


+ 68 
+ 36 

— 6 
~ 22 
-28 

— 18 
+ 12 
+ 42 
+ 62 
+ 71 
+ 63 
+ 36 



— 48 

— 92 
—126 
—133 
-112 
-69 
-20 
+ 26 
+ 65 
+ 94 
+ 99 


+ 66 
+ 24 

- 10 

- 30 
-36 
-29 

- 1 
+ 32 
+ 63 
+ 64 
+ 59 
+ 39 
+ 6 

- 38 
-78 
-106 
-116 
-99 
-63 

- 16 
+ 26 
+ 64 
+ 94 
-i- 95 


+ 60 
+ 22 

— 13 
-37 

— 42 

— 31 

— 6 
+ 28 
+ 59 
+ 74 
+ 68 
+ 49 
+ 11 

— 42 

— 82 
-108 
—116 
-99 
-61 

— 16 
^- 28 
+ 67 
+ 97 
+ 92 


+ 57 
+ 19 

— 17 

— 42 

— 47 

— 31 

— 5 
+ 35 
+ 76 
+ 89 
+ 83 
+ 68 
+ 7 
-56 

— 98 
—126 
-122 

— 98 

— 56 

— 11 
+ 34 
+ 75 
+ 94 
-4- 83 


+ 50 
+ 8 
-26 

— 44 

— 49 

-30 
+ 14 
+ 68 
+ 91 
+111 
+103 
+ 59 
+ 5 

— 61 
—114 
-140 
-127 
-96 

— 62 

— 16 
+ 37 
+ 71 
+ 85 
+ 77 


+ 46 
+ 10 

— 26 

— 55 

— 51 
-29 
+ 15 
+ 69 
+ 97 
+118 
+102 
+ ^ 


-68 
—119 
-138 
-127 
-92 
-56 

— 5 
+ 43 
+ 69 
+ 78 
+ 68 


+ 47 
+ 18 
-25 
-49 

— 45 

— 24 
+ 9 
+ 55 
+104 
+122 
+106 
+ 66 
+ 6 

— 69 
—123 
-147 
—134 

— 96 

— 55 

— 6 
+ 38 
+ 66 
+ 72 
+ 63 


1 .,_, 

2 


8 


4 


6„. _ 


6 ._ 


7 


8 


9 „._ 


10 


11 , 


12 

18 


14,— . 


16,. 


16. 


n....,-.. ... 


18 


19.. 


20.. ■ 


21..... 


n 


28.............. 





Tabl® IV. — Constants of the Harmonic Formula (I) computed from Table III. 



Month, 


Pv 


91- 


P2' 


«2- 


Ps- 


<?8- 


Pa- 


Qe 


January.... . 


-f 0.1341 

+ .(mi 

+ .1208 
+ .2022 
+ .2980 
+ .8622 
+ .2896 
+ .2679 
+ .2476 
+ .1994 
+ .1841 
+ .1898 


+0.5476 
+ .6696 
+ .7496 
+ .7948 
+ .5711 
+ .3860 
+ .2829 
+ .2788 
+ .8162 
+ .4065 
+ .3970 
+ .4364 


+0.2696 
+ .8488 
+ .3949 
+ .3888 
+ .8396 
+ .8269 
+ .3871 
+ .8458 
+ .8061 
+ .2466 
+ .2080 
+ .2169 


II 1! 1 1 1 1 1 M ¿ 


++++++++++++ 


++++ 1 1 1 1 ++++ 


+0.0196 
+ .0268 
+ .0362 
+ .0883 
+ .0313 
+ .0800 
+ .0192 
+ .0088 
+ .0041 
+ .0066 
+ .0125 
+ .0188 


—0.0166 

— .0187 
-.0161 

— .0158 

— .0007 


— .0101 

— .0246 

— .0274 

— .0817 

— .0144 

— .0122 


February 


March 


April 


Ky :.::::;::::::::::::::: 


June 


July..-.. I. 


Auffutt -_, . 


September 


October... 


November...- , 


December 





Table V. — Amplitudes and phase angles of Fourier *'waves** computed from Table IV, using L, M, T. 



Month. 


«1- 


t*2- 


1*8. 


1*4. 


Uv 


t/2. 


u,. 


V^, 


Januttry ^-.-.. .. . . 


0.6688 
.6607 
.7692 
.^02 
.6442 
.6231 
.4047 
.3794 
.4007 
.4^ 
.4376 
,4767 


0.9782 
.9^6 
1.0816 
1.0414 
.9161 
.79^ 
.7604 
.8062 
.8802 
.9409 
.9467 
.9657 


0.1850 
.1100 
.0647 
.0271 

.oom 

.0874 
.0488 
.0269 
.0412 
.0670 
.0974 
.1271 


0.0267 
.0319 
.0392 
.0416 
.0813 
.0300 
.0217 
.0269 
.0278 
.0324 
.0600 
.0707 


/ 

13 46 
3 20 
9 9 

14 17 

27 33 
42 20 
45 39 
42 49 
38 9 
26 8 
24 53 

28 27 


o / 
164 37 
159 32 

167 30 

168 6 

168 14 

166 44 
153 41 
164 36 

169 43 
164 49 

167 36 
166 53 


o f 

11 1 

12 14 
37 40 
53 8 

123 42 

164 49 

159 39 

141 17 

28 6 

37 7 

30 68 

15 26 


o / 
180 17 
126 68 
112 34 
112 25 
91 22 
90 
117 51 
161 16 
171 26 
168 12 
167 69 
164 37 


February -.-.,......... 


SSSuzi:::::::::::::::::::;::: ;:::;": 


April...-,...-., .. .. . 


^!::..;.:.:::::;:::;;:;::::::;;:;: ::;;;: ::'::;: 


June. .-.,--.. 


July 


AuifUit. ,... „..„.. , 




^U»..,:;;;:::::::::;;:::r "■":::"■ 




December. ,. ... . 





HARMONIC FORMULA OP FOURIER AND BESSEL. 
Table VI. 
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JANUARY. 


FEBRUARY. 


Hour. 


Fourier "waves." 


Result- 
ant. 


Depart- 
ures 
from ob- 
served 
values. 


Hour. 


Fourier "waves." 


Result- 
ant. 


Depart- 
ures 
fromol)- 
served 
values. 


First. 


Second. 


Third. 


Fourth. 


First. 


Second. 


Third. 


Fourth. 





+0.134 
+ .271 
+ .390 
+ .482 
+ .541 
+ .564 
+ .548 
+ .495 
+ .407 
+ .293 
+ .158 
+ .012 
-.134 

— .271 

— .390 

— .482 

— .541 

— .564 

— .548 

— .495 

— .407 

— .293 

— .158 

— .012 


+0.259 

— .246 

— .687 
-.943 

— .947 

— .696 

— .259 
+ .246 
+ .687 
+ .948 
+ .947 
+ .696 
+ .259 

— .246 

— .687 

— .948 

— .947 

— .696 

— .259 
+ .246 
+ .687 
+ .943 
+ .947 
+ .696 


40.026 
+ .112 
+ .188 
+ .075 

— .026 
-.112 
-.133 

- .075 
+ .026 


+0,020 

— .004 

— .024 

— .020 
+ .004 
+ .024 
+ .020 


+0.439 
+ .133 

— .188 

- .406 

- ,m 

- .220 
+ . 176 
+ .662 
+1.096 
+1.828 
+1.242 
+ .807 
+ .119 
-.633 
-1.234 
-1.520 
-1.458 
-L124 

— .664 

— .178 
+ .230 
+ .518 
+ .660 
+ .633 


-0.02 

— .04 
+ .01 
+ .04 
-.03 

— .02 
+ .02 
+ .03 

— .02 

— .02 
+ .02 
+ .02 
-.03 


+ .02 
+ .08 

— .02 

— .03 
+ .08 
+ .03 

— .02 

— .04 
+ .02 
+ .04 





+0.038 
+ .208 
+ .864 
+ .494 
+ .591 
+ .647 
+ .660 
+ .627 
+ .668 
+ .440 
+ .297 
+ .184 
- .088 


+0.344 
-.163 

— .^7 
-.922 

- .970 

- .769 

— .344 


+0.028 
+ .098 
+ .107 
+ .060 
— .023 


+0.026 
-.008 

— .029 

— .026 


+0.481 
+ ,185 
-.185 

— .394 

— .399 
-.176 
+ ,236 
+ .727 
+1.174 
+1.429 
+1.877 
+ .982 
+ ,809 

— .467 
-1.127 
—1.502 
-1.636 
—1.284 

— .871 
-.407 
+ .022 
+ .868 
+ .669 
+ .694 


-0.06 

- .06 

- .04 
-.03 
-.02 

- .01 

- .01 
+ .01 
+ .05 
+ .02 
+ .03 
+ .08 
+ .04 
+ .04 
+ .08 
+ .07 
+ .02 
-.01 


+ .Oi 
-.03 

- .07 
-.04 

- .08 


1 . - 


1 


2 — . - 


2 


8 


8 


4 


4 


6 


5 




6 


6 






7 


7 






8 




8 








9 . . 




9 








10- . 






10 








11- - 






11 








12 - . 






12 








13 






13 








14 






14. 










16— 

16 






15 










16 










17 






17 










18 






18 










19 






19 ,.- 










20 






20 










21 - 






21 










22 - 






22 










23 






28 



























MARCK. 


APRIL. 


-.- 


+0.121 
+ .310 
+ .479 
+ .616 
+ .710 
+ .755 
+ .749 
+ .692 
+ .589 
+ .445 
+ .270 
+ .077 
-.121 


+0.396 

— .135 

- .628 

— .964 
—1.023 

- .818 

- .396 


+0.040 
+ .064 
+ .061 
+ .008 
- .040 


+0.036 
+ .006 

— .081 

— .086 


+0.592 
+ .234 

- .129 

- .366 

- .358 

- .096 
+ .889 
+ .824 
+1.226 
+1.427 
+1.389 
+ .984 
+ .270 

- .604 
-1.189 
-L614 
-1.698 
-1.478 
-1.067 

- .644 

- .082 
+ .409 
+ .697 
+ .764 



+0.01 

- .01 

- .02 

- .02 
+ .03 
+ .01 

- .02 

- ,01 

- .01 
+ .04 



- .03 




+ .04 



- .08 

- .01 
+ .08 
+ .02 

- .02 


.01 





+0.203 
+ .401 
+ .672 
+ .706 
+ .790 
+ .820 
+ .795 
+ .716 
+ .687 
+ .419 
+ .222 
+ .011 
-.203 


+0.388 

- .147 

- .642 

- .966 
-1.081 

- .819 

- .888 
— -- 


+0.022 
+ .027 
+ .016 
- .004 
-.022 


+0.089 
+ .006 
-.088 
— .039 


+0.652 
+ .287 
-.067 

- .804 

- .269 
+ .007 
+ .480 
+ .872 
+1.218 
+1.878 
+1.268 
+ .869 
+ .202 

- .569 
-L268 
-1.706 
-1.806 
-1.579 
-1.128 

- .566 


+ .481 
+ .787 
+ .846 


+. 04 

- .01 
+ ,08 



- ,02 

- ,01 
+ ,03 



- .04 

- .01 
+ .02 
+ .08 

- .06 
+ .01 
+ .02 

- .02 
+ ,01 
+ .08 

- .08 

- .08 
+ .04 
+ .08 

- .04 
•- .01 


1 


1 


2 


2 


3 


8 


4 


4 


5 




5 




6 






6 






7 






7 






g 








8 - 








9 








9 








10 








10 








11 








11 








12 


""" i 




12 








13 








13 








14 










14 










15 










16 










16 










16 










17 










17 










18 










18 










19 










19 










20 










20 










21 










21 










22 




i ' 




22 










23 




1" \..,l 




28 














1 














MAY. 






JUNE. 







+0.298 
+ .485 
+ .544 

t:Vá 

+ .629 
+ .571 
+ .475 
+ .846 
+ .198 
+ .0^ 

- .140 

- .298 


+0.840 

- .181 
-.567 

- .861 

- .907 

- .720 

- .340 


+0.007 
+ .002 
-.005 

- .009 

- .007 


+0.081 
+ .015 

- .016 

— .081 


+0.676 
+ .321 

— .044 

— .277 
-.286 

- .077 
+ .267 
+ .630 
+ .904 
+1.016 
+ .916 
+ .687 
+ .066 

- .668 
— L122 


+0.06 

- .08 
+ .02 

— .01 

— .08 


+ .(» 

— .01 

— .02 

— .01 
+ .02 
+ .02 

- .02 
+ .01 



- .02 
+ .08 
+ .02 

- .04 

- .02 
4 .04 
+ .08 

— .04 

- .08 





+0.852 
+ .440 
+ .498 
+ .622 
+ .611 
+ .464 
+ .886 
4 .282 
+ .168 
+ .026 

— .112 

- .240 
-.862 


+0.826 

- .079 

- .468 

- .728 

- .789 

- .644 

- .826 


+0.010 
-.018 
-.086 

- .082 

- .010 


+0.080 
+ .016 

- .016 

- .080 


+0.718 
+ .868 

— .016 
-.268 
-.808 

- .177 
+ .126 

+ .^(m 

+ .616 
+ .700 
+ .626 
+ .887 
-.006 
-.486 
-.940 
—1.248 
-1.806 
-1.111 
-.718 
-.220 
+ .280 
+ .686 
+ .922 
+ .931 


+0.04 


+ .08 
-.04 

- .02 


+ .01 

- .01 


- .01 


+ .08 

- .01 

- .01 

- .02 
+ .01 
+ .02 
+ .01 

- .08 
^ .02 
+ .02 
+ .04 
-.02 

- .06 


1 


l._. 


2 


2 


8 

4 


8 

4 


6 




6-. 




6 - 






6 

7 






7 










8 








8 .- 








9 








9 








10 








10-.—. 








11 








11 








12 








12 








13 








18 








14 'i 










14 










16 1 








—1.487 


16 










16 








-L668 
—1.881 

— .885 

- ,888 
+ .198 
+ .625 
+ .869 
+ .^6 


16. 










17 








17 










xg 








18 










19 








19 











20 








20 










i 21 


'"rT"" :::::::::!::::::::: 


21 










1 22 i 


i., 


22 










1 2a 1 




23 - 










! ' 












1 





38 



HARMONIC FORMULA OF FOURIER AND BESSEL. 
Table VI — Continued. 
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INTRODUCCIÓN- 



Juzgo que el presente folleto "The Harmonic Formula of Fourier and Bessel, and 
Its Aplication to the Study of the Diurnal Variation of the Atmospheric Pressure 
in Manila During the Period 1890-1909'' preparado por el P. Antonio Galán, S, J., 
agregado a la dirección de este Observatorio de Manila, ha de ser bien recibido por 
los meteorólogos, los cuales son de parecer que semejante discusión o estudio debería 
hacerse en todos los observatorios principales, en donde se haya acumulado material 
suficiente en regular número de años, para ser tratado matemáticamente con proba- 
bilidades de resultado científico aceptable. 

Ojalá que pronto, como espero, se pueda encontrar tiempo y oportunidad, para 
tratar de semejante manera como ha hecho ahora el P. Galán con los datos baromé- 
tricos, los reunidos en este Observatorio sobre la temperatura, humedad y otros elementos 
de meteorología periódicos. 

Ojalá también que la estancia en el Observatorio de Manila del P, Galán, destinado 
a continuar la benemérita labor de los PP, Vines, Gangoiti y Gutiérrez Lanza del 
Observatorio de Belén de la Habana, le sea tan provechosa como para mí lo fué mi 
permanencia en la Habana los primeros meses del año 1892. 

José Algué, 
Director de la Oficina Meteorológica de FiUpinas. 

Observatorio de Manila, 19 de Junio de 191 A. 
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FÓRMULA ARMÓNICA DE FOURIER Y DE BESSEL 



La fórmula armónica de Fourier y de Bessel suele presentarse en una de las dos 
forman siguientes: 

(I) y=M-}-p^ eos x-\-q^ sen x-^p^ eos 20:4-^2 sen 2x-\~p.^ eos Sx-{-q^ sen 3a; + 

(II) y=M-^u^ sen (Í7i +^0+^2 sen ([7o+2x) +1^3 sen {U^-{-Sx) -}- 

Fácilmente se ve que estas dos fórmulas son una misma bajo distinta forma, pues 
desarrollando la (II) se tiene: 

y=M'\-u^ sen U^ eos x-{-u^ eos U ^ sen x~\--ii^ sen Un eos 2x-^u^ eos U^ sen 2x-\- 

y haciendo en esta expresión 

u^ sen Ui=p^ Uy eos U^^^q^ 

Uo sen U2=P2 ^h eos U^^q'. 

etc. etc. 
la (II) queda reducida a la (I). 

Las aplicaciones de esta fórmula son muy numerosas e importantes. 

El teorema de Fourier (dicen Thomson y Tait) además de ser uno de los más bellos resultados 
del análisis moderno, puede decirse que es instrumento indispensable para abordar casi todas las 
recónditas cuestiones de la física moderna. Decir que las vibraciones sonoras, la propagación de las 
señales eléctricas a lo largo del hilo del telégrafo y la conducción del calor por la corteza terrestre 
son cuestiones en su generalidad imposibles de discutir sin su auxilio, es dar únicamente una débil 
idea de su importancia. 

En el presente trabajo nos limitaremos a la aplicación que de la fórmula se hace 
en meteorología ; procuraremos declarar lo que cada uno de sus elementos representa, 
y haremos algunas indicaciones sobre su utilidad y el modo práctico de servirse de ella.' 

La mayor parte de los fenómenos que se estudian en meteorología son periódicos ; 
las curvas que los representan son, por lo tanto, curvas periódicas, es decir, curvas que 
representan el movimiento de un punto, que recorre un número indeterminado de veces 
el mismo camino en el mismo tiempo. 

Análogamente a lo que Helmholtz demuestra en acústica, que un mismo sonido 
compuesto se resuelve siempre en los mismos sonidos elementales, probó Fourier que 
cualquiera curva periódica se puede descomponer en una serie de curvas armónicas de 
períodos, 1, |, |, etc., de la curva dada; y que sólo es posible ima comhinacióyi de estas 
curvas elementales, para reproducir una curva determinada. 

Prescindiendo por ahora de la utilidad práctica que esta descomposición de las 
curvas pueda tener en meteorología, vamos a ver cómo la fórmula que estudiamos 
no es otra cosa que la expresión analítica de las curvas periódicas en términos de sus 
componentes armónicas. 

Dícese curva armónica la que representa gráficamente el movimiento armónico de 
un punto. 

^ Uno de los mejores trabajos prácticos que conocemos sobre el empleo de las series de Fourier en 
física matemática es el de W. E. Byerly "An Elementary Treatise on Fourier's Series and Spherical 
Harmonics" publicado por Ginn and Co., New York. 

La obra original del mismo Fourier "Théorie Analytique de la Chaleur" es quizá de tanta 
actualidad como cualquiera otra publicación moderna sobre esta materia. Puede consultarse también 
Mellor '*Higher Mathematics for Students of Chemistry and Physics," Longmans, Green and Co., 
39 Paternoster Row, London. 
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FÓRMULA ARMÓNICA DE FOURIER Y DE BESSEL. 





^ií 1 



Movimiento armónico es aquel a que estaría sometida 
la proyección sobre un diámetro de un punto que recorriera 
la circunferencia con movimiento uniforme. Esta circun- 
ferencia recibe el nombre de circulo de referencia. 

Si el punto M (fig. 1) recorre la circunferencia con mo- 
vimiento uniforme, el movimiento 
de su proyección P sobre el diáme- 
tro BD es armónico; mientras M 
describe una vez la circunferencia 
en el sentido ABCDA, al punto 
P recorre el diámetro en el sentido 
OBODO. 

En gracia de aquellos que no 
estén familiarizados con el estudio de esta clase de movimientos, 
permítasenos recordar algunas definiciones, que hemos de 
necesitar. 

Período es el tiempo que emplea el punto M en recorrer la 

circunferencia, o el punto P en su doble recorrido del diámetro. 

Elongación es la distancia variable del centro O al punto P. 

Amplitud es la elongación máxima, o distancia máxima 

a que puede separarse de O el punto P en su movimiento ; como 

se vé, es el radio del círculo de referencia. 

Los puntos O y A son el origen de los movimientos de 
los puntos P y M respectivamente. El origen de los tiem- 
pos es arbitrario; puede tomarse en cualquier punto de la 
circunferencia. 

Fase es la distancia angular del punto 
del movimiento. Tratándose del punto que 
el diámetro, hay que referirle al círculo de 



móvil al origen 
se mueve sobre 
referencia para 

determinar sus fases ; así en la figura 1, si el punto se encuentra 
en P, su fase está dada por los grados del arco MA o del ángulo 
MOA. 

Época de una fase determinada es el instante en que ocurre 
esa fase; así, en nuestro caso, si el arco MA es de 60"" y han 
transcurrido dos horas desde que el punto se encontraba en 
el origen de los tiempos diremos que las 2*' es la época de la 
fase de 60^ 

Fase inicial es la distancia angular del origen de los tiempos 
al origen de los movimientos. En nuestra figura, si tomamos 
al punto N como origen de los tiempos, la fase inicial del 
movimiento está dada por el arco NA. 

Supongamos, para fijar las ideas, que el punto M emplea 
un día en recorrer la circunferencia^ — es decir, que el período 
de su movimiento es veinticuatro horas. A partir de N, origen 
de los tiempos (fig. 2), dividamos la circunferencia en 24 partes 
iguales ; y rectifiquémosla sobre la recta NX; esta queda dividida 
en 24 partes iguales, cada una de las cuales representa la 
longitud de un arco de IS"", o en tiempo, una hora. Hagamos 
que el círculo de referencia se desplace sobre la recta NX 
paralelamente al diámetro vertical, con un movimiento igual 
al del punto que la recorre. Al fin de cada hora encontraremos 







1 



y 



r\o. 2. 
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el centro O en los puntos sucesivos en que ha quedado dividida la recta, y el punto que 
se mueve sobre el diámetro irá describiendo la sinusoide de la figura. 

Si quisiéramos trazar otra curva armónica de período 4 del de la anterior, las partes 
en que se habría de dividir el círculo de referencia, serían de doble número de grados que 
las anteriores, pues es claro que para que la circunferencia quede recorrida en doce horas, 

el punto habrá de recorrer 30° cada hora; y así para curvas de período i, i, etc., del 

de la primera, los respectivos círculos de referencia habrían de dividirse en partes, cuyo 
número de grados fuera triple, cuádruple, etc. ; por donde vemos que los períodos de estas 
curvas están en razón inversa del número de grados de las partes en que se dividen sus 
círculos de referencia. 

Una ordenada cualquiera en la curva, representa la distancia a que se encuentra de 
O el punto P en la hora correspondiente. Designemos por y,, y/, y/\ etc., estas ordenadas. 
Del triángulo rectángulo POM (fig. 1) deducimos como expresión general de las ordena- 
das — 

OP = y -= OM sen a = u sen ( U-]-x) 

donde u es el radio del círculo de referencia, o amplitud de la curva; U la fase inicial, y 
X representa en cada caso el número de partes de la circunferencia que para ese momento 
ha recorrido el punto M. 

Partiendo ahora de estas consideraciones, obtenemos para las distintas ordenadas de 
las curvas armónicas características de una curva periódica cualquiera las expresiones si- 
guientes : 

Primera Segunda Tercera 

ordenada. ordenada. ordenada. 

Primera curva Ui sen Ui Ui sen (Ui + x) Ux sen {Ui+2x) 

Segunda curva U2 sen U2 u-¿ sen {U2+2x) Uz sen (í72+4íc) 

Tercera curva Uz sen C/3 Ws sen (U^+Sx) u» sen (Us+Qx) 



Sumando algebraicamente las ordenadas del mismo nombre, obtenemos para las orde- 
nadas correspondientes de la curva periódica resultante las expresiones siguientes: 

Resultante. 

Primera ordenada = Ui sen Ui + U2 sen U2 + Us sen Us + 

Segunda ordenada = Ut sen (Ui+x) + Wzsen (í/2+2a;) + tin sen {U.i-}-3x) -f 

Tercera ordenada = Ut sen (Ui+2x) + Ui sen (í72+4£c) -f Un sen (Uti-j-ñx) -f 



que, como fácilmente se vé, no son más que expresiones particulares de la fórmula gene- 
ral — 

Ui sen (Ui+x) + U2 sen (U2+2x) + uz sen (í78 + 3x)-f 

en la que el elemento variable x va tomando los valores que le corresponden según lo expli- 
cado. 

Deducida por este procedimiento la fórmula de Bessel, aparece claramente lo que al 
principio dijimos, que es la expresión analítica de una curva periódica en términos de sus 
componentes armónicas ; cada uno de sus términos representa una de estas componentes. 

Las constantes uy U suelen denominarse : la primera coeficiente armónic lineal; la 
segunda coeficiente armónico angular (o ángulo de posición del punto de partida del 
elemento variable), pues ellas son las que caracterizan el movimiento armónico que se 
estudia; la u da la amplitud, y por tanto la importancia del movimiento; pues, como fácil- 
mente se vé en las figuras precedentes, a medida que decrece u, el movimiento del punto 
P sobre el diámetro va siendo más lento; la U (fig. 2) fija la época de las distintas fases; 
un aumento en su valor hace adelantar la época de todas las fases, y al contrario una dis- 
minución la retrasa. 

Es de suma importancia tener presente esta propiedad de la fase inicial ; si en una 
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estación cuya longitud fuera, por ejemplo, 15° W. de Greenwich se quisiera calcular en 
tiempo medio de Greenwich en vez de en tiempo medio local, se adelanta una hora el origen 
de los tiempos, y por consiguiente hay que disminuir el valor de las fases iniciales, 15"* 
en f/i, 30"^ en U.¿, ib"" en U^, etc. Por esta razón debe hacerse constar en las publicaciones 
el tiempo en que se calcula. 

La cantidad M que aparece en la fórmula representa la media aritmética de las 
observaciones que entran en la discusión. Lo que se desarrolla en la serie trigonométrica 
de la fórmula es la diferencia entre cada observación y la media aritmética de todas ellas ; 
por consiguiente, para que nuestras observaciones (o las ordenadas de las curvas de 
observación) queden representadas según esta fórmula, hemos de agregar la media arit- 
mética al desarrollo de estas diferencias. 

Creemos que lo dicho hasta aquí será suficiente para formarse una idea general del 
mecanismo de la fórmula ; en la aplicación que haremos de ella al estudio de la variación 
diurna de la presión atmosférica en Manila procuraremos concretar más estas ideas. 

APLICACIÓN DE LA FÓRMULA AL ESTUDIO ANALÍTICO DE LA VARIACIÓN DIURNA 
DE LA PRESIÓN ATMOSFÉRICA EN MANILA EN EL PERÍODO 1890^1909. 

En la preparación de los datos que habían de entrar directamente en el cálculo hemos 
procedido de la manera siguiente : 

De las publicaciones oficiales de este Observatorio tomamos las medidas horarias de 
cada mes, calculando además por separado la media correspondiente a la hora cero, o 
primera medianoche del día representativo de cada mes; esta media es el valor medio 
de las lecturas de la primera medianoche de cada día del mes. Como la hora cei^o, o 
primera medianoche de un día, es la última medianoche del anterior, en un mes cual- 
quiera (Enero por ejemplo) las lecturas de las 0*^ serán: Para el día 1.^ la lectura de 
las 12 p. m. del último día de Diciembre del año anterior; la del día 2, la de las 12 p. 
m. del día 1." del mismo Enero, etc., de modo que, los valores que entran en el promedio 
de esta hora son los mismos que entran en el de las 12 p. m. a excepción de la lectura 
de las 12 p. m. del último día del mes; que es sustituida por la correspondiente del mes 
anterior. 

Fácilmente se ve que, teniendo ya calculada la media de las 12 p. m. basta, para 
calcular la de las 0^\ restar la lectura de las 12 p. m. del último día del mes de que se 
trata, de la lectura correspondiente en el mes anterior; dividir esta diferencia por el 
número de días del mes en cuestión, y sumar algebraicamente el resultado a la media de 
las 12 p. m. 

Obtenidas estas 25 medias horarias, se han calculado las medias horarias generales 
de los doce meses para todo el período, agrupando los meses del mismo nombre, y pro- 
mediando las medias obtenidas para cada hora de un mismo mes en los distintos años. 
Estas medias horarias generales se han aproximado hasta las milésimas, y son las que 
forman la Tabla I. 

La Tabla II es la misma Tabla I en la que se ha eliminado el elemento aperiódico. 
Vamos a explicar brevemente en qué consiste esta corrección, y cuál es su fin. 

Variación no cíclica, o aperiódica. — En la explicación de la fórmula vimos que las 
curvas armónicas elementales constan en su integridad de 25 ordenadas, correspondientes 
a las distintas horas del día, de O'' a 24\ ambas inclusive. Las ordenadas extremas son 
iguales para una misma curva; por tanto las ordenadas correspondientes a las Cl*' y 24^ 
en la curva resultante han de ser iguales, por ser la suma algebraica de las ordenadas 
correspondientes en las distintas curvas elementales. 

Ahora se entenderá por qué es necesario calcular la media de las O*' ; es preciso ase- 
gurarse de que la curva que se va a analizar cojno periódica, cumple con el carácter esen- 
cial de periodicidad de que sus ordenadas extremas sean iguales. En caso de no serlo, 
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debemos deducir que se ha. introducido en la marcha del fenómeno algún elemento que 
perturba su periodicidad. Esta perturbación es la que suele denominarse variación no 
cíclica o aperiódica. 

Cuál sea de ordinario el origen de esta variación, se entenderá fácilmente con un 
ejemplo bien familiar. 

En las distintas estaciones del año varía la temperatura de modo que va subiendo 
en primavera y verano, y bajando durante el otoño e invierno; así que en Mayo, por 
ejemplo, aunque no puede decirse en general que la temperatura de un día sea más 
elevada que la del precedente, pero si queremos calcular la variación media diurna de 
la temperatura en ese mes para un año cualquiera, tomando para cada hora del día el 
valor medio de aquella hora, derivado de los treinta y un días, notaremos que la tempe- 
ratura de la primera medianoche del día representativo de Mayo es algo más baja que 
la de la última medianoche del mismo día. 

En el caso de observaciones barométricas nos encontramos con un fenómeno aná- 
logo. Es bien sabido que la presión atmosférica está sujeta a una oscilación anual, 
que se sobrepone, como en el caso de la temperatura, a la oscilación diurna, Al estu- 
diar pues los caracteres propios de la variación diurna mensual de la presión atmosférica, 
debemos comenzar por la eliminación del elemento aperiódico. 

Esta eliminación suele hacerse aplicando a cada hora la corrección — 

D{12~h) 
24 

en la que D se obtiene restando del valor de la hora 24 el de la hora cero, y h representa 
la hora en la que se va a aplicar la corrección. No estará demás alguna discusión sobre 
la legitimidad de esta fórmula. Es lógica esta corrección únicamente si se admite la 
hipótesis de que la diferencia entre los valores de las dos horas extremas se ha intro- 
ducido por una variación uniforme durante las veinticuatro horas. 

El probar si es, o no, legítima esta suposición no carece de dificultad; sin embargo, 
cuando la investigación se limita a las veinticuatro horas, se la da por buena, pues parece 
la única hipótesis aceptable. 

Veamos cómo la fórmula supone implícitamente esta hipótesis. Si la variación se 
supone uniforme durante las veinticuatro horas, ha de estar representada por una recta. 
La fórmula de corrección será por lo tanto la ecuación de la recta — 

^=r:6 + aX (1) 

Tomemos como origen de coordenadas rectangulares la hora cero, y llevemos, como 
de ordinario, las horas al eje X, En la fórmula (1) Ir x representará las horas suce- 
sivas; la y, la corrección que hay que aplicar en cada una. Aplicada la corrección a las 
dos horas extremas, sus valores han de venir a coincidir ; luego hay que añadir a la una 
y restar de la otra la semidiferencia de sus valores. La corrección para estas dos horas, 
llamando D su diferencia será por consiguiente — 

para la hora 24 . ± =6+ 24a 

para la hora O „ ^ T ^b 

Eliminando a b entre estas dos ecuaciones obtenemos — 

24 

Llevando a la ecuación (1) los valores obtenidos para a y b, queda — 

D^D Dil2~x) 
^ 2~2r" 24 
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Los signos quedarán determinados por el que tenga D en cada caso particular. 
Vemos pues que, para la deducción de la fórmula de corrección, es necesaria y suficiente 
la hipótesis de que la variación aperiódica se verifique de un modo uniforme durante 
las veinticuatro horas. 

Un ejemplo hará ver qué influencia ejerce esta corrección en las curvas. 

Tomemos el mes de Noviembre, que es en nuestro caso, el que tiene mayor co- 
rrección. La curva llena (fig. 3) representa la variación sin corregir; la punteada 
representa la misma variación corregida. 

Por la forma y posición de estas curvas se vé claramente, que la eliminación del 
elemento aperiódico se reduce a hacer girar la curva un ángulo pequeñísimo, tomando 
las 12 a. m. como centro, sin que en nada varíen sus ondulaciones características. 

Hecha esta corrección, se han calculado para la Tabla II las medias diurnas, divi- 
diendo por 24 la suma de las 24 medias horarias de cada mes, de O*" a 23^; y restando 
luego en cada mes la media diurna, de cada media horaria, se ha formado la Tabla 
III, que es la que contiene los datos de que nos vamos a servir en el cálculo. 

En esta última evaluación hemos aproximado únicamente hasta las centésimas; 
pero, como en los cálculos anteriores llevamos la aproximación hasta las milésimas, 
la inexactitud de los datos de la última tabla no puede pasar de media centésima, y 
la suma de los valores positivos no se diferencia en valor absoluto de la de los nega- 
tivos más que en dos o tres centésimas a lo sumo, resultado que no se hubiera obtenido 
si en los cálculos anteriores se hubiera aproximado únicamente a centésimas. 




r.g3 

En el estudio analítico que con los datos de esta tabla vamos a emprender hemos 
empleado el tiempo medio local. 

Sabemos que lo único que necesitamos para la determinación de cada una de 
las ondas elementales, en las que vamos a resolver la variación barométrica diurna son : 
Los coeficientes armónicos lineales u, y los angulares U — es decir las constantes de la 
Fórmula 11, por ser ellas, como al principio indicamos, las que caracterizan estas ondas. 

Estas constantes no pueden hallarse directamente; hay que determinarlas valién- 
dose de las constantes p y g de la Fórmula I, que son las que se calculan directamente. 

Las ecuaciones que relacionan estas constantes son, como vimos antes : 



P 
Q 



tisen U 

UCOB U 



(1) 



De aquí se obtienen las relaciones: 



sen U .^ rr P. rr P rr Q 

rr tan ¿7^^^ sen U^- ; eos U^=^ 
eos U q %' u 

De la primera de estas tres ecuaciones podemos deducir la constante angular U en 
función de p y g; y de las dos últimas en función de p o <? y u; esta última se obtiene 
por la ecuación : 
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que se deduce fácilmente de las ecuaciones (1) elevándolas al cuadrado, y sumando 
ordenadamente. Toda la dificultad está en determinar las constantes p y q. 

El método que dio Bessel para determinarlas, y que explicaremos al fin de nuestro 
trabajo, es largo y penoso, y sobre todo, requiere una tensión constante en el trabajo 
para evitar los errores en que fácilmente se puede incurrir a causa, principalmente, 
del considerable número de ángulos que hay que reducir y por la frecuencia con que cam- 
bian de signo las líneas trigonométricas al pasar de uno a otro cuadrante. De este 
método dice al P. Dechevrens, S. J. : *'La méthode naturelle n'est que pour les hábiles.'' 
Este Padre ha publicado dos folletos, el uno titulado '*Méthode simplifiée dite des 
facteurs pour le calcul des series de Fourier et de Bessel appliquées a la Météorologie," 
en el que por medio de tablas, se facilita notablemente el cálculo. El otro folleto lleva 
el título "Note complémentaire á la méthode simplifiée du calcul, etc.'' ' donde explica 
el procedimiento que ha seguido para la formación de sus tablas. 

Las tres primeras de estas tablas sirven para la determinación de los coeficientes 
¡) y q, en series de 4, 6, 8, 10, 12, 24 y 36 observaciones. 

Creemos no carecerá de interés, para quienes no estén avezados a esta clase de 
análisis, resolver prácticamente un ejemplo por el método de Bessel, y ver en qué está 
la utilidad de las tablas del P. Dechevrens ; por esto, aún a riesgo de parecer nimios en 
detallar el cálculo, vamos a seguir paso por paso el método de Bessel en la determinación 
de las constantes p y q para Enero. 

Las ecuaciones que en nuestro caso dan los valores de estas constantes para los 
ocho primeros términos de la fórmula (I) son: - 

12Pi=2 a,n eos mx 12(7 1=2 ^m sen mx 

12^2=2 a„, eos m(2x) 12^0=2 a,,, sen m(2x) 

12^3=2 a,n eos m(Sx) 12^3=2 a„, sen m{Sx) 

12p4=2 a^ eos m(4:X) 12q=^ a,„ sen m{4:X) 

El signo 2 indica que la cantidad a que afecta se ha de desarrollar en una serie 
de términos de su misma forma, que se suceden según la serie de los números natu- 
rales O, 1, 2, etc., así que escribiremos: 

2 a,n eos mx^a^^ eos O+^^i eos x-\-a^ eos 2x^a^ eos 3a;+ 

X representa el valor de una de las partes en que está dividido el círculo de referencia 
de la primera curva elemental. Como el período de esta curva es de veinticuatro horas, 
su círculo de referencia está dividido en 24 partes, y por tanto a;=:15°. Los períodos 
de las otras curvas son |, i, J del de la primera; por consiguiente las divisiones de sus 
respectivos círculos de referencia valen 

2x;Zx\ Ax;o 30^; 45°; 60°; 

a„, representa los valores de la Tabla III correspondientes a las distintas horas; así en 
Enero, por ejemplo— 

ao=46, ai=17, a2=— 20, etc. 

Según esto, para la determinación de p^ en Enero, tendremos — 

12p,--46 eos 0+17 eos 15°-20 eos 30°~45 eos 45°. 

Sólo falta hallar en las tablas trigonométricas naturales el valor de estos cosenos, y 
ejecutar las operaciones indicadas. El trabajo, que sería ímprobo si se hubiera de hacer 
la multiplicación en cada término, puede simplificarse mucho del modo siguiente: 

^ Tipografía della Pace di Filipo Cuggiani, Roma, Via della Pace num. 35. 
* Al ñn de este trabajo daremos la prueba matemática de estas ecuaciones. 
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El valor de Icks cosenos en las distintas horas es: 



o*» 




... eos 0** = 


1 


12" 


1 




. eos 15" = 


.966 


13 


2 




. eos 30^ = 


.866 


14 


3 




. eos 45" = 


.707 


15 


4 




. eos 60° = 


.5 


16 


5 




. eos 75° = 


.259 


17 


6 ,. ... 




. eos 90° = 





18 


7 .... 


.... eos 105° = 


-eos 75° = 


- .259 


19 


8 


.... eos 120° = 


— eos 60° = 


- .5 


20 


9 




— 


- .707 

- .866 

- .966 


?A 


10 







?.?. 


n 







23 



= -1 

= ~~ .966 
= - .866 
= ~ .707 
= ~ .5 
= - .259 
= O 
= . ,259 
= .5 
= .707 
= .866 
= .966 



Cada uno de estos números aparece como factor común de varias horas. Tendre- 
mos por tanto : 

12 p, =(0-12)1+ (1^11-13+23). 966+ (2-^10-^ 14+22). 866 

+ (3__9_i5+21) .707+ (4-~8-~16+20) .5+ (5-7-^17+19) .259 

Sustituyendo en vez de las horas que aparecen en los paréntesis, los valores correspon- 
dientes tomados de la Tabla III tenemc^ por fin : 

12 p, =(46-15) 1+ (17-79+63+59) .966+ 

Para la determinación de q^ obtenemos por el mismo procedimiento: 



0^ sen 0^=0 

1 .., , sen 15'*= .259 

2 sen 30°= .5 

3 sen 45**= .707 

4 sen 60°= .866 

5 sen 75^= .966 

6 sen 90" =1 

7 sen 105'* = sen75'*= .966 

8 - = .866 

9 = .707 

10 = .5 

11 = .259 



12V 

13 . 

14 . 

15 . 

16 - 

17 . 

18 . 

19 . 

20 . 

21 . 

22 . 

23 . 



= O 
= " .259 
= - .5 
= ™ .707 
= - .866 
= - .966 
= -1 
= - .966 
= - .866 
= - .707 
= ~ .5 
= - .259 



Teniendo en cuenta los factores comunes de varias horas, escribiremos : 

12 q, = (1+11-^13-23) .259+ (2+10-14-22) .5+ (3+9-15-21) .707 

4. (4+8-^16 -20) .866+ (5+7-17-19) .966+ (6-18) 1 

y sustituyendo en vez de las horas sus valores queda determinada q^. 

Para la determinación de Pg y Q2 procedemos del mismo modo, teniendo en cuenta 
que las divisiones del círculo de referencia valen 30°; diremos pues: 



P2 



o^..... coso** 

1 eos 30° 

2 eos eo'' 

3 eos 90" 

4 eos 120^ = -eos 60'' 

5 ........ eos 150" = -eos 30° 

6 

7 

8 

9 

10 

11 



1 


12" 


.866 


13 


.5 


14 





15 


.6 


16 


.866 


17 


1 


18 


.866 


19 


.5 


20 





21 


.5 


22 


.866 


23 



= 1 

= .866 
= .5 
= O 
= -- .5 
= - .866 
= -1 
= - .866 
= - .5 
= O 
= .5 
= .866 
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12 p, = (0-^6+12-18) 1+ (1^5^-7+11+13-- 17-. 19+23) .866 

+ (2-4-^8+10+14-- 16--20+22) .5 



sen 0°= O 



1 sen 30°= .5 

2 sen 60''= .866 

3 sen 90°= 1 

4 sen 120°=sen 60°= .866 

5 = .5 

6 = O 

7 =- .5 

8 =- .866 

9 =-1 

10 =- .866 

11 =- .5 



12\ 

13 . 

14 . 

15 . 

16 . 

17 . 

18 . 

19 . 

20 . 

21 . 

22 . 

23 . 



= O 
= .5 

= .866 
= 1 
= .866 
= .5 
= O 
= - .5 
= - .866 

= - .866 

= - .5 



12 92 = (l+5~7-ll+13 + 17-19-23).5+ (2+4-8-10+14 

+16-20™22).866+ (3-9 + 15-21)1 

El círculo de referencia de la tercera curva está dividido de 45'' en 45° ; por consi- 
guiente: 



P, 



0^ eos 0°= 1 

1 eos 45°= .707 

2 eos 90°= O 

3 eos 135° = -eos 45° = ~ .707 

4 =-1 

5 =-. .707 

6 = O 

7 = .707 

8 = 1 

9 = .707 

10 = O 

11 =- .707 



12\ 

13 . 

14 . 

15 . 

16 . 

17 . 

18 . 

19 . 

20 . 

21 . 

22 . 

23 . 



= -1 
= - .707 
= O 
= .707 
= 1 
= .707 
= O 
= - .707 

= - .707 
= O 
= .707 



12 P3= (0-4+8.-12+16-20) l+(l-3-5+7+9-ll-13 + 15 + 17-19-21+23).707 



Qa 



O*" sen 0**= O 

1 sen 45°= .707 

2 sen 90°= 1 

3 sen 135°=sen 45°= .707 

4 - = O 

5 :. =- .707 

6 =-1 

7 = - .707 

8 = O 

9 -, = .707 

10 = 1 

11 = .707 



12'' 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 



= O 
= ™ .707 

= - .707 
= O 
= .707 
= 1 
= .707 
= O 
= - .707 
= -~l 
= - .707 



12^3= (1+3-5-7+9+11-13-15+17+19-21-23). 707+ (2-6+10-14+18-22)1 
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Las divisiones del círculo de referencia de la cuarta curva valen 60"^ 
pues : 



tendremos 



P4 



0^' 

1 

2 

3 


.... eos 120° = 


... eos 0^= 1 
. eos 60""= .5 
-eos 60° = ~ .5 
=--1 


12^ 
13 
14 
15 


4 




=- .5 


16 


5 




= .5 


17 


6 




. = 1 


18 


7 




„.. = .5 


19 


8 




=- .5 


20 


9 ..... 






21 


10 




=™ .5 


22 


11 




= .5 


23 



= 1 

= .5 
= ^ .5 
= -1 

= - .5 
= .5 
= 1 
= .5 
= ~ .5 
= -1 
= - .5 
= ,5 



12 p^==(0~-3 + 6-^9+12--15+18-™21)l+ (1-2-^4+5+7-8-10 

+ 11+13-14-16+17+19- 



20-22+23) .5 



Qé 



0»" -, 

1 

2 

3 .. 


sen 0'' = 

sen 60° = 

S3n 120°=sen 60° = 



.866 
.866 



- .866 

- .866 


.866 
.866 


- .866 

- .866 


12'' 
13 
14 
15 


4 





16 


5 ...... 


__ 


17 


6 


— 


18 


7 





19 


8 





20 


9 ........ 





21 


10 





22 


11 


— 


23 



= o 

= .866 
= .866 
= O 
= - .866 
= - .866 
= O 
= .866 
- .866 
= O 
= - .866 
= - .866 



12 q,= (1+2-4-5+7+8-10-11 + 13 + 14—16-17 + 19+20-22-23) .866 

Hemos descendido a tantos pormenores porque pretendemos hacer ver que las tablas 
del P. Dechevrens, para la determinación de estos coeficientes, son de entera confianza, 
por estar fundadas en el método mismo de Bessel, y tienen la ventaja de dar verificadas 
ya las operaciones más expuestas a errores. Las fórmulas que da el P. Dechevrens son : 

'p,= (0-12) .083+ ( 1+23-11-13) .08+ (2+22-10-14) .072+ (3+21-9-15) .059 

+ (4+20-8-16) .042+ (5+19-7-17) .022 
q,= (6-18) .083+ (1+11-13-23) .022+ (2+10-14-22) .042+ (3 + 9-15-21) .059 

+ (4+8-16-20) .072+ (5+7-17-19) .08 

>,=:(0+12-6-18).083+ (1 + 11 + 13+23-5-7-17-19) .072+ (2+10+14 

4-22-4-8-16-20) .042 
q„= (3+15-9-21) .083+ (1 + 5+13 + 17-7-11-19-23) .042+ (2+4 

+ 14+16-8-10-20-22). 072 

P3=(0+8+16-4-12-20).083+ (1 + 7+9+15+17 

+23-3-5-11-13-19-21) .059 
(?,,= (2+10+18-6-14-22) .083+ (1+3+9+11+17 

+19-5-7-13-15-21-23) .059 

fp,= (0+6+12+18-3-9-15-21). 083+ (1+5+7+11+13+17+19+23-2 

-4-8-10-14-16-20-22) .042 
I 9,=-- (1+2+7+8+13+14+19+20-4-5-10-11-16-17-22-23) .072 
Comparándolas con las halladas por el método de Bessel, vemos que son idénticas; 
los coeficientes de los paréntesis son los mismos de Bessel divididos por 12. 



I 



II 



III 



IV 
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Por medio de estas tablas determinamos las constantes que forman la Tabla IV. 

Con estas constantes se pasa a determinar los coeficientes armónicos, por medio de 
las fórmulas que, según antes vimos, los relacionan; damos a continuación las que he- 
mos empleado, y el procedimiento que hemos seguido en el cálculo. 

^ = tg Í7 .\ log p— log q = log tg U 

u = ^|7 •'• I^g ^ = lo? P— log sin U 
En el ejemplo que estamos resolviendo para Enero tendremos, según esto: 

logp, =1,127429 
log q, = 1. 738463 

fog^;Zio'^Qr3"r^"3^^ = log tg Í7, .-. U, = 13° 45^ 37'' 
log p, = 1. 127429 
log sin f7, = 1.376320 



log p, — log sin U, = 1. 751109 = log u, /. ti,= 0.5638 

Hay que tener presente en qué cuadrante se halla U. Puede seguirse esta regla, 
para averiguarlo fácilmente. 

Si p y q son positivas, U pertenece al primer cuadrante ; si p es positiva, y q negati- 
va, al segundo ; si ambas son negativas, al tercero ; finalmente, si p es negativa, y q posi- 
tiva, al cuarto. 

Así, en este mismo mes, al determinar Í7. se hallará directamente en las tablas 
trigonométricas el ángulo 15"^ 23'; pero como Po ^s positiva, y Q2 negativa, hay que re- 
ducirle al segundo cuadrante, lo que da 

U, = 164° 37' 

De este modo se han ido deduciendo para cada mes los valores de u y U, que damos 
en la Tabla V. 

Aunque en rigor basta este cálculo para discutir las curvas elementales, nos ha pare- 
cido calcular por separado para cada mes estas cuatro curvas. 

Los valores de las ordenadas correspondientes a cada hora, juntamente con los de la 
resultante, aparecen en la Tabla VI, en la que hemos añadido la última columna con las 
diferencias entre la variación resultante calculada y la observada. Lo insignificante de 
estas diferencias nos parece una buena comprobación de los cálculos. 

Una vez calculadas estas ondas, juzgamos no carecería de interés añadir las 12 
láminas, en cada una de las cuales aparecen las curvas armónicas elementales de un 
mes con su resultante. Así se podrá apreciar, más fácilmente que por los valores de 
las tablas, la importancia de cada una de las curvas. En ellas se ve claramente por qué 
basta calcular los cuatro primeros términos de la fórmula; la serie es rápidamente con- 
vergente, de modo que la cuarta sinusoide se confunde casi con una recta. 

Aunque, para la verificación de estos últimos cálculos, da también tablas de P. 
Dechevrens, juzgamos que sería, si no más fácil, mucho más breve calcular por el método 
natural. 

Por vía de ejemplo, vamos a calcular las cuatro componentes de Enero. La expre- 
sión general de las ordenadas de la componente de veinticuatro horas es 

u sen (U.-h^) 

por consiguiente, en nuestro caso, tendremos como expresión general 

0.564 sen (13° 46' +x) 

Para determinar en las distintas horas el águlo expresado por el paréntesis, basta re- 
cordar que, en esta primera onda, el punto que recorre el círculo de referencia, ha 

127182 5 
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46'= .238 


7\ 


. senll8M6' = 


sen 61° 14' = 


.877 


46'= .481 


8 


.sen 133° 46' = 


sen 46° 14' = 


.722 


46'= .692 


9 . 


.sen 148° 46' = 


sen 31° 14' = 


.519 


46'= .855 


10 , 


.sen 163° 46' = 


sen 16° 14' = 


.280 


46'= .960 


11 


. senl78M6' = 


sen 1° 14' = 


.022 


46' = 1 


12 . 


. sen 193' 46'= - 


-sen 13° 46' = - 


-.238 


14'= .971 











descrito un arco de IS"" al fin de la primera hora, 30° al fin de la segunda; 45"^ al fin 
de la tercera, etc., de modo que se obtendrá para las distintas horas : 

0\.. sen 13° 

1 ..-- , sen 28" 

2 sen 43° 

3 sen 58" 

4 .-...- sen 73^ 

5 sen 88" 

6 sen 103" 46' = sen 16° 

De las 12^ en adelante se van repitiendo los valores en el mismo orden, pero con 
signo contrario. Basta ahora multiplicar cada una de las funciones naturales, que apa- 
recen en último miembro, por el coeficiente 0.564, para obtener los valores que repre- 
sentan las distintas ordenadas de esta onda. 
Segunda onda. — La expresión general es: 

u, sen {U.-\--2x) 
que en nuestro caso se reduce a 

0.978 sen (164^^ m'-^2x) 
El elemento variable 2x aumenta en 30"^ cada hora; por tanto 



O". .... sen 164" 37'= sen 15' 23'= .265 

1 ...... sen 194 37'= -sen 14^^ 37'= --^.252 

2 ... sen 224 37'= -sen 44 37'.= -.702 

3 sen 254 37'= --sen 74 37'=™-.964 



4" .... sen 284^^ 37'=--sen 75° 23'= -.968 

5 ..... sen 314" 37'= -sen 45° 23'= -.712 

6 ..... sen 344^ 37' = -sen 15° 



23' 



.265 



Desde las 6'' comienzan a repetirse los valores, pero con signo contrario ; de modo que, 
para obtener las ordenadas de esta curva, basta multiplicar estos valores por 0.978, e 
ir repitiendo esta serie de productos hasta la hora 23, teniendo cuidado de cambiar de 
signo en cada repetición. 

Las dos ondas restantes se obtienen por análogo procedimiento teniendo en cuenta 
el valor del elemento variable en cada una. 

Tercera onda, 

O'' .......^ sen ir^ 1' = 

1 .,..,..-....^...,. sen 56° 1' = 

2 ...... sen 101" 1'= sen 78' 59' = 

3 ..... sen 146 1'= sen 33' 59' = 

4 ... sen 191' l'=^~-sen 11" 1 = 









Cuarta onda. 




191 


0" 


- sen 130' 


17'= -sen 49° 


43'= -.763 


829 


1 . 


. sen 190" 


17'= -sen 10° 


17'= -.179 


982 


2 . 


.. sen 250° 


17'= -sen 70° 


17'= -.941 


569 


3 . 


-. sen 310° 


17'= -sen 49° 


43'= -.763 


191 











Este es el modo natural de proceder en el cálculo; en realidad es mucho más sen- 
cillo, si nos fijamos en que los ángulos van siendo complementarios de los precedentes, 
desde las 6" en la primera onda ; desde las 3" en la segunda, y desde las 2" en la tercera! 
De modo que basta determinar los ángulos hasta las 5" inclusive en la primera onda • 
hasta las 2" en la segunda, y hasta la 1- en la tercera, y sin cuidarse de los restantes! 
tomar a la vez, en las tablas de las funciones trigonométricas naturales, el valor 
del seno de cada uno de estos ángulos, y el de su complementario, que se encuentra 
al extremo de la misma línea. Por ejemplo, para la hora cero, en la primera onda de 
Enero, hallamos el valor del seno de 13' 46', que es 0.238, y al otro extremo de la 
misma linea, encontramos el valor 0.971 correspondiente al seno de su complementario 
76^ 14', que escribiremos al frente de la hora 6. 

Por medio de las tablas del P. Dechevrens, para hallar la ordenada correspondiente 
a las O" de Enero, por ejemplo, hay que aplicar la siguiente fórmula: 

(0-12).083+(l-13).08+(2--14).072+(3-15).059+(4-16).042+(5-17).022 

+ (7-19) X -.022+ (8-20) X -.042+ (9-21) 
X-.059+(10~22)x-.072+ (11-23) X-.08 
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sustituir en vez de las horas que aparecen en los paréntesis, sus valores, y ejecutar las 
operaciones indicadas. En nuestro caso, hechas ya las restas indicadas en los paréntesis, 
tendríamos 

.31X.083+.8X.08+1.05X.072+1.1X.059+1.04X042 

+.89X.022-.84X.022--.87X.042-.79X.059™.58X.072~-.2X.08=0.1341 

Por el método natural, una vez hallados los valores de las líneas trigonométricas, 
obtenemos el mismo resultado con una sola multiplicación: 

0.238 X0.564=:0.1342 

Claro está que el método del P. Dechevrens es menos expuesto a errores, por 
reducirse a sencillas operaciones aritméticas; pero creemos que basta alguna práctica 
en el cálculo trigonométrico, para ganar tiempo por el método natural. 

Épocas de las fases críticas (máximos y mínimos). — La determinación de las 
épocas de las fases críticas no ofrece dificultad. Teniendo presente que una ordenada 
cualquiera tiene por expresión 

u sen {U~\-x) 

la condición necesaria y suficiente, para que este valor sea máximo, es que el seno que 
aparece en la expresión general, valga la unidad, lo que deja finalmente como condición 
para el máximo. 

U+x=90'' 6 450^ 

según que U sea menor o mayor que 90"^. Esto supuesto, recordemos que el elemento 

variable x tiene por valores para las horas sucesivas O, 1, 2, 3, h, en la primera 

curva; O'', 15°, 2x15°, 3X15°, /¿Xl5°, del mismo modo en la segunda, 

tercera y cuarta tiene respectivamente por valor para una hora cualquiera /i, las ex- 
presiones /iX30°, /iX45°, /ix60°. Según esto, vamos a determinar las épocas del pri- 
mer máximo para las cuatro curvas de Enero. 
Para la primera tenemos como condición: 



13" 46'+ftX18"=90».-.fc?5:^^°^í«: = ^«"^y-'-S» 6- 



para la segunda: 

para la tercera: 
cuarta : 



164° S7'+hxS0^==^450\\h=^^^^^ =9^ 31" 



11° l'+;ix45°=:90°.'./i=P 46" 



130° 17'+/tX60°=450°.-./i=5»^ 20™ 
Una vez determinada la época del primer máximo, quedan determinadas las restantes, 
por distar entre sí doce horas en la primera curva, seis en la segunda, cuatro en la 
tercera, y tres en la cuarta. Las épocas del primer máximo, para las cuatro curvas 
en los distintos meses, aparecen en la siguiente tabla: 

Época del primer máximo. 



Meses. 



Enero .. 
Febrero 
Marzo .. 
Abril ... 
Mayo--- 
Junio._. 



Primera 


Sesrunda 


Tercera 


Cuarta 


onda. 


onda. 


onda. 


onda. 


H, m. 


H. m. 


H. m. 


H. m. 


5 5 


9 81 


1 46 


5 20 


5 47 


9 41 


1 44 


5 24 


6 23 


9 45 


1 10 


5 87 


5 3 


9 44 


49 


6 38 


4 10 


9 44 


7 15 


6 59 


3 11 


9 49 


6 20 






Meses. 



Julio 

Agrosto 

Septiembre 
Octubre — . 
Noviembre 
Diciembre . 



Primera 
onda. 



2 57 
8 9 

3 25 

4 16 
4 20 
4 26 



SeflTunda 
onda. 



ff, m. 
9 52 



Tercera 


Cuarta 


onda. 


onda. 


H. m. 


H. m. 


6 27 


6 32 


6 52 


4 49 


1 23 


4 38 


1 11 


4 42 


1 19 


4 42 


1 40 


4 46 
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DISCUSIÓN. 

En lo que precede, hemos resuelto la oscilación diurna de la presión en Manila en 
una serie de oscilaciones armónicas simples, las dos primeras de las cuales, de períodos de 
veinticuatro y doce horas, son mas importantes sin comparación que las de períodos de 
ocho y seis horas; de modo que los principales caracteres de la oscilación diurna de la 
presión, están representados casi exclusivamente por las dos primeras. 

Pero este análisis de las curvas no es de suyo otra cosa que un artificio de cálculo; 
toda curva periódica puede someterse a él. Por lo tanto, las curvas armónicas componentes 
pudieran tener una significación puramente matemática, sin corresponder a nada real. 

Helmholtz al tratar de esta clase de análisis dice : '*No es más que una ficción mate- 
mática, admirable por lo que facilita el cálculo, pero no guarda relación necesaria con algo 
real.** Sin embargo, hablando en general, no puede negarse la utilidad práctica de este 
análisis. 

Un fenómeno periódico no puede resultar sino de causas periódicas, y de períodos 
submúltiplos del período del fenómeno ; de lo contrario, la periodicidad de este quedaría 
perturbada. Por medio de este análisis, se descompone la curva en las únicas elementales, 
que son posibles, de períodos submúltiplos del suyo ; por tanto, es sumamente probable que 
las causas parciales que influyen en el fenómeno aparezcan, o totalmente, o en sus caracteres 
principales, en la manera de variar de las curvas elementales. Si deseáramos averiguar si 
dos fenómenos están, o no, relacionados entre sí, el estudio de las armónicas características 
de ambos, y especialmente su modo de variar durante el año, es el único procedimiento 
por donde podemos llegar a cerciorarnos de si existe, o no, esa relación; el que en las cur- 
vas que los representan no aparezcan rastros de relación mutua, no nos dice mucho ; pues 
la influencia que el uno pudiera tener en el otro, puede fácilmente estar contrarrestada por 
otras causas que influyan a la vez. 

Aunque de la comparación de dos fenómenos por medio de este análisis no sacáramos 
más que la persuasión de que no están relacionados entre sí, bien se ve que no se podría 
dar por perdido nuestro trabajo. 

Viniendo al caso particular del barómetro, vamos a traducir lo que el Dr. Hann siente 
de este estudio:' 

Por mi parte, estaba persuadido de que todos los esfuerzos que se hicieran para explicar la 
oscilación diurna del barómetro por medio de las variaciones diarias de los elementos meteorológicos 
en una región cualquiera, como lo intentaron Kreil, Blanford, Renou y otros, no podía llevar a con- 
clusión alguna; he publicado una serie de escritos, en los que doy una descripción precisa del fenómeno 
según se manifiesta en toda la tierra, lo mismo al nivel del mar que a todas las alturas para las 
que existen observaciones, y me he esforzado en presentar los resultados en una forma tal, que pudiera 
servir de base a una teoría físico-matemática. Con este objeto, he representado todos los resultados 
de las observaciones por funciones periódicas, y he calculado las amplitudes y épocas de las fases de 
las ondas individuales, con sus períodos propios, las cuales, al combinarse, producen el resultado com- 
plejo de la curva barométrica diaria, tal como se nos presenta a la observación directa * * * 
Este es el único medio con que nos ha sido posible hacer progresos notables en la comprensión y expli- 
cacion del fenómeno. 

Lo que principalmente nos proponemos en este trabajo es comparar los resultados 
obtenidos para Manila con los que se han obtenido para otras regiones; seguiremos pues 
en la discusión una forma tal, que se preste a esta comparación. 

Primrm midch^üno de los hechos más fundamentales que se establecen en meteoro- 
logía es, quv Uxs capas inferiores de la atmósfera reciben su calor de la superficie terrestre 
por conducción, y lo transmiten luego por convección a las inmediatas. Este movimiento 
'l^/''^^^'''''^^^'^ ^^ principalmente ascendente; el aire se expansiona lateralmente en la alta 

;^' J^tniher contributions to the foundation of a theory of the dail^ b¡^rometHc osdl^^^ 
nal of the Meteorological Society. London, 1B99, p. 40. 
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atmósfera, y determina en consecuencia un descenso en el barómetro. Lo contrario sucede 
con el enfriamiento de la superficie terrestre ; el aire en contacto con ella se enfría, y con- 
trae, dando lugar a corrientes descendentes y ascenso en el barómetro. Por lo tanto, 
en virtud de esta variación de temperatura, han de producirse en la presión atmosférica 
un máximo y un mínimo, hacia las horas del mínimo y del máximo de temperatura. 

El cambio de calor entre la superficie terrestre y las capas inferiores de la atmósfera 
es de carácter periódico, siendo su principal período de veinticuatro horas. Decimos su 
principal período, porque esta variación no es simplemente armónica; de modo que si se 
analizara por separado, no quedaría representada por un sencillo movimiento armónico, 
sino que, además de una componente fundamental de período de veinticuatro horas, resul- 
tarían otras menos importantes de períodos de doce, ocho, seis, horas. 

La primera onda en el análisis de la curva barométrica es de período de veinticuatro 
horas, y presenta invariablemente el máximo y el mínimo hacia las horas de mínima y 
máxima temperatura. Este carácter general de la curva hace ya sospechar, que representa 
el carácter principal de la variación de presión, causada por el cambio de calor entre la 
superficie terrestre y las capas inferiores de la atmósfera. 

Tanto la amplitud como las épocas de las fases críticas de esta onda presentan en 
Manila caracteres especiales, que no podemos menos de hacer notar. La amplitud es 
menor, y la época del máximo ocurre más temprano que en la mayor parte de las estaciones, 
en las que se ha verificado este análisis. Hann hizo notar estas mismas propiedades para 
Mauricio (isla) Singapore y Batavia, en su publicación "Untersuchungen über den 
Táglichen Gang des Barometers.'* 

Mr, Eliot incluye estas tres estaciones en la tabla siguiente, que copiamos del t.X de 
^'Indian Meteorological Memoirs," p. 73, con varias de las estudiadas por él en la India, 
para comparar con ellas a Trivandrum, única estación de la India que presenta estos 
caracteres. A estas cuatro estaciones puede añadirse Manila como aparecerá por los va- 
lores, de la tabla. 



Estación. 



Latitud 
norte. 



Trivandrum.- - 8 

Madras --. 13 

Bombay 

Rangoon 

Chittagongr 22 

Calcuta (All pore) 22 

Kurrachee-.- 24 

Aden.- 12 



31 

4 

18 54 

16 46 

21 

32 

47 

45 



Valor 
medio 
anual de 
la ampli- 
tud u\. 



0.413 
.595 
.472 
.685 
.582 
.667 
,497 
,779 



Época del 
máximo. 



H, m. 



6 56 

7 5 
7 26 



Valor 
medio 
anual de 
la fase 
inicial 
Uu 



Estación. 



21 2 I Allahabad . 

359 38 i I Maí?pur ... 

330 53 I j Lahore .... 

355 3 i Mauricio. 



Latitud 
norte. 



25 26 
21 9 

31 34 

20 6S 



347 45 !i Sinprapore ,„. 1 17 

346 5 i Batavia _.„.,._. 6 US 

343 49 i Manila _ 14 35 

338 24 II I 

ij ! 



Valor 
medio 
anual de 
la ampli- 
tud MI. 



0.774 
.767 
.667 
.310 
.531 
.626 
.544 





Valor 




medio 


Época del 


anual de 


máximo. 


la fase 




inicia! 




Uu 


//. m. 


o ' 


7 38 


mB 33 


7 4 


343 52 


8 32 


322 2 


2 54 


46 35 


4 18 


25 35 


i 4 19 


1 25 17 



4 16 
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Mr. Eliot termina la discusión de esta onda con las siguientes palabras : 

La primera componente de la oscilación diurna en Trivandrum presenta, por tanto, propiedades 
características que difieren grandemente de las otras estaciones de la India, pero que se asemejan de 
una manera notable a las de Batavia y Singapore, y también (aunque no tanto) a Mauricio. Los 
principales caracteres son: su pequeña amplitud, y la temprana época del máximo. Esto probable- 
mente se debe en gran parte, sino del todo, a su posición marítima y su baja latitud. 

Esta semejanza, que hace notar Mr. Eliot entre Trivandrum, Batavia y Singapore, 
aparece mucho más marcada en Manila, tanto en la amplitud, como en la fase inicial y 
época del máximo ; es pues Manila una confirmación más de estas propiedades que se sos- 
pechan para estaciones costeras de baja latitud. 

Vai'iación amial de u^. — La variación anual de u^ y la variación anual de la absorción 
del calor en las capas inferiores de la atmósfera guardan entre sí relación de causa a 
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efecto, como lo han hecho notar Hann C^Untersuchungen, etc/^ y Angot ("Annales du 
Bureau central Meteorologique de France/' 1887, p. B. 297), y puede verse confirmado 
en los trabajos de Mr. Eliot ('Indian Meteorological Memoirs,'' vol. XII, p, 347). 

Breves consideraciones harán ver que se verifica lo mismo en Manila. Como puede 
verse en la Climatología publicada por el P. Coronas en el t. II de '*E1 Archipiélago 
Filipino," los meses de Febrero, Marzo y Abril son los más secos, y de nubosidad 
mínima; Julio, Agosto y Septiembre por el contrario, son los más lluviosos, y de nubo- 
sidad máxima. Las circunstancias pues, son sumamente favorables en los tres primeros 
meses para la absorción del calor en las capas inferiores de la atmósfera, y desfavorables 
en los tres últimos, por interceptar la gran cantidad de nubes buena parte de la energía 
calorífica, y por haberse de convertir parte de la que llega a la superficie terrestre en 
calor latente de evaporación. 

Una ojeada a los valores de u,, hace ver que los tres meses en que alcanza sus valores 
máximos son Febrero, Marzo y Abril, en los que la absorción de la energía calorífica en 
las capas inferiores de la atmósfera ha de ser máxima; y llega a sus valores mínimos 
en Julio, Agosto y Septiembre, época en que la absorción es mínima. 

Mr, Eliot hace notar en las publicaciones antes citadas, que en el tránsito de la 
estación lluviosa a la estación seca suele verificarse un máximo relativo en los valores 
de w,. Este máximo relativo aparece efectivamente en Manila en Octubre, mes que puede 
considerarse como divisorio entre las dos estaciones. 

Segunda onda. — Dos son los caracteres que principalmente llaman la atención en 
esta onda; primero, ser notablemente más importante que las otras tres en todos los 
meses, y en todas las regiones para las que se ha verificado este análisis; segundo, la 
constancia en la época de sus fases criticas. 

Véanse en el siguiente cuadro las épocas del primer máximo halladas por Mr. Eliot 
para varias estaciones de la India: 



latitud 
NorUí. 



Hstat'i*'»». 



Knero. í *^*^' Marxo. Abril. ' Mayo. Junio. 



;/, m. lí. m. //. m. If. m. II. m. //. m. 



8 ;í1 Trivandrum 

10 m ; Trichinop^íly 

\:\ 4 '' Madras 

IH f»l Hondmy 

24 47 Kurrach<?t? . . 

.i4 t(» í.«?h 

Valor rniHlio para la India Tropica! 

14 :í5 Manila . 



lio 

50 



9 3« 
9 59 



9 45 

9 4:í 



9 37 9 



9 m 



9 41 9 58 



9 4H 
♦.* 42 



9 49 

9 54 

9 57 

9 58 



52 
48 
54 



9 58 
9 52 



42 
52 
9 53 
10 2 
lü 10 
9 50 



Julio, :Agosto. 



//. m. H. m. 



9 47 9 47 9 54 



9 ;n 9 41 



10 22 
9 55 



9 45 9 44 9 44 9 49 9 52 9 51 



Sep- 


tiem- 


bre. 


H. 


m. 


9 


29 


9 


48 


9 


39 


9 


45 


9 


49 


9 


47 


9 


41 


9 


41 



Octu- 


No- 


bre. 


vtem- 
bre. 


//. m. 


H. m. 


9 16 


9 14 


9 38 


9 39 


9 28 


9 27 


9 25 


9 23 


9 28 


9 29 


9 49 


9 49 


i 9 29 


9 28 


! 9 30 

1 


9 25 



Diciem 


bre. 


H 


. m. 


9 


22 


9 


47 


9 


33 


9 


36 


9 


35 


9 


49 


9 


36 


9 


26 



Ksta onda t s, por su doble oscilación diurna, la que produce las mareas caracterís- 
ticas del barómetro. Es por tanto de sumo interés averiguar su origen. 

Líis discusiones en (|ue entran Hann, Angot y Lord Kelvin para averiguar el origen 
de esta onda, son demasiado largas para que podamos reproducirlas aquí en todos sus 
detalles. \'anu)s a intentar resumirlas de la manera más breve y clara que nos sea 
posible. 

¥A fenómeno, (jue esté representado por esta onda, ha de tener estas dos propiedades: 
Primera, influir más que ninguno otro en el barómetro; segunda, ser independiente de 
circunstancias locales. 

Desde luego, no se ve qué otro origen pueda tener la causa de este fenómeno, que 
la influencia del sol sobre la atmósfera terrestre. Pero no puede deberse, dice Lord 
Kelvin, a la atracción del sol que origina las mareas, pues en ese caso encontraríamos 
una influencia mucho mayor de la misma especie, debida a la luna, siendo así que, en 
realidad la marea barométrica lunar es nula, o alonas perceptible. Parece pues cierto 
que ha de atribuirse a la variación de la temperatura. Ahora bien ; esta variación en las 
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capas inferiores de la atmósfera está sujeta a todas las variaciones posibles a que dan 
lugar las circunstancias locales de cada estación, mientras que la acción periódica de 
los rayos solares que, día tras día se verifica en las capas superiores, ha de producir 
necesariamente movimientos periódicos de gran regularidad. 

Esto explica fácilmente la regularidad de esta onda. Lo que no aparece tan claro 
es el porqué de su doble oscilación diaria y su notable amplitud. La oscilación de la 
temperatura en las capas superiores de la atmósfera ha de ser, en su forma general, 
análoga a la oscilación de la temperatura en las inferiores, pues las causas que la pro- 
ducen son las mismas ; la insolación durante el día, y la radiación durante la noche. La 
primera hace que el máximo se verifique poco después que el sol adquiere su altura 
máxima ; la segunda produce el mínimo poco antes de la salida del sol. ¿ Cómo pues esta 
oscilación térmica, de un solo máximo y un solo mínimo, ha de poder dar origen a la 
doble oscilación barométrica que estudiamos? 

Para solucionar la dificultad, Hann acude al único medio, que como indicamos antes, 
existe para averiguar la relación que guardan entre sí dos fenómenos. Resuelve la 
variación diurna de la temperatura en sus componentes armónicas, y halla que, de las 
dos ondas principales que completan el fenómeno, la primera (de período de veinticuatro 
horas) aparece constantemente con una amplitud mayor que la segunda (de período de 
doce horas) inversamente a lo que sucede con la presión. He aquí algunos ejemplos: 

Latitud 6° O', Pacífico ecuatorial 0°.87 OMl 

Latitud 16° O', Pacífico tropical r.39 O .31 

Latitud 33° O', Pacífico subtropical 0.97 O .25 

Latitud 48° 12', Zona templada (Viena) 2°.86 0°.51 

A pesar de esta discordancia en las componentes armónicas de los dos fenómenos, 
Hann no duda en asignar como causa de la doble oscilación diaria del barómetro esta 
doble oscilación diaria de la temperatura. 

Al aparecer esta afirmación del eminente meteorólogo, varios científicos opusieron a 
su hipótesis estas dos dificultades : Primera, la variación térmica que se analiza es la que 
se observa en las capas inferiores; la que se trata de estudiar es la de las superiores, 
y puede dudarse de si la marcha de esta segunda es precisamente del mismo carácter 
que la primera. Segunda, la doble oscilación diaria de la temperatura es únicamente un 
resultado del cálculo, una mera ficción matemática, que no puede servir de base a la ex- 
plicación de un fenómeno real. En realidad no existe variación diaria de temperatura 
con dos máximos y dos mínimos; y si, a pesar de esto, obtenemos una doble onda para 
la variación de la temperatura, por empeñarnos en representarla por una serie de senos, los 
resultados que obtenemos son matemáticamente obligados, y no pueden servir para ex- 
plicar un fenómeno de observación ; la causa de este debe buscarse en algo real y natural. 

En cuanto a la primera dificultad, queda resuelta con lo que antes indicamos, acerca 
de la forma de la variación de la temperatura en las capas superiores; atendiendo a las 
causas que la producen, ha de ser de la misma forma general que la que observamos en 
las inferiores. 

La solución de la segunda pide alguna mayor explicación. Es sabido que, sí una 
masa limitada de fluido, líquido o gaseoso, entra en oscilaciones simples o pendulares, 
las amplitudes de estas oscilaciones obedecen a las condiciones propias del fluido (sus 
dimensiones, temperatura, etc.), de modo que, si la masa del fluido recibe un impulso 
único, se moverá siempre con aquella clase de ondulaciones, que le son posibles en virtud 
de sus condiciones propias. Estas son las que suelen denominarse oscilaciones libres. 

Si el impulso se verifica periódicamente, se fuerza a la masa a entrar en oscilaciones 
del mismo período que el del impulso, aunque estas nuevas oscilaciones no coincidan con 
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la forma de oscilación a que pertenecen las ondas libres. Si este impulso periódico no 
fuera de carácter simplemente armónico— es decir, si no pudiera quedar representado 
por una simple sinusoide, podría resolverse, por el teorema de Fourier, en sus compo- 
nentes armónicas de períodos 1, |, h etc., del suyo, de modo que el efecto debido a este 
impulso sería el mismo que producirían oblando a la vez varios impulsos simplemente 
armónicos, cada uno de los cuales estuviera representado por una de las componentes 
armónicas. Nótese bien que está completamente fuera de propósito la cuestión de si 
se dan, o no, en la práctica esos impulsos componentes; basta que el impulso real pueda 
considerarse como resultante de ellos. 

Si el período de alguno de estos impulsos componentes coincide con el de las oscila- 
ciones libres de la masa del fluido, la oscilación forzada resultante adquiere una ampli- 
tud desproporcionalmente grande; de modo que, aunque la amplitud de este impulso sea 
relativamente pequeña comparada con la de los otros impulsos componentes, el efecto 
por él producido es mucho mayor; por donde se ve claramente que, si se analizara por 
medio de la serie de senos la oscilación compleja resultante, la curva armónica corres- 
pondiente a las oscilaciones libres sería de gran amplitud, mientras que la correspon- 
diente en el análisis del impulso sería de amplitud pequeña, y no por esto dejaría de ser 
la segunda, causa de la primera. 

Viniendo ahora a nuestro caso : La curva que representa la variación diurna de la 
temperatura en las capas superiores de la atmósfera, ha de ser necesariamente asimé- 
trica, por serlo sus dos causas, la insolación y radiación, que dependen de la duración 
del día y de la noche. Si pues, al analizarla, vemos que queda reproducida casi exclu- 
sivamente por la combinación de las dos componentes armónicas de períodos de veinti- 
cuatro y doce horas, podemos considerar la variación real que observamos, como la re- 
sultante de dos variaciones armónicas de la temperatura; una de período de veinti- 
cuatro horas, y otra de doce; no importa que no existan de hecho dos causas, que den 
origen a estas dos variaciones; basta, para nuestro propósito, que el fenómeno se pro- 
duzca del mismo modo que si en realidad existieran. 

Ahora bien; el Dr. Margules, siguiendo las indicaciones de Lord Kelvin, emprendió 
el difícil y laborioso trabajo de calcular las oscilaciones que en la atmósfera terrestre 
ha de causar su calefacción y enfriamiento periódico.' En la primera de sus publica- 
ciones probó que el período de oscilación Ubre de nuestra atmósfera es de doce horas 
próximamente. 

Estos trabajos y conclusiones del Dr. Margules han sido acogidos con entusiasmo, 
y han dado ocasión a tan interesantes publicaciones como la de Paul Jaerísch "Zur 
Theorie der hydrodynamischen Cileichungen in sphárischen Koordinaten'' (Meteorologis- 
che Zeitschrift. 1907, pp. 481-498) ; la de W. Trabert **D¡e Theorie der táglichen Luft- 
druckschwankung von Margules, und die tagliche Oszillation der Luftmassen*' (Met. Zeit., 
1903, pp. 481--501) ; la de W. J, Humphreys *'0n the diurnal variations of atmospheric 
pressure*' (Bulletin of the Mount Weather Observatory, 1912, vol 5, part 2, pp. 132-156) 
y la más completa quizá de todas de Horace Lamb **Atmaspheric Oscillations" (Proc. Roy. 
Soc, 191K vol. 84, pp. 551^572).= 

' **tUH»r áh* Schwin^ii^on periodisch erwarniter Luft." SitzunjBfsb. Wien., Akad. Marzo, 1890. 
**Luftb<»weirunirc*n iti einer rotirenden Spharoidschale." Ibid. Abril, 1890; Enero y Diciembre, 1893. 
Puede verse una traducción al iiiKlés en **The Mechanics of the Earth's Atmosphere'' by Cleveland 
Abbe. Smithsonian misceUaneous collections — 843— pág". 296. 

' Estando ya en hi imprenta nuestro trabajo se ha recibido en este Observatorio una publicación 
de Prof, Frank H. Hi|yreIow **I.a termodinámica de la circulación, y la radiación de la atmósfera 
terrestre" (Boletín de la Oficina Meteorológica Arprentina No. 3). En la pá^na 15 propone dos di- 
ficultades contra la teoría: (1) La?i ondas semidiurnas de presión desaparecen dentro de los 2,000 
metros de la superñcie: C2) el término de calor (Q,~Q,*) no estaba incluido en ía ecuación g^eneral 
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Resumiendo pues brevemente, podemos explicar la doble oscilación y notable am- 
plitud de la segunda onda del modo siguiente: Al recibir nuestra atmósfera un impulso 
térmico (si es lícita esta denominación) se produce en ella el mismo efecto que si re- 
cibiera simultáneamente dos, que pudieran estar representados por sinusoides de perío- 
dos de veinticuatro y doce horas respectivamente. La amplitud o intensidad del segun- 
do es menor que la del primero, pero la oscilación que en la atmósfera causa, ha de ser 
de gran amplitud, por coincidir su período con el de las oscilaciones libres, propias de 
nuestra atmósfera; mientras que las oscilaciones producidas por el primero quedan no- 
tablemente reducidas, por no coincidir su período con el de las oscilaciones libres. 

La hipótesis de que la segunda onda se deba a la absorción y radiación de calor en 
la alta atmósfera, está confirmada por la comparación de la variación anual de estas, y 
la de la amplitud 2^.. Bastarán a nuestro propósito breves indicaciones.' 

Uno de los factores más importantes en la influencia que la absorción y radiación 
pueden tener en los cambios de presión, es la presencia o ausencia de nubes en la atmós- 
fera. La cantidad de calor, que absorbe una atmósfera seca y despejada, se emplea en 
dilatarla, dando lugar a movimientos ascendentes y laterales, que se traslucen en el 
barómetro por cambios de presión. Al contrario; la energía calorífica que llega a una 
atmósfera cargada de nubes, se convierte en gran parte en calor latente de evaporación, lo 
que no puede afectar en nada a la presión. Por tanto, en la estación seca, los cambios de 
presión debidos a estas causas serán mayores que en la estación lluviosa, por la notable 
diferencia de nubosidad entre estas dos estaciones. 

La Tabla V hacer ver claramente esta diferencia en los valores, de u^, que son cons- 
tantemente mayores en los meses de la estación seca que en los de la lluviosa, y en los 
que resalta sobre todo la separación perfectamente bien definida de las dos estaciones 
en el salto de Mayo a Junio que en estos valores se verifica. 

Vamos a terminar la discusión de esta onda comparando sus principales caracteres 
en Manila con los que presenta en las mismas estaciones que comparábamos en la primera 
onda. 

El cuadro que sigue está copiado del Tomo X de 'Indian Meteorological Memoirs,'' 

de movimiento empleada por Margules, Jaerisch y Gold. No es la primera vez que Mr. Bigelow 
impugna esta teoría; véase por ejemplo su publicación ^'Studies on the diurnal periods in the lower 
strata of the Atmosphere^' (Monthly Weather Review, 1905, p. 93) y el capítulo 9, página 458 de su 
''Report on the internacional cloud observations," 1898-99. No parece se pueda suponer que científicos 
como Lord Kelvin y Hann ignoraran que la oscilación semidiurna va disminuyendo con la altura 
hasta hacerse nula o casi imperceptible; y a pesar de esto no dudaron en establecer su teoría. 
Deberían pues suponer que la oscilación semidiurna reaparece en la alta atmósfera. Algo de esto 
parece indicar Mr. Ernest Gold en varios pasajes de su interesante memoria, 'The International 
Kite and Balloon ascents" que obtuvo el primer premio ofrecido en 1908 por la Sociedad Meteorológica 
Alemana. Para abreviar nos limitaremos a traducir dos líneas de la página 120 "De aquí que cam- 
bios de presión en la superficie terrestre van unidos a cambios en dirección opuesta en la estratoesfera, 
que dan lugar a corrientes entre la estratoesfera y la tropoesfera." 

La segunda dificultad pudiera inducir a creer que a los tres eminentes matemáticos que nombra 
se les había pasado por alto la influencia que el término de calor libre pudiera tener en la ecuación 
general de movimiento, siendo así que el Dr. Margules después de introducir en la ecuación diferencial 
de movimiento el término de calor libre termina así: "Para nuestro propósito será más conveniente 
considerar las variaciones de presión como consecuencia de las variaciones de temperatura, no como 
consecuencia del flujo variable de calor. Por tanto volvemos a la ecuación (4a)'* — donde ciertamente 
ya no entra el término de que habla Mr. Bigelow. 

' Véase una larga y detallada discusión sobre el particular en el Tomo XII de "Indian Meteorologi- 
cal Memoirs, Calcutta," páginas 328-347. 

Pueden consultarse también Langley: "Researches on solar heat and its absorption by the earth's 
atmosphere** publicado en el No. 15 de "Professional Papers of the Signal Service,'* y en el Vol. 
I de "Annals of the Astrophysical Observatory of the Smithsonian Institution," Cleveland Abbe 
"Preparatory Studies for Deductive Methods in Storm and Weather Predictions.*' 
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página 74, en el que hemos reducido a milímetros los valores que en el original están 
dados en pulgadas para facilitar la comparación. 



Ettación. 

Trivandrum -..« .-„ . 

Madras . . , . . , 


1 

p 
1 

mm. 

1. 100 

1.115 

,977 

h(m 

.946 

i Am 

.797 

.mm 

.918 
1. 010 

.628 

.733 

.9% 

.969 

.919 


B 

& 

o ' 

164 14 

158 5 
166 51 

153 3 

148 5 

147 50 

155 40 

159 59 
150 13 

154 44 

138 18 
163 44 

156 18 

159 56 

160 5 


a . 
11 

* s 

//. m. 

9 32 
9 44 
9 46 

9 54 

10 4 

10 4 
9 49 
9 40 
9 59 
9 51 

10 23 

9 32 

9 47 

9 40 

9 40 


11 

^ ^ Mes en que 

g. ,-j ocurre el 
o * fft máximo. 

Ill 

1.4 ! Febrero ...... 

1.2 Marzo 

1.5 Febrero 

1.4 Marzo 

1.3 1 Abril 

1.3 ! Marzo .. 

1.6 \ Febrero .... 

1. 5 t Noviembre . 

1.2 j Marzo. 

1.4 i ... Id - 

1.3 i Febrero 

1.2 Octubre ..- 


Mes en que 
ocurre el 
mínimo. 

Julio 

id....... 

id 

.-id--..- 

Junio 

Julio 

Junio. 

.id 

id 

Julio 

.....id 

id-.-...- 


Épocas extre- 
mas del máxi- 
mo de U2. 

Más 1 Más 
tem- 1 tarde, 
prano.! 


Mes de la 
época más 
temprana. 


Mes de la 

época más 

tardía. 


9 14 
9 27 


H. m. 
9 47 
9 57 


Noviembre _ 

id 

id 

(Noviembre . 
IDiciembre -- 

Noviembre , 

id 

Octubre 

id 

Noviembre . 
Octubre 

id 

NoviftTnhre 


Julio. 
Id. 


Bombay „_....,.. ... „. ., ,.. 

Rangroon ..,..„ , . .. . - 

Chittag'onft . » . -~ ... - _.....» 


9 15 
9 37 

9 39 


9 58 
10 12 

10 19 


Id. 

1 Agosto. 

1 Marzo. 
lAfirosto. 
Julio. 


Calcutta, . - . , . 


9 41 1 10 24 
9 28 10 10 
9 20 9 58 


K u rrac h ee » , ^. ... 

Aden * ..- *. , .«■ . - 


Junio. 
Abril. 


Allahabad....,-.... ...... ^..- 

Ma^pur ..,...-....,,....„......, ,.. ....... 

.Lahore .,.„,„...,._. ..-,.. 

Mauricio 

Sing-apore . ... ......_.-...w,.. 

Bata via.,,.,.. ., 

Maaila.,. ..^.^ 


9 37 

9 28 

9 56 
9 20 
9 31 
9 25 

9 25 


10 17 
10 10 

10 36 

9 45 


Agosto. 

Julio. 
/Febrero. 
IJulio. 

Febrero. 


^ o i/Marzo 

^•^ lAbril 

1.2 \ Octubre.... 

, , /Marzo ---... 
^•^ |lAbrü 


id 

}-- 


9 58 id 

9 54 id. 

o Ko iJ Noviembre _ 
^ ^2 IDiciembre - 


Julio. 

Febrero. 

IJulio. 
/Agosto. 



Por este cuadro puede apreciarse de un modo general la constancia en los caracteres 
principales de esta onda en las distintas regiones. Resalta sobre todo la concordancia 
casi perfecta de la relación entre el máximo y mínimo de la amplitud, y la época del año 
en que se verifica el mínimo. 

Nótanse sin embargo algunas discrepancias, que indican no estar del todo libre esta 
onda de la influencia de circunstancias locales. Así por ejemplo, en algunas regiones, 
la amplitud presenta dos máximos y dos mínimos en el año, mientras que en otras, Ma~ 
ni la entre ellas, desaparece esta doble variación. En estaciones que, como Batavia, Mau- 
ricio y Adden, presentan en Octubre o Noviembre el máximo absoluto, se verifica otro 
relativo en Marzo o Abril. 

La razón de estas variíiciones es bastante obvia. La onda segunda, tal como resul- 
ta del análisis, ha de ser una onda compleja, resultante de la interferencia de dos, con el 
mismo período de doce horas. Una de estas (que denominaremos onda semidiurna prin- 
npal) es la producida por las variaciones de temperatura en las capas superiores de la 
atmósfera, según lo anteriormente explicado; la otra {onda semidiurna secundaria) re- 
presenta las oscilaciones (lue en la presión ha de producir la componente semidiurna de 
la temperatura en las capas inferiores. La primera de estas dos presenta, según Hann 
y Angot, una doble variación durante el año, siendo máxima en los equinoccios, y míni- 
ma en los solsticios. La segunda, más irregular y de menor amplitud que la primera, 
por obedecer a influencias locales, ha de pasar por un solo máximo (en la estación seca) 
y un solo mínimo (en la lluviosa). 

Fácilmente se ve cómo de la interferencia de estas dos ondas puede resultar la 
desaparición completa, o casi completa, de la doble variación anual en la onda que estu- 
diamos; pues el equinoccio de otoño tiene lugar en la mayor parte de las regiones en la 
estación lluviosa, y el solsticio de invierno hacia el comienzo de la estación seca, de modo 
que el máximo y mínimo correspondiente a estas épocas en la onda principal, quedan con- 
trarrestadcM, o en todo o en gran parte, por coincidir con las épocas de amplitud míni- 
ma y máxima de la onda secundaria. Lo contrario sucede con el solsticio de verano y 
equinoccio de primavera, por coincidir el primero con la estación lluviosa, y el segundo 
con la seca. 

M. Angot, en la publicación ya citada, p. B311 y siguientes, propone un método bas- 
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tante ingenioso para la separación de estas dos ondas. No juzgamos oportuno detener- 
nos en la explicación de todo el procedimiento, que nos ocuparía demasiado espacio. 
Nos limitaremos a dar en la tabla de la pág. 25, los resultados que por este medio he- 
mos obtenido para Manila.^ 

No podemos comparar nuestros resultados con los de otras regiones, pues no sabe- 
mos que se haya verificado este análisis más que por el mismo Angot, y éste única- 
mente ha tomado uno o dos meses para las nueve regiones que somete a este cálculo; 
juzgamos sin embargo útil haber comprobado para Manila este carácter que se atri- 
buye a la onda semidiurna principal. 

Tercera onda. — Dos son las únicas hipótesis posibles acerca del origen de esta onda : 

1.^ Que provenga de una acción independiente de período de ocho horas. 

2.^ Que sea parte del efecto producido por las causas que originan la primera onda ; 
pues como indicábamos en la discusión de esta, las causas de la primera onda no pue- 
den representarse por un movimiento simplemente armónico. 

Lo primero que ocurre contra la primera hipótesis es, que no se conoce acción física 
independiente de período de ocho horas. A esto se añade que la amplitud de la tercera 
componente es muy pequeña, lo que dice mejor con el carácter de un residuo, que de un 
efecto total debido a una causa independiente. 

En confirmación de estas razones generales viene la propiedad que en todas las re- 
giones se nota en esta onda ; la variación anual de su amplitud, y de la épKXía de sus fases 
críticas sigue una ley parecida a la primera, con la particularidad de que el sentido de 
esta variación es casi perfectamente inverso, lo que hace sospechar que la tercera onda 
es una corrección de los valores de la primera. 

Dada la pequeña influencia que, tanto esta onda como la cuarta, han de tener en la 
variación barométrica, creemos no vale la pena detenerse más en su discusión. No 
queremos sin embargo omitir una particularidad de la tercera onda que hace notar el Dr. 
Hann; esto es, que hacia los equinoccios, la época de las fases experimenta una inver- 
sión notable, y que esta inversión concuerda exactamente con la época de los mínimos de 
la amplitud. Esta propiedad aparece bien marcada en Manila; la inversión de que ha- 
bla el Dr. Hann se verifica de Abril (0'^ 49"^) a Mayo (7^^ 15"^) y de Agosto (6*' 52"0 a 
Septiembre (1*^ 23"") ; los mínimos en la amplitud se verifican también de Abril a Mayo, 
y de Agosto a Septiembre, como puede verse en la Tabla V. 



^ La fórmula a que finalmente se llega por este método de Angot, es 

Qmm 926^11^'^^^^^^ (2m + 64"-2c) 

en la que h representa la presión media; 5 la declinación del sol; r el radio de la órbita terrestre; A 
la latitud del lugar; m las distintas horas expresadas en grados, de modo que será cero para las O"; 15' 
para la 1**; 30° para las 2'', etc.; e la ecuación del tiempo. Los valores que hemos tomado para estas 
cantidades son: 

/i = 759,10 X=:14° 35' 

2e ^^«; ^ 

r' 

Enero +4 48' 0.9033 

Febrero -f 7*^ 01' 0.9722 

Marzo +4" 19' 1.0043 

Abril 4-0° 02' 0.9597 

Mayo -r 46' 0.8734 

Junio +0" 08' 0.8197 

Julio -f 2" 41' 0.8424 

Agosto +1" 50' - 0.9208 

Septiembre -2" 32' 0.9852 

Octubre -^7" 00' 0.9806 

Noviembre -7^ 18' , 0.9193 

Diciembre ~~r 56' 0.8746 
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CONCLUSIÓN. 

De todo lo hasta aquí expuesto se deduce, que la teoría moderna de la doble osci- 
lación diaria del barómetro supone ser dos las causas del fenómeno. La una esta 
completamente libre de influencias locales, y es precisamente la que introduce en la curva 
barométrica esa doble oscilación, siempre tan misteriosa para los que han tratado de 
investigar su origen. Está representada por la que hemos denominado onda semidmrna 
prmoipal, y es debida a la variación de la temperatura en las capas superiores de la 
atmósfera. La otra varía según las circunstancias locales de cada región, por deberse 
a la variación de la temperatura en las capas inferiores. La curva que la representa 
suele denominarse onda térmica. En la tabla de la página 26 damos sus valores para 
todos los meses; están deducidos restando de la variación observada (Tabla III) la onda 
semidiurna principal. 

En las Láminas XIII y XIV pueden verse las curvas que representan la onda 
térmica, y la onda semidiurna principal, para varios meses. 

Estas son sin duda las dos ondas de que habla M. Angot en su "Traite élémentaire 
de Météorologie'' (p. 104) al explicar la variación diurna del barómetro. 

Sin saber de qué ondas se trata no es fácil darse perfecta cuenta de esta expli- 
cación. Hemos juzgado pues será de interés traducirla aquí, pues al mismo tiempo que 
se entenderá mejor teniendo a la vista las curvas, servirá como discusión de estas. 

'*La variación diurna del barómetro no está aún explicada de una manera completa. 
Las investigaciones más recientes parecen indicar que resulta de la superposición de dos 
oscilaciones distintas; una, que presenta dos máximos y dos mínimos por día, cuyo pe- 
ríodo es, por consiguiente, de medio día, y que se denomina, por esta razón, oscilación 
n. onda semidiurna; otra, que presenta un solo máximo y un solo mínimo en veinticuatro 
horas, y que se denomina onda diurna. 

**Los dos máximos diarios de la onda semidiurna serán iguales entre sí, lo mismo 
que los dos mínimos. En toda estación del año, y para toda la tierra, los dos máximos se 
verificarán hacia las 10*' de la mañana y de la tarde; los dos mínimos hacia las 4^ de 
la mañana y de la tarde. La diferencia entre los máximos y los mínimos de la onda 
semidiurna — es decir, la amplitud de esta onda,' varía poco, para una misma estación, 
durante el año; es máxima en los equinoccios, y mínima en los solsticios. Esta ampli- 
tud no dependerá de condiciones topográficas, sino únicamente de la latitud — cerca de 
2 mm. en el ecuador, disminuye primero lentamente hasta los trópicos, luego muy rá- 
pidamente, de modo que llega a ser muy pequeña en latitudes altas. Es probable que 
esta onda semidiurna esté causada por la acción del calor solar en toda la masa de la 
atmósfera; pero no se ha podido explicar todavía convenientemente el modo como se 
produce. 

^'Mientras que la onda semidiurna es la misma para todas las regiones situadas 
en la misma latitud, la onda diurna tiene una amplitud muy variable, que depende de 
las condiciones topográficas de la región, y que puede ser muy diferente aún para esta- 
taciones vecinas. En general es mucho más débil en estaciones costeras que en estaciones 
continentales, y su amplitud está en razón directa de la amplitud de la variación diurna 
de la temperatura. Las horas del máximo y del mínimo varían notablemente según 
los países o las estaciones del año; el máximo se verifica más temprano en verano que 
en invierno; su época media se acerca a las 9^'; el mínimo, al contrario, cuya época 
media es cerca de las 1(>\ se produce más temprano en invierno que en verano. El 
período de bajada dura pues, por término medio solamente siete horas, mientras que 



* M. AiiíTot toma aquí como amplitud de la onda, la amplitud de una oscilación completa; otros 
autores, y el mi.^mo Anj^ot en la publicación que hemos citado varias veces, toma como amplitud de la 
onda la amplitud de una semioscilación, o la ordenada máxima de la curva. 
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el período de subida se extiende a diez y siete horas ; la subida no es uniforme durante 
este tiempo; rápida inmediatamente después del mínimo, se detiene durante la noche, y 
crece de nuevo con más rapidez entre la salida del sol y la hora del máximo. 

"Mientras que la teoría de la onda semidiurna no está hallada aún ' la onda diurna 
se explica, por el contrario, completamente por el efecto local de la variación diurna 
de la temperatura. Durante la noche, a medida que se produce el enfriamiento en las 
capas inferiores de la atmósfera, estas se contraen, y para llenar el vacío que pro- 
duciría esta contracción, el aire afluye por todos lados, sobre todo por la parte superior ; 
la masa de aire situada sobre la región que se considera se hace por lo tanto mayor, 
y la presión aumenta. Desde la salida del sol las capas inferiores se calientan rápi- 
damente y tienden a elevarse; pero la calefacción se produce primero en las capas 
que están en contacto con la superficie terrestre; las capas de aire, comprendidas 
entre la superficie y las capas superiores que no se calientan aún, no aumentan al principio 
de volumen de una manera apreciable; lo mismo que sucede cuando se calienta el 
aire en una cámara cerrada, la elevación de temperatura se manifiesta entonces prin- 
cipalmente por un aumento de presión; el primer efecto de la calefacción diurna es por 
consiguiente un aumento de presión en las capas que están en contacto con el suelo ; el 
barómetro comienza pues a subir. 

''Un poco más tarde entre 7^ y 11^ de la mañana, según los países y las estaciones 
del año, la calefacción de las capas inferiores es suficiente para determinar por fin un 
movimiento ascendente general del aire sobre la estación; el aire se dilata; pero desde 
que la columna de aire situada sobre la región caliente alcanza mayor altura que la 
que cubre las regiones circunvecinas, el aire se extiende hacia los lados, desde la región 
más caliente a los meridianos más distantes del sol ; la cantidad de aire que insiste sobre 
la estación disminuye por lo tanto, y el barómetro baja. Hacia el fin del día, cuando 
la temperatura comienza a disminuir rápidamente en las capas inferiores, la corriente 
ascendente pierde actividad, cesa y acaba por ser reemplazada durante la noche por un 
movimiento inverso; el barómetro sube entonces, rápidamente primero mientras la tem- 
peratura disminuye con rapidez; después más lentamente hasta la salida del sol.'' 

Nótese que la explicación que de la onda térmica o diurna da M. Angot coincide 
en gran parte con la primera teoría de la variación diurna del barómetro propuesta por 
Raimond, y perfeccionada luego por Kreil y Renou, y que era insuficiente cuando se la 
aplicaba a la variación diurna total de la presión. 

Concluyamos pues, que el análisis de las curvas que representan los diversos ele- 
mentos meteorológicos de cada región, aunque un tanto largo y penoso, puede ser de 
gran utilidad, tanto para el estudio simultáneo de varios fenómenos, cuyas relaciones 
mutuas queramos investigar, como para profundizar más en las causas que en cada uno 
de ellos puedan intervenir. 

* No sabemos si M. Angot tendría noticia de las investigaciones del Dr. Margules, pues ni aquí 
ni en el trabajo que publicó en los Anales del Buró Central de Francia, hace mención de ellas, 
siendo de tan capital importancia en la explicación de esta onda. El Dr. Hann en publicaciones pos- 
teriores a las de Angot, se lamenta de que los trabajos del Dr. Margules no sean más conocidos, y de 
que no hubieran llamado aún la atención de la **Meteorologische Zeitschrift." 



DEDUCCIÓN DE LAS ECUACIONES DE BESSEL, PARA LA 

DETERMINACIÓN DE LOS COEFICIENTES 

DE LA FÓRMULA (I). 



Supongamos que se nos da un número cualquiera n de observaciones. 

Representemos por a^, a^, a^, etc., las observaciones correspondientes a las horas O, 
1, 2, etc. 

Vimos en la página 43 que estas observaciones están representadas por las ecua- 
ciones siguientes:^ 

ao=p-\-pi+p2+p3 

ai=p-\-pi eos x + Qi sen x+ps eos 2x + q2 sen 2x+ps eos 3x-f ^3 sen 3x + 

a>=:p+pi eos 2x + qi sen 2x-\-p2 eos 4x + g2 sen 4x-f pa eos ^x + q^ sen 6x + 

az=.p+px eos 3a:-f gi sen Zx-^p-z eos Gx + g^ sen Qx+p^ eos Qx + g» sen 9x + 

an— i=:p+pi eos (n— l)x+gi sen {n — l)x+p2 eos 2(n— l)x+g2 sen 2{n—l)x-\-p'i eos 3(n— l)a^ 

+^3 sen 3(n— l)x + 



Si el número de los coeficientes que pretendemos determinar, fuera el mismo que 
el de nuestras observaciones, tendríamos un sistema de tantas ecuaciones como incóg- 
nitas, y podríamos resolverle fácilmente por los métodos ordinarios del álgebra. Pero 
ordinariamente el número de las observaciones excede con mucho al número de coefi- 
cientes que necesitamos ; así, en el análisis de la curva barométrica, teníamos veinticuatro 
observaciones, y empleábamos únicamente ocho de estos coeficientes. 

En este caso el sistema es más que determinado; y en la imposibilidad de hallar 
los valores de p y g que, sustituidos en las ecuaciones, reproduzcan exactamente nues- 
tras observaciones, hemos de contentarnos con averiguar cuáles son los valores más pro- 
bables dé ?} y (?, que dan para los segundos miembros valores compatibles con las 
observaciones. 

En la teoría de los mínimos cuadrados se prueba, que estos valores más probables 
son aquellos que hacen mínima la suma de los cuadrados de los errores residuos. 

Error residuo es lo que se aparta del valor observado el valor obtenido al sustituir 
en las ecuaciones los valores más probables de las incógnitas. 

Según esta definición, designando por v,,, v,, v^, etc., los errores residuos, y supo- 
niendo que p, Pj, q^, etc., son los valores más probables de estos coeficientes, podemos 
escribir : 

Vo= — CN, + p + pi+ P2 + P8 + 

Vizr— tti+p+pi eos x + qt sen x+p^ eos 2X'i-q.i sen 2x-\-p,. eos Sx-fga sen 3x + 

V2=--(i2+p+pt eos 2x+g, sen 2x+p2 eos áx-j-q, sen ix-\-pn cos6x-{-q3 sen 6ír + 

Vgrn— Oa+p-fpi COS Sx+qi sen 3x-fp2 eos 6x+qi¿ sen 6x4- P» eos 9x+qt sen 9x + 



i;n— i=an— i+p+pi eos {n—l)x + qi sen {n — l)x+p¿ eos 2(n — l)x + g2 sen 2(n~l)x 

+p^ eos 3(n— l)x+g3 sen 3(n — l)x+,...:. 

La condición que hace más probables los valores de p y g es 

vlA'VÍ+vl+vl+ + v^n i=zmínimo. 

y escribiendo en vez de t;|, vj, etc., sus valores, queda 

(~-a,+p+p,+?>.,+ )2+ (-™a,+p+p, eos x+q, sen x+ .)^ 

+ (— ag+p-fPi eos 2x-j-q^ sen 2x-^ .)2=mínimo. 

Ahora bien ; sabemos por el cálculo diferencial que, cuando una función está en un 
mínimo, sus derivadas con respecto a cada una de las variables son cero. Por consi- 

^ La p, que aparece en estas ecuaciones, representa la constante M que vimos completaba la fór- 
mula de Bessel; nótese además que las ecuaciones de la página 43 están dadas aquí según la Fórmula (I) . 
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guíente, derivando la expresión anterior con respecto a p, p„ g„ etc., y representando 
por P, P,, Q„ etc., estas derivadas, obtenemos el sistema siguiente: 

P=2(— .ao4-/?4»p,+p2+ _)-\-2{^at+p+pi eos x^^qt sen x^- ...) 

+2(— Os+P+Pi eos 2a;+-g, sen 2a; + )+ —O 

Piz=:2(-~ao+p+pi+pn +...„ .)+2 eos x{—(h+p+pi eos x + qi sen x + ) 

+2 eos 2a:(— Oa+P+Pi eos 2x + q, sen 2x+ )+ = 

Qi=2 sen x(—an+p-\-pt eos o^ + g* sen x-\- ) 

+ 2 sen 2x(-^€h-^p+pt eos 2a; + gi sen 2x + )4_ z^o 

F2=2(— a«+p+pt+p2 4- )+2cos 2a;(-»ai+P+pi eos x + qt sen ír + ) 

+2 eos 4ír(~-a2+p+pi eos 2x + qt sen 2x+ )+ — O 

Q = 2 sen 2a;( — ai+p+p, eos x + qi sen a:+ ) 

+2 sen 4a7(--a2+p + Pi eos 2x + qt sen 2a; + )+ =0 

Este es el sistema que hemos de resolver para hallar los valores de p y g, que serán 
los más probables; pues en esta hipótesis se han obtenido. las ecuaciones precedentes. 
Después de dividir por 2, podemos darle la siguiente forma: 



iP=~ao + np -\~ pi 
— a, 



— ttn-l 

— at eos x 

— a2 eos 2x 

— as eos Zx 



+ q^ 



+ eos X 
+ eos 2x 
+ eos 3x 
4- eos 4x 



+ eos(íi-l)a; 



sen X 
+ sen 2x 
+ sen 3x 
4- sen áx 



+ P2 



-fP 



~an-icos(ti~l)a^; 



-f eos a; 
-f eos 2x 
-f eos 3¿f 



+ sen(n~-l)í;c 

+ eos' ar 
+ eos" 2x 
+ eos^ 3x 



+ eos 2x 
+ eos 4x 
+ eos 6x 
+ eos 8x 



+ ^2 



+ cos(n~l)a; 



+ P2 



+ eos 2(n~l)a; 



sen 2x 
+ sen 4»; 
+ sen 6x 
+ sen 8a? 



+ - 



+ cos"(n-l)a- 



+ sen 2(n-l)5c 
eos X sen x 
+ eos 2x sen 2a; 
+ eos 3a: sen Sx 



+g.- 



+ eos X eos 2a: 

H~ eos 2a: eos 4x 
-f eos 3a- eos 6x 



iQi-- 



at sen a: 
a 2 sen 2a: 
«8 sen 3a^ 



4-p 



— íiu-t sen(ti-l)a; 



sen X eos a: 
-f sen 2a: eos 2a: 
-|r sen 3a: eos Bx 



í 4-cos(n~l)a: eos 2(n-l)a: 

sen a* -fp 

+ sen 2a: 
-f sen 3a: 



+ sen(n~l)x' 



+ eos(ti-l)a: sen(n-l)a: 
eos X sen 2a: 
+ eos 2a: sen 4x 
+ eos 3a: sen 6x 



+ ■ 



-f eos ( n-1 ) X sen 2 ( ti-1 ) x 



+gi 



I sen" X 
I + sen^ 2x 
! 4- sen* 3a: 



-f sen(n-l)xeos(ti-l)x i + sen'(w~l)x 

-j-Pi I sen X eos 2a: +g« i sen x sen 2a: 

+ sen 2a: eos 4a: \ + sen 2x sen 4a: 

+ sen dx eos 6a' I + sen 3a: sen 6x 



+ - 



-rzO 



= 



-f sen (w~l ) X eos 2 (n-1) íe 



+ sen {n-1 ) x sen 2 (n-1 ) x 



(A) 
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Este sistema según la notación explicada en la página 47, puede escribirse del 
modo siguiente: 

— S am+np+pi Seos mx+Qi S sen mx+p¿ Seos 2mx+q2 S sen 27nx + = 

— S am eos mx+p S eos mx+pi S eos* mx + qi S eos mx sen mx+pi S eos mx eos 2mx 

+ ^2 S eos mx sen 2mx + = 

—2 am sen mx+p 2 sen mx+pi 2 sen wa; eos ma? + ^i 2 sen^ mx+pi 2 sen twa; eos 2mx 

+^3 2 sen m^ sen 2ma;+ =0 

—-2 am cos 2mx+p 2 eos 2m^-\-px 2 eos 2ma; eos m^ + qi 2 eos 2mx sen mx+pt 2 eos'* 2ma; 

+ g3 2 eos 2ma; sen 2mx+ .....=0 

--2 am sen 2mx+p 2 sen 2mx+pi 2 sen 2mx eos mít?+gi2 sen 2nix sen wx+/?2 2 sen 2fnx eos 2ma; 

+^2 2 sen^ 2mx+ =0 

Cuando nx (producto del número de observaciones por el número de grados de las 
partes en que están divididos los distintos círculos de referencia) es múltiplo de 27r, el 
sistema anterior se reduce al siguiente: ^ 

— S a„, eos mx+Pi S cos^ mx=0¡ 
—2 ürn sen mx-{~q^ 2 sen^ mx=o{ ,-ri\ 
—2 an, cos 2mx-\-P2 S cos^ 2mx=ol 
— 2 a„, sen 2ma;+g2 2 sen^ 2mx=o] 



Pero como luego probaremos 



2 cos^ mx~^^ 



2 sen^ mx=J^^ 



2 sen^ 2wíi;=| 



luego, sustituyendo estos valores en el sistema (B) se obtiene finalmente 

Íp_S a^ 
n 
^~Pi=2 a™ eos mx 
. ^9i=2 am sen ma; 
—2 a,„ eos 2mX'{- ^P2=0l ^^2=^ ^m eos 2mx 
—2 a« sen 2mx+ p Í2=0 I p 0^2=2 a», sen 2ma; 



Estas son, en su forma general, las ecuaciones de Bessel que buscábamos ; en nuestro 
caso particular de veinticuatro observaciones, ^=12, lo que reduce estas ecuaciones a 
las que escribimos en la página 47. Queda probado también que la cantidad cons- 



ignemos juzgado conveniente dejar para el fin esta demonstr ación, para no interrumpir el 
raciocinio. 

127182 6 
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tante p o M, que completa la fórmula de Bessel es, como indicamos antes, la media 
aritmética de las observaciones; así en nuestro caso, tenemos: 



P= 



24 



Fáltanos únicamente probar el paso del sistema (A) al sistema (B) y las igual- 
dades (C). 

Para lo primero, nos basta probar que, cuando nx es múltiplo de 27r, 5sen mx=0, 
2 eos ma>=0, 2 sen 2ma;==0, 2 eos 2mx=0, etc. 

Tomemos la fórmula bien conocida en trigonometría: 

2 sen z sen y=coB (z-^y) -cob (z+y) /.sen z= ^^^ iz~-y)^coB (^+3/) ^) 

Por una serie de aplicaciones de la fórmula (1) en la que escribiremos sucesiva- 
mente en vez de z, x, 2x, 3x, , (n-l)x, y en vez ie í/,| obtenemos: 

X Sx 

cos^-cos--- 

sen x=^ — = =- 

2 sen I 

3a; 5a! 
cos^l^ - eos- 
sen 2a;= ^ 



2 sen I 

, ¿i 

hx Ix 

cos-^-cos-g- 

sen Zx— = 

2 sen I 

¿i 



sen (n~2)a;= 

sen (w— 1)0?= 
y sumando ordenadamente queda : 



eos -^ — ;r — — — eos -^i ~ — — 



2sen| 

eos -^— -•■---- — COS-^^^ — '--^ 



2sen^ 



cos^ — eos "^ ' 



o o 

sen a^+sen 2íi:+sen 3a?+ +sen (ti~l)a;=2 sen wa;= — = ^ 



2 sen^ 



cos^— eos nx cos|— sen nx sen| 

di £á éd 



2sen| 
si nx es múltiplo de 27r, sabemos que 

eos na:=l ; sen naj=0 

y queda finalmente : 

X x 

COSg- cos^ 

2 sen mx = - — ™^_=o 
2 sen? 
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Por análogo procedimiento se probaría que 2 sen 2mx=0, 5 sen Zmx=0, etc., así para 

este último, escribiríamos en la fórmula (1) Sx, 6a;, 9íc, , {n~l)Sx, en vez de 

3a? 
z,y~ en vez de y, lo que daría : 

cos^^cos (^^-"3)^ cos|~~cos Snx cos?|~sen Snx sen^ 

2 sen Smx= ? ^— ? = ^ 2 2^ 

2 sen^ 2 sen^ 

y como Snx es múltiplo de 27r, resulta: S sen 3mx=0 

Para probar que 2 eos ma;=0, 2 eos 2ma;=0, etc., nos serviremos de la fórmula trigo- 
nométrica 

2 eos z sen 2/=sen (2:+|/)-sen (;2-2/)/.cos g^^-gLJ^+^)""Q^^ (■gz:l)(2) 

2 sen y 

y escribiendo en vez de z, O, x, 2x, , (n-l)a:, y en vez de y, |, obtendremos por 

análogo procedimiento : 

sen|+sen| 
cosO= — ^ — ^ 

3a; X 

sen "^ — sen-^ 

2 2 

cosa:= =~— — -^ 

2 sen| 

sen ^ — sen ^ 

eos 2x= — ' 

2 sen| 

sen'?! — sen^ 
eos 3a: = 



2 sen| 



eos (n-2)a;= 



o^^ (2w^3) X ^^„ (2n^5) a? 
sen-^^ — ~ — ^— — sen -^= — ~ — '— 

Zi ^ 

2 sen| 
sen-^^ — ^ — sen-^ — ^ 



cos()t-l)a;=- 



Sumando ordenadamente obtenemos : 



2 sen? 



8en?+8en-<2»^-l>^ 



2 — 2 
1+cos a;+cos 2a;+cos Bx-{-.: +cbs (w— l)£c=5 eos mx=— 



2 sen? 



sen^+sen nx eos^— eos nx sen| 
2 sen ^ 
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y teniendo en cuenta que nx=%r, queda: 

X X 

sen ^ — sen ^ 

% eos mx:=z = ===:0 

2seng 

Probado lo que precede, fácilmente se ve, que los demás términos del sistema (A) 
de la forma 

2 eos mx sen mx ; 2 sen mx eos 2mx ; 2 sen mx sen 2mx ; 2 eos mx sen 2mx ; 

2 eos mx eos 2nix 

se reducen también a cero; pues recordando que el producto del seno por el coseno de 
un mismo arco es igual a la mitad del seno del arco duplo, se obtiene: 

2 eos mx sen ma;=i2 sen 2mx=0 

lo mismo podemos probar para las cuatro expresiones restantes valiéndonos de las 
bien conocidas fórmulas de trigonometría : 

sen X eos t/=|sen (x-{-y) +|sen (x—y) 
eos X sen i/=isen (x-^y) — Isen (x—y) 
eos X eos y=^coB (a;+2/)+icos (x—y) 
sen X sen i/==|cos (x— ?/)— -feos (a:-fí/) 

pues escribiendo en ellas en vez de x e y, mx y 2mx, obtenemos : 

2 sen mx eos 2ma;=¿2 sen 3ma;— ^2 sen mx . 

2 eos mx sen 2mx=|2 sen 3ma:+|2 sen mx 
2 eos mo? eos 2mx=^% eos 3ma7+i2 eos mo? 
2 sen mx sen 2ma;=i2 eos mx—^t eos 3ma; 

y por lo probado anteriormente, los segundos miembros son cero. 

Pasemos a probar las ecuaciones (C) ; un fácil desarrollo trigonométrico hace ver que: 

2 cos^ a?=l+cos 2x 
según esto tenemos: 

2 cos^ O==l+cos O 
2cos^ ír==l+cos 2x 
2 cos^ 2a?=l+cos áx 
2 cos^ 3a:=l+cos Bx 



2 cos^ (n--^ 1)0^=1+ eos 2(n— l)aj 

y sumando ordenadamente 

2 2 cos^ ma?=n+l+cos 2x+cob 4a;-f-cos 6a?+ -feos 2(ti-"l)a;=n+2 eos 2mx. 

pero según lo anteriormente demonstrado, 

2 eos 2mx==^0 
luego 

2 cos^ ma:=| 

Por análogo procedimiento probaríamos este mismo valor para un término cualquiera 
de la forma 2 cos'i mkx. 

Para probar que 2 sen^ mx-^% nc^ valdremos de la conocida igualdad» 
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y escribiremos 

sen^ 0=1 —eos- O 

sen^ a:=l— eos" x 

sen^ 2a:=:l— eos2 2a; 

sen^ 3ír=l~-eos2 3a? 



sen^ (n™l)a;=:l— cos^ (n~-l)x 
sumando ordenadamente tendremos : 

S sen^ mx=n—^ eos"' mx 

pero acabamos de probar que 2 eos- mx='~^ luego: 

2 sen^ mx^^n—"^ = ^ 
2 2 

No juzgamos oportuno detenernos más en dar a conocer todo el método de Bessel; 
quien en ello estuviera interesado, podría consultar la explicación original del mismo Bes- 
sel, traducida al inglés por Mr. R. H. Scott en un apéndice al "Quarterly Weather Report 
of the Meteorological Office, London, 1870.'' 

O 



